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ABSTRACT 


Adaptive  nulling  is  applied  to  the  problem  of  generating  a  therapeutic  thermal 
dose  distribution  in  electromagnetic  hyperthermia  treatment  of  cancer.  A  system 
design  concept  for  implementing  adaptive  hyperthermia  is  introduced.  With  the 
proposed  design  concept,  it  may  be  possible  to  maximize  the  applied  electric  field 
at  a  tumor  position  in  the  target  body  and  simultaneously  minimize  or  reduce 
the  electric  field  at  target  positions  where  undesired  high  temperature  regions  (hot 
spots)  occur.  In  a  clinical  situation,  either  a  gradient  search  algorithm  or  sample 
matrix  inversion  algorithm  would  be  used  to  rapidly  form  the  adaptive  null  (or  nulls) 
prior  to  any  significant  tissue  heating.  Auxiliary  short-uipole  field  probes  are  used  in 
effecting  the  iesired  electric  field  nulls.  An  adaptive  null  formed  at  the  surface  of  the 
target  has  a  finite  width  and  extends  into  the  target.  This  finite  null  width  can  allow 
for  noninvasively  positioned  auxiliary  probes  in  an  adaptive  hyperthermia  system. 
The  allowed  spacing  or  resolution  between  a  deep  null  and  focus  is  fundamentally 
equal  to  the  hyperthermia  antenna  half-power  beamwidth  in  the  tissue.  A  closer 
spacing  between  the  desired  null  position  and  focus  is  achieved  by  reducing  the  null 
depth  via  the  signal-to-noise  ratio  at  the  auxiliary  probe. 


Analysis  of  an  annular  phased-arrav  antenna  embedded  in  an  infinite  homo¬ 
geneous  medium  shows  the  potential  merit  of  combining  adaptive  nulling  with  con¬ 
ventional  near-field  focusing  used  in  hyperthermia.  The  analysis  is  based  on  a  well- 
known  moment-method  theory  for  conducting  thin- wire  antennas  in  a  homogeneous 
conducting  medium.  The  theory  and  software  used  to  compute  the  moment-method 
received  voltage  at  a  short-dipole  probe  due  to  a  transmitting  dipole  array  are  doc¬ 
umented.  Computer  simulation  results  are  presented  for  a  fully  adaptive  eight- 
element  annular  array  operating  at  120  MHz  and  embedded  in  homogeneous  muscle 
tissue.  Multiple  simultaneous  nulls  are  used  in  adaptively  generating  a  desired  ra¬ 
diation  pattern  in  the  muscle  tissue.  The  power  received  at  a  movable  short-dipole 
E-field  sensor  is  then  used  as  the  power  source  in  a  second  computer  simulation 
that  performs  a  transient  thermal  analysis  of  an  elliptical  target  body  surrounded 
by  a  constant-temperature  water  bolus.  The  computer  simulations  show  that  adap¬ 
tive  nulling  can  prevent  undesired  high-temperature  regions  from  occurring  while 


simultaneously  heating  a  deep-seated  tumor  site. 
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1.  INTRODUCTION 


The  successful  treatment  of  deep-seated  malignant  tumors  within  a  patient  is  often  a  difficult 
task.  The  objective  of  the  treatment  is  to  reduce  in  size  or  completely  remove  the  tumor  mass 
by  one  or  more  modalities  available  at  the  treatment  facility.  Common  modalities  are  surgery, 
chemotherapy,  and  x-ray  therapy  [1].  A  modality  used  alone  or  in  conjunction  with  one  of  the  above 
modalities  is  “tissue  heating,”  or  hyperthermia  [1,2, 3, 4],  Hyperthermia  can  be  considered  as  a  form 
of  high  fever  within  the  body;  a  controlled  thermal  dose  distribution  is  required  for  hyperthermia 
to  have  a  therapeutic  value.  Typical  localized-hyperthermia  temperatures  required  for  therapeutic 
treatment  of  cancer  are  in  the  43-45°C  range.  Normal  tissue  should  be  kept  at  temperatures 
below  43° C  during  the  treatment.  The  most  difficult  aspect  of  implementing  hyperthermia,  with 
either  radio- frequency  (RF)  waves  or  acoustic  (ultrasound)  waves,  is  producing  sufficient  heating  at 
depth.  Multiple-applicator  RF  hyperthermia  arrays  are  commonly  used  to  provide  a  focused  near¬ 
field  main  beam  at  the  tumor  position.  A  focal  region  should  be  concentrated  at  the  tumor  with 
minimal  energy  delivered  to  surrounding  normal  tissue.  As  the  hyperthermia  antenna  beamwidth 
is  proportional  to  the  wavelength,  a  small  focal  region  suggests  that  the  RF  wavelength  be  as  small 
as  possible.  However,  due  to  propagation  losses  in  tissue,  the  RF  depth  of  penetration  decreases 
with  increasing  transmit  frequency.  One  of  the  major  problems  in  heating  a  deep-seated  tumor 
with  a  hyperthermia  antenna  is  the  formation  of  undesired  “hot  spots”  in  surrounding  tissue. 
This  additional  undesired  heating  often  produces  pain,  burns,  and  blistering  in  the  patient,  which 
requires  terminating  the  treatment  immediately.  The  patient  does  not  receive  anesthetics  during 
the  hyperthermia  treatment  in  order  to  provide  direct  verbal  feedback  of  any  pain.  Thus,  techniques 
for  reducing  hot  spots  are  necessary  in  hyperthc.mia  treatment. 

Several  journals  have  published  special  issues  on  acoustic  and  electromagnetic  hyperthermia 
treatment  of  cancer  [5,6].  Electromagnetic  transmitting  arrays  in  the  frequency  band  60-2000  MHz 
are  used  to  localize  heating  of  malignant  tumors  within  a  target  body.  Many  studies  have  been 
conducted  to  produce  improved  therapeutic  field  distributions  with  hyperthermia  phased  arrays  [7- 
20].  Phase  control  can  be  used  to  synthesize  therapeutic  RF  radiation  patterns  without  adaptive 
control  of  the  transmit  array  weights  [7-12].  Array  transmit  weights  can  be  adaptively  controlled 
to  maximize  the  tumor  temperature  (or  RF  power  delivered  to  the  tumor)  while  minimizing  the 
surrounding  tissue  temperature  (or  RF  power  delivered  to  the  surrounding  tissue)  [13-20].  All  other 
studies  require  invasive  techniques  to  optimize  the  radiation  pattern  [13-20].  This  report  addresses 
the  potential  benefit  of  using  adaptive  nulling  with  noninvasive  auxiliary  sensors  to  reduce  the  field 
intensity  at  selected  positions  in  the  target  body  while  maintaining  a  desired  focus  at  a  tumor. 

It  is  useful  to  describe  certain  features  of  a  clinical  hyperthermia  system  used  to  treat  deep- 
seated  tumors.  A  hyperthermia  annular  phased-array  antenna  system  (Model  BSD-2000,  SIGMA- 
60  applicator  [21],  BSD  Medical  Corporation,  Salt  Lake  City,  Utah)  is  shown  in  Figure  1.  By 
fully  surrounding  the  patient  with  an  annular  phased  array,  it  is  possible  to  obtain  constructive 
interference  (or  signal  enhancement)  deep  within  the  target  volume.  This  hyperthermia  system  uses 
a  60-cm  array  diameter  with  eight  uniformly  spaced  dipole  elements  operating  over  the  frequency 
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Figure  1.  A  clinical  hyperthermia  phased- array  antenna  system. 
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band  60-120  MHz  [22].  The  eight  dipoles  are  fed  as  four  active  pairs  of  elements.  There  are  four 
high-power  amplifiers  which  drive  the  dipole  pairs  with  up  to  500  W  average  power  per  channel. 
Each  of  the  four  active  channels  has  an  electronically  controlled  variable-phase  shifter  for  focusing 
the  array.  Temperature  and  electric-field  probe  sensors  (both  invasive  and  noninvasive)  are  used 
to  monitor  the  treatment.  A  cool-water  (5-40°C)  bolus  between  the  patient  and  the  phased  array 
is  used  to  prevent  excess  heating  of  the  skin  surface.  The  bolus  is  filled  with  circulating  distilled 
water,  which  has  a  very  low  propagation  loss.  Computer  simulation  models  for  this  hyperthermia 
antenna  system  with  inhomogeneous  targets  have  been  developed  [7,9].  To  the  author’s  knowledge, 
the  hyperthermia  system  in  Figure  1  has  not  been  used  to  form  adaptive  nulls,  although  the  system 
may  have  sufficient  temperature  and  electric-field  sensing  capability  to  implement  adaptive  nulling 
with  an  appropriate  transmit-weight  control  algorithm.  A  candidate  adaptive  nulling  algorithm 
would  be  a  gradient  search  based  on  minimizing  the  signal  received  by  electric-field  sensors  at  the 
desired  null  positions  while  maximizing  the  signal  at  the  focal  point.  The  reason  for  considering 
this  algorithm  is  that  the  BSD-2000  system  is  currently  capable  of  measuring  the  applied  electric- 
field  amplitude  but  not  the  phase.  Note:  the  E-field  detectors  used  in  the  BSD-2000  system  are 
semiconductor  diode  detectors.  A  potentially  better  approach  to  performing  adaptive  hyperthermia 
is  to  use  a  channel-covariance-matrix-based  algorithm  to  control  the  weights.  If  the  electric-field 
phase  and  amplitude  are  measured  (by  modification  of  the  E-field  probe  hardware  to  form  in-phase 
and  quadrature  signals)  the  sample  matrix  inversion  (SMI)  algorithm  [23]  could  be  used  to  form 
the  nulls.  The  concept  of  adaptive  nulling  as  it  applies  to  hyperthermia  will  now  be  described. 

Recently,  the  author  has  been  investigating  a  technique  called  focused  near-field  adaptive 
nulling  [24-30] .  for  which  the  general  concept  is  shown  in  Figure  2.  A  calibration  probe  antenna 
is  used  to  focus  the  main  beam  at  approximately  one  diameter  D  of  the  phased  array.  Additional 
(auxiliary)  probe  antennas  are  located  in  the  sidelobe  region  of  the  phased-array  quiescent  radiation 
pattern.  Adaptive  nulls  are  formed  in  the  direction  of  these  auxiliary  probes.  Computer  simulations 
show  that  a  radar  system  can  be  tested  by  using  the  focused  near-field  nulling  technique  [29].  Ex¬ 
perimental  measurements  of  focused  near-field  adaptive  nulling  are  available  [30].  During  the  radar 
system  testing  study,  similarities  between  the  radar  application  and  electromagnetic  hyperthermia 
were  noticed  by  the  author.  These  similarities  are  that  1)  phased-array  near-field  focusing  is  used 
to  form  the  main  beam  and  2)  uncontrolled  or  high  sidelobes  are  potentially  deleterious  to  system 
performance.  For  the  radar  system,  interference  signals  entering  through  uncontrolled  sidelobes  can 
reduce  the  signal-to-noise  ratio  (SNR);  adaptive  nulling  [31]  is  commonly  used  to  counteract  this 
possible  degradation.  In  the  hyperthermia  system,  high-transmit  antenna  sidelobes  can  give  rise 
to  hot  spots  in  the  target  tissue.  These  high-temperature  regions  could  possibly  be  alleviated  with 
adaptive  nulling  of  the  electric-field  sidelobes  in  the  target.  Much  of  the  software  developed  for  the 
radar  system  testing  analysis  may  be  readily  modified  to  handle  adaptive  hyperthermia  analysis. 

There  is  continued  interest  at  Lincoln  Laboratory  in  applying  techniques  used  in  Department 
of  Defense  (DoD)-sponsored  research  to  other  disciplines,  such  as  medical  applications.  This  study 
applies  radar  system  testing  technology  (currently  in  the  form  of  software  simulation)  to  the  medical 
application  of  hyperthermia.  Figure  3  depicts  the  similarities  between  radar  system  near-field 
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Figure  2.  General  concept  for  focused  near-field  adaptive  nulling. 


testing  and  electromagnetic  hyperthermia  treatment  of  cancer.  In  Figure  3(a),  a  radar  phased-arrav 
antenna  system  is  positioned  in  an  anechoic  test  chamber  with  radiating  sources  that  simulate 
typical  radar  signals.  Notice  that  the  test  antenna  is  focused  in  the  near  field  at  one  aperture 
diameter.  Computer  simulations  have  shown  that  the  near-held  signals  (jammer,  radar  clutter, 
and  target)  received  and  processed  by  an  adaptive  radar  system  are  equivalent  to  fielded  radar 
conditions  [29’.  In  Figure  3(b).  a  hyperthermia  (transmit)  phased  array  is  used  to  therapeutically 
illuminate  a  target  body  by  focusing  its  electromagnetic  energy  on  a  malignant  tumor  deep  within 
the  body.  Nulls  are  formed  adaptively  to  reduce  the  electromagnetic  energy  delivered  to  potential 
hot  spots.  The  nulls  are  shown  invasive  to  the  target;  however,  noninvasive  null  positions  on  the 
surface  of  the  target  are  also  possible  as  shown  in  this  report. 

An  example  of  a  deep-seated  tumor  is  cancer  of  the  prostate  [1].  The  tumor  volume  often 
has  a  decreased  blood  flow  which  aids  in  heating  the  tumor,  compared  to  normal  tissue  for  which 
heat  is  carried  away  by  normal  blood  flow.  Note:  a  blood-flow  thermal  analysis  [32.33]  is  beyond 
the  scope  of  this  report;  we  would  like  to  show  initially  that  with  a  computer  simulation  of  RF- 
illuminated  phantom  muscle  tissue  [34j  adaptive  nulling  can  improve  the  thermal  dose  distribution. 
The  thermal  analysis  used  in  this  report  is  similar  to  the  method  presented  in  Zhang,  et  al.  [35], 
where  blood  flow  effects  are  ignored.  In  practice,  undesired  high-temperature  regions  away  from 
the  focus  can  occur  inside  the  target  volume.  For  example,  scar  tissue,  which  has  a  decreased 
blood  flow  rate,  will  tend  to  heat  up  more  rapidly  than  normal  tissue.  In  the  proposed  adaptive 
hyperthermia  concept,  electric-field  nulls  are  used  to  reduce  the  power  delivered  to  potential  hot 
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spots.  Computer  simulations  show  that  noninvasive  field  probes  can  be  used  to  eliminate  hot  spots 
interior  to  the  phantom  target  tissue. 
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Figure  3.  Two  applications  of  focused  near-field  adaptive  nulling:  (a)  Radar  system 
testing  and  (b)  Hyperthermia  treatment  of  cancer. 


An  antenna  analysis  code  (called  WIRES)  capable  of  analyzing  a  wude  variety  of  antenna 
or  radar  cross  section  problems  has  been  developed  by  J.H.  Richmond  [36-38].  The  software  can 
analyze  complex  geometries  involving  multiple-connected  and/or  isolated  wires  and  is  well-suited  to 
analyzing  low-frequency  phased  arrays  radiating  into  an  infinite  homogeneous  conducting  medium. 
WIRES  is  a  moment-method  code  that  uses  the  electric  field  integral  equation  (EFIE)  to  enforce 
the  boundary  condition  of  the  tangential  electric  field  being  zero  at  the  surface  of  the  antenna  of 
interest.  The  moment-method  basis  and  testing  functions  used  in  this  code  are  piecewise  sinusoidal. 
For  a  complete  description  of  the  theory  and  capabilities  of  the  WIRES  code,  the  reader  is  referred 
to  the  user’s  manual  [37],  The  WIRES  software  has  been  modified  by  the  author  to  include  near- 
field  focusing  and  adaptive  nulling  capabilities.  This  modified  code  has  been  used  to  analyze  the 
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near-field  and  far-field  adaptive  nulling  performance  of  thin- wire  phased  arrays  in  free  space  [27,29). 
A  new  version  of  the  thin-wire  code  that  can  analyze  adaptive  hyperthermia  arrays  in  an  infinite 
homogeneous  conducting  medium  has  recently  been  written  and  is  documented  in  this  report. 

This  report  is  organized  in  the  following  manner.  Section  2  discusses  the  concept  for  a 
noninvasive  adaptive  hyperthermia  system  and  describes  in  detail  the  theory  used  in  performing 
the  focused  near-field  adaptive  nulling  simulation.  Pertinent  theory  for  wave  propagation  in  a 
conducting  medium  is  reviewed.  A  brief  description  of  the  transient  thermal  analysis  approach  used 
to  compute  the  temperature  distribution  in  a  target  is  given.  Section  3  presents  simulation  results 
for  a  fully  adaptive  eight-element  hyperthermia  ring  array  operating  at  120  MHz.  The  calculated 
received  RF-power  distributions,  before  and  after  adaptive  nulling,  at  a  short-dipole  field  probe 
with  both  the  field  probe  and  the  hyperthermia  array  embedded  in  homogeneous  muscle  tissue, 
are  presented.  The  RF  power  distribution  is  then  fed  into  transient  thermal  analysis  software 
which  computes  the  two-dimensional  temperature  distribution  in  an  elliptical  muscle-tissue  target 
surrounded  with  a  constant-temperature  water  bolus.  The  temperature  distribution  obtained  by 
using  multiple  independent  adaptive  nulls  clearly  indicates  the  elimination  of  undesired  hot  spots. 
Conclusions  are  made  in  Section  4.  Performance  degradation  due  to  insufficient  number  of  auxiliary 
null  positions  is  described  in  Appendix  A.  Appendix  B  contains  data  files  and  source  code  for  the 
focused  near-field  adaptive  nulling  computer  simulations. 
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2.  THEORY 


2.1  NONIN VASIVE  ADAPTIVE  HYPERTHERMIA  SYSTEM  CONCEPT 

The  concept  of  a  noninvasive  adaptive- nulling  hyperthermia  system  [39,40]  is  shown  in  Fig¬ 
ure  4.  Theoretically,  to  generate  the  desired  field  distribution  in  a  clinical  adaptive  hyperthermia 
system,  receiving  sensors  are  positioned  as  close  as  possible  to  the  focus  (tumor  site)  and  to  where 
high  temperatures  arc  to  be  avoided  (such  as  near  the  spinal  cord  and  scar  tissue).  As  simulations 
in  this  report  will  show,  for  an  annular  array  configuration  the  receiving  sensors  can  be  located 
noninvasively  on  the  surface  (skin)  of  the  target.  Initially,  the  hyperthermia  array  is  focused  to 
produce  the  required  field  intensity  at  the  tumor.  An  invasive  probe  is  required  to  achieve  the 
optimum  focus  at  depth.  To  avoid  undesired  hot  spots,  it  is  necessary  to  minimize  the  power  re¬ 
ceived  at  the  desired  null  positions  and  to  constrain  the  array  weights  to  deliver  a  required  amount 
of  transmitted  or  focal  region  power.  The  adaptive  array  weights  (with  gain  g  and  phase  <j>)  are 
controlled  by  either  the  SMI  algorithm  or  a  gradient  search  algorithm  to  rapidly  (within  seconds) 
form  the  nulls  before  a  significant  amount  of  target  heating  takes  place.  With  this  adaptive  tech¬ 
nique,  it  should  be  possible  to  avoid  unintentional  hot  spots  and  maintain  a  therapeutic  thermal 
dose  distribution  at  the  tumor. 

Figure  5(a)  shows  a  specific  example  of  an  eight-element  hyperthermia  ring  phased  array  with 
a  target  cross  section  at  the  prostate  level  [9].  Figure  5(b)  shows  an  elliptical  phantom  target  which 
is  used  to  model  the  prostate-level  cross  section.  A  noninvasive  adaptive  nulling  system  is  achieved 
by  placing  auxiliary  sensors  1.2.---,  Naux  on  the  target  skin  as  shown.  The  null  zones  centered  at 
each  auxiliary  probe  naturally  extend  into  the  elliptical  target  region  to  eliminate  undesired  hot 
spots.  (An  example  of  a  two-dimensional  near-field  null  zone  is  shown  in  Figure  10(b)  of  Fenn  [28]). 
The  width  of  each  null  zone  is  directly  related  to  the  depth  of  each  null.  The  depth  of  each  null 
(sometimes  referred  to  as  the  amount  of  cancellation)  is  directly  related  to  the  SNR  at  the  sensor 
position.  A  low  SNR  produces  a  small  amount  of  nulling,  and  a  high  SNR  produces  a  large  amount 
of  nulling.  The  resolution  or  minimum  spacing  between  the  focus  and  null  position  is  normally  equal 
to  the  half-power  beamwidth  of  the  antenna.  The  resolution  is  enhanced  somewhat  by  using  weak 
nulls  whenever  the  separation  between  the  null  and  focus  is  closer  than  the  half-power  beamwidth. 

Figure  6  summarizes  the  thermal  simulation  for  a  hyperthermia  ring  array  transmitting  into 
a  homogeneous  elliptical  target  region.  The  tumor  site  is  assumed  to  be  at  the  center  of  the  ellipse. 
The  thermal  distribution  in  Figure  6(a)  contains  two  undesired  hot  spots  to  the  left  and  right  of 
the  focus.  In  Figure  6(b).  adaptive  nulling  is  applied  and  the  hot  spots  are  eliminated.  The  details 
of  the  computer  simulations  are  contained  in  Section  3. 

In  the  near-field  nulling  technique  described  here,  it  is  assumed  that  the  hyperthermia  phased- 
array  antenna  is  focused  (as  it  normally  is)  in  the  near  field  and  that  a  main  beam  and  possibly 
sidelobes  are  formed  in  the  target.  In  this  report,  it  is  assumed  that  phase  focusing  is  used  to 
produce  the  desired  quiescent  main  beam.  Figure  7  shows  a  primary  (calibration)  probe  antenna 
located  at  a  desired  focal  point  of  an  array.  The  array  can  maximize  the  signal  received  by  the 
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Figxm  5.  (a)  Transverse  cross  section  at  the  prostate  level.  An  eight- element  hyperther¬ 

mia  rtnq  array  of  dipole  elements  is  used  to  irradiate  a  target  body  with  sufficient  power 
to  elevate  the  temperature  of  the  target  body  to  a  therapeutic  level,  (b)  Elliptical  phantom 
target  used  to  model  the  prostate  (target)-level  cross  section.  Shown  are  four  noninvasive 
auxiliary  E-field  probes  that  are  used  in  forming  null  zones  with  widths  controlled  by  the 
individual  null  depths.  The  null  zones  extend  into  the  body  to  reduce  the  internal  electric 
field  at  specified  locations. 
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Figure  fi.  Simulated  two-dimensional  thermal  pattern  at  time  t  =  20  minutes  in  elliptical 
phantom  muscle-tissue  target  surrounded  with  10°  C  constant-temperature  water  bolus. 
The  incident  RF  power  distribution  is  at  120  MHz  and  the  initial  temperature  of  the 
phantom  is  25°  C  Temperature  contour  levels  are  given  in  2°  C  steps,  (a)  Before  nulling: 
hot  spots  on  the  left  and  right  sides  of  the  target  are  present,  (b)  After  nulling:  hot  spots 
art  eliminated 
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calibration  probe  by  adjusting  its  phase  shifters  so  that  the  observed  element-to-element  phase 
variation  is  removed.  One  way  to  do  this  in  a  numerical  simulation  is  to  choose  a  reference  path 
length  as  the  distance  from  the  focal  point  to  the  phase  center  of  the  array.  This  distance  is 
denoted  rp  and  the  distance  from  the  focal  point  to  the  nth  array  element  is  denoted  .  The 
voltage  received  at  a  probe  (located  at  the  focal  point)  due  to  the  nth  array  element  is  computed 
here  using  the  method  of  moments  as  described  in  Section  2.3.  To  maximize  the  received  voltage  at 
the  focal-point  probe  output,  it  is  necessary  to  apply  the  phase  conjugate  of  the  observed  signals, 
due  to  the  array  elements,  at  the  transmit  array.  The  resisting  near-field  radiation  pattern  will 
have  a  main  beam  and  sidelobes.  The  main  beam  will  be  pointed  at  the  array  focal  point,  and 
sidelobes  will  exist  at  angles  away  from  the  main  beam.  Auxiliary  probe  antennas  can  then  be 
placed  at  the  desired  null  positions  in  the  quiescent  sidelobe  region.  These  sidelobes  are  where 
tissue  hot  spots  are  likely  to  occur:  they  are  nulled  by  the  adaptive  nulling  algorithms  described 
below. 
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Figurt  7.  Adaptive  transmit  phased-array  antenna  near-field  focusing  concept. 
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2.2  ADAPTIVE  TRANSMIT  ARRAY  FORMULATION 


Consider  the  hyperthermia  array  and  probe  geometry  shown  in  Figure  7.  Typically,  the 
hyperthermia  array  contains  A’  identical  antenna  elements.  The  input  signal  to  each  of  the  N 
array  elements  is  obtained  from  me  weighted  signal  distributed  by  the  power  divider.  The  number 
of  adaptive  channels  is  denoted  M ,  and  for  a  fully  adaptive  array  M  —  Ar.  Note:  the  report 
will  concentrate  on  fully  adaptive  arrays  (arrays  in  which  all  transmit  array  elements  have  adaptive 
weight  control):  however,  the  computer  simulation  program  is  developed  for  both  fully  adaptive  and 
sidelobe  canceller  arrays.  In  this  report,  ideal  transmit  weights  (a  complex  voltage  gain  vector)  are 
assumed  in  the  compute:  simulation,  with  w  =  (w i.ti’2.  ■  •  •  ,w\)T  denoting  the  adaptive  channel 
weight  vector  as  shown  in  Figure  4.1  (Superscript  T  means  transpose).  To  generate  adaptive  nulls, 
the  transmit  weights  (phase  and  gain)  are  controlled  by  either  the  SMI  algorithm  or  a  gradient 
search  algorithm  The  SMI  algorithm  has  the  flexibility  to  operate  in  either  open-  or  closed-loop 
feedback  modes  41  :  the  gradient  search  algorithm  operates  only  in  a  feedback  mode. 

2.2.1  Sample  Matrix  Inversion  Algorithm 

With  the  SMI  algorithm  23  .  the  fundamental  quantities  required  to  fully  characterize  the 
incident  field  for  adaptive  nulling  purposes  are  the  adaptive  channel  cross  correlations.  To  imple¬ 
ment  this  algorithm  it  is  necessary  to  know  the  complex  received  voltage  at  the  auxiliary  probes. 
For  example,  the  moment-method  formulation  (described  in  Section  2.3)  allows  computation  of  the 
complex-received  voltage  at  the  auxiliary  probes. 

The  block  diagram  in  Figure  8  shows  how  the  SMI  algorithm  is  used  in  the  hyperthermia 
analysis.  The  diagram  indicates  four  performance  measures  used  to  quantify  the  computer  sim¬ 
ulations:  electric-field  distribution,  covariance  matrix  eigenvalues,  adaptive  transmit  weights,  and 
interference  cancellation.  The  calculation  of  these  performance  measures  is  described  in  detail 
below. 

Referencing  Figure  7.  let  a  spheric.  '  wavefront  be  incident  at  the  2th  probe  antenna,  due  to 
each  array  element  (radiating  one  at  a  time  with  a  unity-amplitude  reference  signal),  which  results 
in  a  set  of  probe- received  complex  voltages  denoted  ?  !■(>.  •  •  • .  r\.  The  cross  correlation  R’mn  of 

the  received  voltages  due  to  the  rnth  and  nth  adaptive  transmit  channels  at  the  ft h  probe  is  given 
b\ 


R  mr ,  ~  t  rn  t  n  1 


(1) 


’The  effects  of  transmit-weight  quantization  and  transmit-weight  random  errors  are  addressed  in 
the  sof*  "are 
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Figure  8.  Implementation  of  the  SMI  algorithm  for  an  adaptive  hyperthermia  system. 
Four  performance  measures  are  used  in  evaluating  the  adaptive  response  of  the  system: 
electric -field  distribution,  covariance  matrix  eigenvalues,  adaptive  transmit  weights,  and 
interference  cancellation. 
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where  ’  means  complex  conjugate  and  E{-)  means  mathematical  expectation.  Because  vm  and 
vn  represent  voltages  of  the  same  waveform  but  at  different  times,  R'mn  is  also  referred  to  as  an 
autocorrelation  function.  Note:  for  convenience,  in  Equation  (1)  the  superscript  i  in  vm  and  in  t;n 
has  been  omitted. 

In  the  frequency  domain,  assuming  the  transmit  waveform  has  a  band-limited  white  noise 
power  spectral  density  (as  commonly  assumed  in  jamming  of  a  radar  system),  Equation  (1)  can  be 
expressed  as  the  frequency  average 

*™^j'2Vm(fK(f)df  ,  (2) 

where  B  =  fi  -  f\  is  the  nulling  bandwidth  and  /  is  the  frequency.  It  should  be  noted  that 
v m(f  )  takes  into  account  the  transmit  wavefront  shape.  wThich  is  spherical  for  the  hyperthermia 
application.  For  the  special  case  of  a  continuous  wave  (CW)  transmit  waveform,  as  normally  used 
in  hyperthermia,  the  cross  correlation  reduces  to 

Rmn  =  vrn(fo)Vn(fo)  i  (3) 

where  fQ  is  the  transmit  frequency  of  the  hyperthermia  array. 

Let  the  channel  or  interference  covariance  matrix  be  denoted  R.  (Note:  in  hyperthermia, 
interference  is  used  to  refer  to  the  signals  received  at  the  auxiliary  probes.  The  undesired  hot 
spots  can  be  thought  of  as  interfering  with  the  therapy.)  If  there  are  JauT  independent  desired 
null  positions  or  auxiliary  probes,  the  JauI-probe  covariance  matrix  is  the  sum  of  the  covariance 
matrices  observed  at  the  individual  probes;  that  is. 

Jan; 

R=Y.K-I  (4) 

2=] 

where  R  is  the  sample  covariance  matrix  observed  at  the  zth  probe  and  I  is  the  identity  matrix 
used  to  represent  the  thermal  noise  level  of  the  receiver.  (It  should  be  noted  that  both  Nauz  and 
Jaux  can  be  used  interchangeably  to  describe  the  number  cf  auxiliary  nulling  sensors  or  received- 
interference  signal  positions.  The  symbol  Jaux  is  perhaps  more  descriptive  and  can  equivalently  be 
thought  of  as  the  number  of  positions  where  jamming  signals  are  nulled). 

Prior  to  generating  an  adaptive  null,  the  adaptive  channel  weight  vector,  w,  is  chosen  to 
synthesize  a  uesired  quiescent  radiation  pattern.  vVhen  nulling  is  desired,  the  optimum  set  of 
transmit  weights  to  form  an  adaptive  null  (or  nulls),  denoted  wa ,  is  computed  [23]  by 

wa  =  R  ~  ^  xv  q  .  (5) 
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where  _1  means  inverse  and  wq  is  the  quiescent  weight  vector.  During  array  calibration,  the 
normalized  quiescent  transmit  weight  vector,  with  element  1  radiating,  is  chosen  to  be  wq  = 
(1, 0. 0.  ■  •  • ,  0)T;  that  is,  the  transmit  channel  weight  of  element  1  is  unity  and  the  remaining  transmit 
channel  weights  are  zero.  Similar  weight  settings  are  used  to  calibrate  the  remaining  transmit 
elements.  For  a  fully  adaptive  annular  array  focused  at  the  origin  in  homogeneous  tissue,  the 
normalized  quiescent  weight  vector  is  simply  wq  —  (1, 1, 1, • • • , 1)T.  Commonly,  the  weight  vector 
is  constrained  to  deliver  a  required  amount  of  power  to  the  hyperthermia  array  or  to  the  tumor. 
For  simplic:ty  in  thj  software  used  to  analyze  the  hyperthermia  array,  the  weights  are  constrained 
such  that 


M 


Y. 

m  =  l 


j 


2 


=  1 


(6) 


where  uy,  is  the  transmit  weight  for  the  nth  element.  Note:  in  the  computer  simulations,  the 
electric  field  due  to  the  normalized  weight  vector  is  scaled  appropriately  to  deliver  the  required 
amount  of  power  to  the  tissue  so  that  a  desired  focal-region  temperature  level  is  achieved  after  t 
minutes. 

The  summation  of  power  received  at  the  probes  is  given  by 


p  =  wf  Rw  .  (7) 

where  f  means  complex  conjugate  transpose.  The  interference-plus-noise-to-noise  ratio,  denoted 
INR.  is  computed  as  the  ratio  of  the  auxiliary  probe  array  output  power  [defined  in  Equation  (7)] 
with  the  transmit  signal  present  to  the  probe  array  output  power  with  only  receiver  noise  present; 
that  is. 


/.Vi? 


wj  Rw 
w^w 


(8) 


The  adaptive  array  cancellation  ratio,  denoted  C,  is  defined  here  as  the  ratio  of  the  summation  of 
probe-received  power  after  adaptation  to  the  summation  of  probe-received  power  before  adaptation; 
that  is. 


C  =  —  (9) 

Substituting  Equation  (7)  in  Equation  (9)  yields 


15 


Wa^Rwa 


(10) 


WqRwq 

Next,  the  covariance  matrix  defined  by  the  elements  in  Equations  (2)  or  (3)  is  Hermitian 
(that  is,  R=  R*).  which,  by  the  spectral  theorem,  can  be  decomposed  in  eigenspace  [42]  as 

M 

’  (U) 

Jt=  l 


where  X^.k  —  1.2.  •■•.A/  are  the  eigenvalues  of  R,  and  e*.,/c  =  1,2  are  the  associated 

eigenvectors  of  R.  The  interference  covariance  matrix  eigenvalues  (Ai,  A2,  •  •  • ,  Am)  are  a  convenient 
quantitative  measure  of  the  use  of  the  adaptive  array  degrees  of  freedom  [43].  The  amplitude  spread 
between  the  largest  and  smallest  eigenvalues  is  a  quantitative  measure  of  the  dynamic  range  of  the 
interference  (hot  spot)  signals. 

2.2.2  Gradient  Search  Algorithm 

Under  conditions  where  only  the  probe- received  voltage  amplitude  is  measured  (as  in  the 
BSD-2000  system  previously  mentioned),  it  is  appropriate  to  consider  a  gradient  search  algorithm 
to  minimize  the  interference  pow’er  at  selected  positions.  The  gradient  search  is  used  to  iteratively 
control  the  transmit  weights  so  that  the  RF  signal  received  by  the  probe  array  is  minimized.  The 
transmit  array  weights  (gain  and  phase)  are  adaptively  changed  in  small  increments  and  the  probe 
array  output  power  is  monitored  to  determine  weight  settings  that  reduce  the  output  power  most 
rapidly  to  a  null.  The  mathematical  formulation  for  the  gradient  search  is  developed  in  a  straight¬ 
forward  manner  [44]  and  will  not  be  considered  further  here.  Only  the  sample  covariance  matrix 
inversion  algorithm  has  been  implemented  in  the  focused  near-field  nulling  simulation  presented  in 
this  report. 

2.3  MOMENT-METHOD  FORMULATION 

This  section  describes  a  method  of  moments  formulation  to  compute  the  probe-received  volt¬ 
ages  in  Equation  (2)  due  to  the  transmitting  hyperthermia  phased-arrav  antenna  in  am  infinite 
homogeneous  conducting  medium.  The  medium  is  described  by  the  three  parameters  n,  e,  and  a, 
which  are  discussed  in  the  next  section.  The  formulation  given  here  is  analogous  to  that  developed 
under  array-receiving  conditions  for  an  adaptive  radar  [29].  The  software  used  to  analyze  a  hyper¬ 
thermia  array  is  based  on  the  receive-array  analogy  but  the  theory  presented  below  is  given  in  the 
context  of  a  transmit  an  ay.  Referring  to  Figure  9,  assume  that  each  transmit  element  is  fed  with 
a  generator  having  a  known  impedance  Zi .  Let  v°  c3  represent  the  open-circuit  voltage  at  the  jth 
probe  due  to  the  nth  transmit-array  element.  Here,  the  jth  probe  can  denote  either  the  focal  point 
calibration  probe  or  one  of  the  auxiliary  probes  used  to  null  a  sidelobe.  The  number  of  auxiliary 
probes  is  denoted  by  Jaux .  The  receive  probe  is  assumed  to  be  terminated  in  an  impedance  ZT. 
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Figure  9.  Transmit  array  and  receive  antenna  probe.  Open-circuit  mutual  impedance 
between  array  elements  is  denoted  Zmn.  The  open-circuit  voltage  at  the  probe  is  computed 
from  the  array  terminal  currents  and  from  Z]n.  the  open-circuit  mutual  impedance  between 
the  nth  array  element  and  the  jth  probe  antenna. 


172668-10 


Figure  10.  Equivalent  circuit  model  for  receive  antenna  probe. 
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Next,  let  Z  denote  the  open-circuit  mutual  impedance  matrix  (with  dimensions  N  x  N  for  the 
jV-element  array).  The  open-circuit  mutual  impedance  between  array  elements  m  and  n  is  denoted 
Zm, n-  It  is  assumed  that  multiple  interaction  between  the  hyperthermia  array  and  the  probe  sensor 
can  be  neglected.  Thus,  the  hyperthermia  array  terminal  current  vector  i  can  be  computed  in 
terms  of  the  transmit  weights  w  as 

i  =  [Z  +  Zi/]-1m  (12) 

Next,  let  Zl  be  the  open-circuit  mutual  impedance  between  the  jth  probe  and  the  nth  array 
element.  The  induced  open-circuit  voltage  at  the  jth  receive  probe,  due  to  the  nth  array 
element  transmit  current  in,  can  then  be  expressed  as 


ro  c  -  7]  ■  i 

Ln,j  ~  ^ n  ln 


In  matrix  form,  the  induced  open-circuit  probe-voltage  matrix  is 


V  probe  ~~  %  probe, array  l 


(13) 


(14) 


or 


Vprob(  ~  Z probe.array  +  Zl& ]  tV  ,  (15) 

where  Zprobe.array  is  a  rectangular  matrix  of  order  JauT  x  N  for  the  open-circuit  mutual  impedance 
between  the  probe  array  and  the  hyperthermia  array.  Note:  the  jth  row  of  the  matrix  Zprobe.array  is 
written  as  (Zj ,  Z32.  •  ■  ■ ,  Z\- ) .  where  j  =  1,  2.  •  •  • ,  Jaux.  The  receive  voltage  matrix  is  then  computed 
by  the  receiving  circuit  equivalence  theorem  for  an  antenna  [50].  The  receive-antenna  equivalent 
circuit  is  depicted  in  Figure  10,  where  it  is  readily  determined  that 


' probe 


=  V 


Zr 


probe 


Zr 


(16) 


where  Ztn  is  the  input  impedance  of  the  probe.  It  should  be  noted  that  the  matrix  is  a 

column  vector  of  length  Jaux  and  ujec  is  the  jth  element  of  the  matrix.  The  probe-receive  current 
matrix  is  given  by 


*pro6€ 


=  V 


probe  7 


1 

7Z 


(17) 
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The  jth  element  of  the  column  vector  is  denoted  ijec,j  =  1, 2,  •  •  ■ ,  «7aui.  Finally,  the  power 
received  by  the  jth  probe  is 

Pyc=l-Re(vYc-iTc’)  ,  (18) 

where  Re  means  real  part.  Substituting  Equations  (16)  and  (17)  into  Equation  (18)  yields 


Pj 


o.c.|2  RejZr) 

2  J  \Zin  +  Zr\2 


(19) 


The  total  interference  power  received  by  the  auxiliary  probe  array  is  given  by 


prec  =  £  prr  ■  (20) 

>=i 

The  incident  electric  field  E  is  related  to  the  open-circuit  voltage  voc'  by  the  effective  height 
h  of  the  probe  antenna  [48]  as 


voc  =  hE 


(21) 


If  the  length  Lp  of  the  probe  antenna  is  approximately  0.1A  or  less,  the  current  distribution  is 
triangular  and  the  effective  height  is  h  =  0.5 Lp  [48].  Thus,  for  a  short-dipole  probe  the  open-circuit 
voltage  can  be  expressed  as 

voc  =  yE  (22) 

It  then  follows  from  Equation  (22 j  that  the  E  field  for  a  short-dipole  probe  at  position  ( x,y,z )  is 
given  by 


E(x.y.z) 


2voc(x,y.z) 

Lp 


(23) 


Finally,  the  quiescent  and  adapted  E-field  radiation  patterns  are  computed  using  the  quiescent  and 
adapted  weight  vectors  wq  E.nd  wa,  respectively,  in  Equations  (15)  and  (23). 

The  moment-method  expansion  and  testing  functions  are  assumed  to  be  sinusoidal.  The  open- 
circuit  mutual  impedances  in  Equation  (15)  between  thin- wire  dipoles  in  a  homogeneous  conducting 
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medium  are  computed  based  on  subroutines  from  a  well-known  moment-method  computer  code  [38]. 
In  evaluating  ZJn  for  the  jth  auxiliary  probe,  double  precision  computations  are  used. 

As  mentioned  previously,  the  array  is  calibrated  (phased  focused)  initially  using  a  short  dipole 
at  the  focal  point.  To  accomplish  this  numerically,  having  computed  vTy£cut,  the  receive  array,  weight 
vector  w  will  have  its  phase  commands  set  equal  to  the  conjugate  of  the  corresponding  phases  in 
vfocua-  Transmit  antenna  radiation  patterns  are  obtained  by  scanning  (moving)  a  dipole  probe 
with  half-length  l  in  the  near-field  and  computing  the  receive  probe-voltage  response. 

The  received  voltage  matrix  for  the  jth  probe  (denoted  rjec)  is  computed  at  K  frequencies 
across  the  nulling  bandwidth.  Thus,  Uj£c(/i), fjec(/2),  •  •  • , Wjec(//c)  are  needed.  In  this  report,  the 
impedance  matrix  is  computed  at  K  frequencies  and  is  inverted  K  times.  The  probe  covariance  ma¬ 
trix  elements  are  computed  by  evaluating  Equation  (2)  numerically,  using  Simpson’s  rule  numerical 
integration.  For  multiple  auxiliary  probes,  the  covariance  matrix  is  evaluated  using  Equation  (3). 
Adaptive  array  radiation  patterns  are  computed  by  superimposing  the  quiescent  radiation  pattern 
with  the  weighted  sum  of  auxiliary-channel-received  voltages. 

2.3.1  WAVE  PROPAGATION  IN  CONDUCTING  MEDIA 

To  gain  insight  into  the  effect  of  a  lossy  medium  on  the  propagation  of  an  electromagnetic 
wave,  it  is  useful  to  review  certain  fundamental  equations  which  govern  the  field  characteristics.  In 
a  conducting  medium,  Maxwell’s  curl  equations  in  time-harmonic  form  are 


V  x  H  =  J  +  jutE 

(24) 

Vx£=  -j^nH  , 

(25) 

where  E  and  H  are  the  electric  and  magnetic  fields,  respectively.  J  is  the  conduction  current 
density,  u;  =  27 -f  is  the  radian  frequency,  e  is  the  permittivity  of  the  medium,  and  n  is  the 
permeability  of  the  medium.  The  permittivity  is  expressed  as  e  =  ere0,  where  er  is  the  dielectric 
constant  (relative  permittivity)  and  eQ  is  the  permittivity  of  free  space.  Similarly,  /j  =  HrHo,  where 
Hr  is  the  relative  permeability  and  /i0  is  the  permeability  of  free  space.  For  a  medium  with  electrical 
conductivity  cr,  J  and  E  are  related  as 

J  —  aE  (26) 

Substituting  Equation  (26)  into  Equation  (24)  yields 

VxH  =  {a  +  j^t)E  (27) 

From  Equations  (24)  and  (25).  the  vector  wave  equation  in  terms  of  E  is  derived  as 
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V2E-~i2E  =  0 


(28) 


It  is  readily  shown  [49]  that 


1  = 


(29) 


The  quantity  a/ioe  is  referred  to  as  the  loss  tangent.  It  is  common  to  express  the  complex  propa¬ 
gation  constant  as 


-)  =  a  +  j3 


(30) 


where  a  is  the  attenuation  constant  and  0  is  the  phase  constant.  The  constants  a  and  0  are  found 
by  setting  Equation  (29)  equal  to  Equation  (30)  and  then  squaring  both  sides,  equating  the  real 
and  imaginary  parts,  and  solving  the  pair  of  simultaneous  equations,  with  the  result 


-  r  /,  - 

v/2  V 


a  =  ^-{J  l  +  (_)2_i} 


1/2 


3  = 


TT‘V1  +  (^  +  1) 


1/2 


(31) 

(32) 


The  wavelength  A  in  the  lossy  dielectric  is  then  computed  from 


The  intrinsic  wave  impedance  r?  is  given  by 


(33) 


77  = 


\ 


a  -r 


/t  _ 

\  1 


(34) 


The  instantaneous  powder  density  of  the  electromagnetic  field  is  given  by  Poynting's  vector, 
denoted  P, 


P  =  -ExH‘ 

2 


(35) 
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which  has  units  of  (W/m2).  The  time-average  power  flow  density  is  equal  to  the  real  part  of  the 
complex  Poynting's  vector.  The  time-average  power  dissipation  per  unit  volume  Pd  (W/m3)  is 
derived  from  Maxwell’s  equations,  with  the  result 

Pd=l-E-  r  =  \a\E\2  (36) 

The  specific  absorption  rate  (SAR)  [2]  is  the  power  dissipated  or  absorbed  per  unit  mass  (W/kg) 
of  the  medium  (tissue),  or 


=  ,  (37) 

P  2  p 

where  p  is  the  density  of  the  medium  in  kg/m3. 

It  is  convenient  to  have  a  simple  equation  for  computing  the  propagation  loss  between  any 
two  points  in  the  near  field  of  an  isolated  transmitting  antenna.  Thus,  mutual  coupling  effects  are 
ignored  for  the  time  being.  Consider  a  time-harmonic  source  radiating  a  spherical  wave  into  an 
infinite  homogeneous  conducting  medium.  For  an  isotropic  radiator,  and  suppressing  the  eJU,(  time 
dependence,  the  electric  field  as  a  function  of  range  r  can  be  expressed  as 

E(r)  =  E0—  ,  (38) 

T 

where  Ea  is  a  constant. 

For  a  source  at  the  origin,  the  amplitude  of  the  electric  field  at  range  r i  is  given  by 

p-OTl 

I  £(r0!  =  E0 -  ,  (39) 

n 


and  at  range  r2  by 


\E(r2)\  =  E0 


e 


ar  2 


r2 


(40) 


The  total  propagation  loss  between  ranges  ri  and  r2  is  found  by  taking  the  ratio  of  Equations  (40) 
and  (39).  or 


\E{t2)  _  ^e-o(r2-ri ) 

|£(n)i  r2 
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The  field  attenuation  Aa  in  dB  from  range  n  to  range  r2  due  to  the  lossy  dielectric  is  simply 


AQ  =  201og10(e-Q^-r>>)  . 

Similarly,  the  1/r  attenuation  loss  Ar  in  dB  is 


(42) 


At  =  20  log10  — 
r2 


(43) 


2.3.2  THERMAL  MODELING  OF  AN  INHOMOGENEOUS  TARGET 

A  thermal  analysis  computer  program  called  the  transient  thermal  analyzer  (TTA),  developed 
by  Arthur  D.  Little.  Inc.,  has  been  used  at  Lincoln  Laboratory  for  space  payload  analysis  since  1969. 
Recent  documentation  for  TTA  exists  [51,52].  The  important  details  of  the  analysis  are  described 
in  this  section.  The  TTA  software  is  used  in  the  next  section  to  accomplish  the  thermal  modeling  of 
homogeneous  muscle  tissue  surrounded  by  a  constant-temperature  water  bolus.  Thermal  modeling 
techniques  have  been  described  in  detail  [53]. 

The  TTA  program  uses  the  finite-difference  technique  to  solve  a  set  of  nonlinear  energy  balance 
equations.  Consider  a  system  of  interconnected  nodes  that  model  an  inhomogeneous  volume  for 
which  the  temperature  T,  of  the  zth  node  is  to  be  determined.  The  heat-balance  equation,  which 
is  solved  by  TTA.  is  expressed  as 

-AW +  A/^=  0  ,  (44) 


where  is  the  net  outward  heat  flow  from  node  i  in  the  direction  of  node  j,  Pi(t )  is  the  power 
into  node  i  at  time  t,  and  M ,  is  the  thermal  mass  (mass  times  specific  heat)  of  node  i.  Figure  11 
shows  an  electric  circuit  analog  which  is  used  to  model  the  two-dimensional  thermal  characteristics 
of  the  material  volume.  Power  P,  in  watts  is  delivered  to  the  z'th  node.  Capacitor  Cj  (with  units 
Joules /0C)  is  used  to  model  the  thermal  capacitance  at  the  zth  node.  Resistor  (with  units 
°C/\Y)  is  used  to  model  the  heat  resistance  between  nodes  z  and  j. 

With  a  spacing  of  A /  between  nodes  (assuming  cubic  cells),  the  values  of  Rij,  Ci  and  P*  are 
computed  as 


Ri,j  = 


1 

k~\l 


(45) 


where  fctJ  is  the  thermal  conductivity  (with  units  W/m°C)  between  nodes  z  and  j\ 


23 


1 72008-1 1 


R  =  RESISTOR 
C  =  CAPACITOR 
P  =  POWER 


R  =  —  [;C  /  Wl.  k  =  THERMAL  CONDUCTIVITY 
kAi  1  1 

C  =  pCP  (S(f  [J/SC],  p  =  DENSITY, CP  =  SPECIFIC  HEAT 
P  =  SAR  p  |^/)3  IWl.SAR  =  —  |E|2  (Specific  Absorption  Rate) 

'  2  p 

a  =  ELECTRICAL  CONDUCTIVITY 
)E;  =  MAGNITUDE  OF  ELECTRIC  FIELD 


Figure  11.  Equivalent  circuit  model  used  in  thermal  analysis.  The  inhomogeneous  body 
of  interest  is  modeled  with  a  grid  of  nodes  interconnected  by  resistors.  Each  node  has  an 
associated  power  source  and  capacitor. 
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(46) 


C,  =  piCptiAl)3 

where  Cp,  is  the  specific  heat  at  the  2th  node  and  p,  is  the  density  (kg/m3)  at  the  2th  node;  and 
P,  =  (SAR)lPl(Alf  ,  (47) 

where  (5.47?),  is  the  SAR  for  the  2th  node,  which  is  given  by 

(SAR),  =  ^\Et\2  .  (48) 

lpt 

where  o,  is  the  electrical  conductivity  of  the  2th  node  and  \E,[  is  the  magnitude  of  the  electric 
field  delivered  by  the  hyperthermia  array  to  the  it h  node.  Note;  in  substituting  Equation  (48) 
into  Equation  (47).  the  density  p,  cancels.  Thus,  an  equivalent  approach  to  computing  the  power 
delivered  to  the  ith  node  is  written  in  terms  of  the  time-average  power  dissipated  per  unit  volume 
of  the  2th  node  (denoted  Pdl)  as 

P  =  Pdt(  A/)3  (49) 

A  block  diagram  showing  how  TTA  is  used  in  the  hyperthermia  analysis  is  given  in  Figure  12. 
First,  the  method  of  moments,  controlled  by  the  SMI  nulling  algorithm,  is  used  to  compute  the 
electric  field  radiation  pattern  throughout  a  homogeneous  region  (muscle  tissue).  This  power  dis¬ 
tribution  is  then  read  into  the  TTA  program,  which  computes  the  temperature  distribution  inside 
an  elliptical  muscle-tissue  target  surrounded  with  a  constant-temperature  water  bolus. 
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Figure  12.  The  moment-method  simulation's  use  with  the  transient  thermal  analysis.  RF 
power  distribution  over  the  target  region  is  converted  into  a  time -dependent  temperature 
distribution. 
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3.  SIMULATION  RESULTS 


To  demonstrate  the  effectiveness  of  focused  near-field  adaptive  nulling  in  reducing  undesired 
hot  spots,  several  computer  simulations  are  presented.  The  E-field  simulations  are  for  the  signal 
received  by  a  short-dipole  probe  due  to  a  transmitting  phased  array  embedded  in  an  infinite  homo¬ 
geneous  lossy  dielectric  (muscle  tissue).  However,  the  thermal  simulations  are  for  an  elliptical  target 
(muscle  tissue)  surrounded  by  a  constant-temperature  water  bolus.  Because  the  FF  wavelengths 
in  the  target  and  water  bolus  are  similar,  the  El-field  simulations  are  believed  to  give  a  reasonable 
approximation  to  the  field  distribution  inside  the  elliptical  target. 

The  E-field  calculation  in  the  assumed  infinite  homogeneous  medium  introduces  additional 
field  attenuation  not  present  in  a  clinical  hyperthermia  system  with  an  annular  array  transmitting 
through  a  water  bolus  into  a  patient.  (As  mentioned  earlier,  the  water  bolus  has  very  little  RF 
propagation  loss.)  In  addition,  the  transmit  array  weights  are  normalized  according  to  Equation  (6). 
Thus,  in  this  report  no  attempt  is  made  to  compute  the  absolute  E-field  strength  in  volts/meter 
in  the  elliptical  target.  Instead,  the  peak  power  in  the  elliptical  target  is  adjusted  (by  a  scale 
factor)  to  produce  a  desired  maximum  focal-region  temperature  ( Tmax )  after  t  minutes.  Note:  an 
approximate  absolute  scale  factor  could  be  computed  by  making  an  initial  computer  simulation 
with  an  infinite  homogeneous  water  bolus  and  then  matching  the  target  boundary  field  to  the 
infinite  homogeneous  muscle  tissue  simulation. 

The  computer  simulation  model  is  related,  in  part,  to  the  hyperthermia  annular  phased-array 
antenna  system  shown  in  Figure  1.  This  hyperthermia  system  uses  a  60-cm  array  diameter  with 
eight  uniformly  spaced  dipole  elements  (fed  as  four  pairs)  which  operate  over  the  frequency  band 
60-120  MHz  21  .  The  simulated  array  is  assumed  to  be  fully  adaptive,  and  with  eight  adaptive 
elements  seven  independent  nulls  can  be  formed  while  simultaneously  focusing  on  a  tumor.  For 
comparison,  the  BSD-2000  hyperthermia  system  could  be  used  to  form  up  to  three  independent 
nulls. 

In  this  section,  we  assume  that  the  adaptive  radiation  pattern  null-width  characteristics  in 
a  homogeneous  target  will  be  similar  to  the  characteristics  observed  in  an  inhomogeneous  target. 
The  null  width  is  directly  related  to  the  wavelength  and  there  is  only  a  5  percent  change  in  wave¬ 
length  between  the  assumed  muscle  tissue  and  water  bolus.  With  this  assumption,  the  transmit 
array  is  embedded  in  homogeneous  tissue,  which  allows  direct  use  of  the  thin-wire  moment-method 
formulation  in  Section  2.3.  After  computing  the  two-dimensional  electric-field  distribution  in  the 
homogeneous  medium,  we  then  consider  only  an  elliptical  portion  of  the  homogeneous  region  and 
use  the  ellipse  as  the  homogeneous  target.  In  the  thermal  analysis,  the  elliptical  target  is  sur¬ 
rounded  with  a  constant  10:C  water  bolus.  The  electric  field  amplitude  is  scaled  to  produce  a 
46  C'  peak  temperature,  at  time  t  =  20  minutes,  at  the  center  of  the  elliptical  phantom.  The  initial 
temperature  of  the  phantom  is  assumed  to  be  25°C  (room  temperature). 

All  computer  simulations  presented  in  this  section  are  at  a  120-MHz  operating  frequency 
with  four  auxiliary  nulling  sensors:  that  is,  A'aux  =  4.  Degraded  system  performance  due  to  an 
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insufficient  number  uf  auxiliary  probes  is  considered  in  Appendix  A.  The  parameters  used  in  the 
electrical  and  thermal  analyses  are  summarized  in  Table  1.  These  parameters,  obtained  in  part 


TABLE  1 

Parameters  Used  in  Electrical/ Thermal  Analysis 


1715M-5 


PARAMETER 

PHANTOM  MUSCLE 
TISSUE 

DISTILLED  WATER 

DIELECTRIC  CONSTANT 
@  100  MHz 

73.5 

80.0 

ELECTRICAL  CONDUCTIVITY 
@  100  MHz 

0.5  S/m 

0.0001  S/m 

DENSITY 

970.0  kg/m3 

1000.0  kg/m3 

SPECIFIC  HEAT 

3516.0  J/kg  °C 

4200.0  J/kg  °C 

THERMAL  CONDUCTIVITY 

0.544  W/m  °C 

0.6019  W/m  °C 

from  Sullivan  [9j.  are  for  a  frequency  of  100  MHz:  we  assume  that  similar  values  of  the  parameters 
will  exist  at  120  MHz,  Notice  that  the  relative  dielectric  constants  of  phantom  muscle  tissue  and 
distilled  water  are  very  similar;  however,  the  electrical  conductivities  are  vastly  different.  The 
relevant  thermal  characteristics — density,  specific  heat,  and  thermal  conductivity  [35] — are  veri¬ 
similar  for  phantom  muscle  tissue  and  distilled  water. 

3.1  ELECTRIC  FIELD  FOR  ARRAY  IN  HOMOGENEOUS  TISSUE 

Substituting  the  values  /  =  120  MHz.  a  =  0.5  S/m.  and  eT  —  73.5  into  Equation  (29) 
yields  =  10.0  -r  j23.8  for  the  muscle  tissue.  With  3m  =  23.8  radians/m,  the  wavelength  in 
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the  phantom  muscle  tissue  is  Am  =  26.5  cm.  The  attenuation  constant  for  the  muscle  tissue  is 
am  =  10.0  radians/m.  Similarly,  for  distilled  water  -yw  =  0.0021  +  j22.5,  so  the  wavelength  is 
Au.  =  27.9  cm.  The  attenuation  constant  for  the  distilled  water  medium  is  aw  =  0.0021  radians/m. 
The  propagation  loss  in  the  phantom  muscle  tissue  is  201og10e-10  0,  or  -0.87  dB/cm.  Similarly, 
the  propagation  loss  in  the  distilled  water  is  found  to  be  -0.0002  dB/cm.  Thus,  the  total  loss 
due  to  propagation  through  15  cm  of  distilled  water  is  0.003  dB.  For  15  cm  of  muscle  tissue 
the  corresponding  loss  is  13.1  dB.  The  wave  impedance  in  the  muscle  tissue  is  computed  from 
Equation  (34)  a s  T]m  —  33.9  +j  14.2  17.  and  similarly  in  the  distilled  water  r)u  =  42.1  +j  0.004  17. 

The  geometry  used  in  the  simulations  is  shown  in  Figure  13.  A  60-cm-diameter  ring  array 
of  eight  dipoles  uniformly  surrounds  a  fictitious  elliptical  target  zone  with  major  axis  30  cm  and 
minor  axis  20  cm.  The  length  of  each  perfectly  conducting  center- fed  dipole  array  element  at 
120  MHz  in  the  infinite  homogeneous  muscle  tissue  is  A/2,  or  13.25  cm.  The  array  focus  is  assumed 
at  the  origin,  and  four  auxiliary  short-dipole  sensors  with  length  1.27  cm  (0.05A)  are  positioned  in 
(x.y.z)  coordinates  at  (15  cm.  0.  0).  (-15  cm.  0,  0),  (0,  0,  10  cm),  and  (0,  0,  -10  cm);  that  is, 
the  auxiliary  E-field  sensors  are  located  every  90°  in  azimuth  on  the  perimeter  of  the  target.  In 
rectangular  coordinates,  each  dipole  is  oriented  in  the  y  direction  and  the  feed  terminals  of  each 
dipole  are  located  at  y  =  0.  The  moment-method  computer  simulations  were  run  on  a  Sun  3/260 
workstation.  The  total  CPU  time  for  a  complete  moment-method  run  is  19.2  minutes.  This  CPU 
time  includes  computing  the  quiescent  and  adaptive  radiation  patterns  on  a  41  by  41  grid  of  points. 
The  CPU  time  without  radiation  pattern  calculations  is  33  seconds. 

The  two-dimensional  radiation  pattern  in  the  plane  y  =  0,  before  nulling,  at  120  MHz  with 
uniform  amplitude  and  phase  illumination,  is  shown  in  Figure  14.  The  calculated  data  are  collected 
on  a  41  by  41  grid  of  points  over  a  square  region,  with  side  length  76.2  cm,  centered  at  the  focus. 
The  spacing  between  data  points  is  1.905  cm.  or  0.072 A,  and  the  contour  levels  are  displayed  in  10- 
dB  steps.  The  E-field  data  are  computed  for  the  case  of  a  1.27-cm  short-dipole  observation  probe. 
The  positions  of  the  eight  dipole  radiators  are  clearly  evident  by  the  — 20-dB  contours  surrounding 
each  element.  The  radiation  pattern  is  symmetric  because  of  the  symmetry  of  the  array  and  the 
assumed  homogeneous  medium. 

Finer  contour  levels  (1-dB  steps)  for  the  quiescent  radiation  pattern  are  displayed  in  Figure  15. 
Here,  it  is  evident  that  the  focused  main  beam  of  the  ring  array  is  increasing  in  amplitude  as 
the  observation  point  moves  closer  to  the  focus.  Away  from  the  main  beam  region,  the  pattern 
amplitude  is  seen  to  increase  as  the  observation  position  moves  away  from  the  focus. 

A  quiescent  radiation  pattern  cut  at  z  =  0  is  shown  in  Figure  16.  The  large  amplitude  that 
occurs  at  ±30  cm  is  due  to  the  E-field  probe's  close  proximity  to  the  transmitting  elements.  The 
large  attenuation  that  occurs  from  the  ar.av  diameter  to  the  focus  is  due  to  the  1/r  attenuation 
loss  and  the  loss  in  the  uniform  homogeneous  muscle  tissue.  The  radiation  pattern  cut  for  x  =  0 
in  Figure  17  is  identical  to  the  pattern  in  Figure  16  due  to  the  symmetry  of  the  array.  In  both 
Figures  16  and  17  the  boundary  of  the  fictitious  elliptical  target  zone  is  indicated.  The  increasing 
radiation  pattern  amplitude  near  the  left  and  right  sides  of  the  elliptical  target  (Figure  16)  will  be 
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Figure  13.  Geometry  for  eight-element  ring  array  and  four  E-field  auxiliary  sensors. 
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Fig ure  14.  Simulated  two-dimensional  quiescent  radiation  pattern  at  120  MHz  for  uni¬ 
formly  illuminated  eight-element  nng  array  in  infinite  homogeneous  conducting  medium 
(phantom  muscle  tissue:  eT  =  73.5 .0  =  0.5/  A  fictitious  elliptical  target  region  is  indi¬ 
cated  by  the  dashed  curve.  Radiation  contour  levels  are  in  10-dB  steps. 


31 


E-FIELD  BEFORE  NULLING 


-20  0  20 

x  (cm) 


Figure  15.  Simulated  two-dimensional  quiescent  radiation  pattern  at  120  MHz  for  uni¬ 
formly  illuminated  eight-element  ring  array  in  infinite  homogeneous  conducting  medium 
(phantom  muscle  tissue:  er  =  73.5 .o  =  0.5 ).  A  fictitious  elliptical  target  region  is  indi¬ 
cated  by  the  dashed  curve.  Radiation  contour  levels  are  given  in  1-dB  steps. 
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shown  to  produce  hot  spots  in  the  thermal  distribution.  As  shown  in  the  next  section,  because 
the  top  (anterior)  and  bottom  (posterior)  of  the  elliptical  target  are  not  as  strongly  illuminated  as 
the  left  and  right  sides,  no  quiescent  hot  spots  occur  at  the  top  or  bottom.  From  Figure  16,  the 
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Figure  16.  Simulated  one- dimensional  quiescent  radiation  pattern  cut  (z  —  0)  at 
120  MHz  for  uniformly  illuminated  eight-element  ring  array  m  infinite  homogeneous  con¬ 
ducting  medium  ( phantom  muscle  tissue:  er  =  73.5, er  =  0.5/ 


ring-array  half-power  beamwidth  in  the  target  region  is  approximately  13  cm,  or  approximately 
one-half  the  wavelength  (26.5  cm)  in  the  phantom, muscle  tissue.  2  The  adaptive  nulling  resolution 


2Sullivan  [54]  has  recently  published  simulations  and  measurements  of  SAR  patterns  at  70-110 
MHz.  using  a  homogeneous  35  x  25-cm  elliptical  target  surrounded  with  a  layer  of  fat  and  a  water 
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Figure  17.  Simulated  one- dimensional  quiescent  radiation  pattern  cut  (x  =  0)  at 
120  MHz  for  uniformly  illuminated  eight-element  ring  array  in  infinite  homogeneous  con¬ 
ducting  medium  (phantom  muscle  tissue:  er  =  73.5,  a  =  0.5/ 


or  closest  allowed  spacing  between  a  deep  adaptive  null  and  the  main  beam  is  equal  to  the  half¬ 
power  beamwidth  of  the  antenna  [55] .  Thus,  the  closest  allowed  null  position  is  13  cm  from  the 
focus.  Since  the  target  width  is  30  cm,  two  nulls  can  be  formed  at  ( x  =  ±15  cm,  z  =  0)  without 
disturbing  the  focus.  However,  if  two  deep  nulls  are  formed  at  (x  =  0,  z  =  ±10  cm)  the  focus  will 
be  compromised.  In  practice,  the  water  bolus  would  restrict  the  placement  of  short-dipole  sensors 
to  the  surface  of  the  target.  Thus,  only  weak  nulls  can  be  formed  at  (x  =  0,  z  =  ±10  cm)  so  that 
the  focus  will  not  be  affected  by  the  adaptive  nulling  process. 

Next,  adaptive  radiation  patterns  are  computed  with  four  auxiliary  dipole  sensors;  their 
positions  are  shown  in  Figure  13.  The  value  of  the  receiving  gain  for  auxiliary  dipoles  1  and  2  is 
adjusted  to  produce  a  greater-than-35-dB  SNR.  This  amount  of  SNR  results  in  greater  than  35  dB 
of  nulling  in  the  direction  of  auxiliary  dipoles  1  and  2.  In  contrast,  the  gain  values  for  auxiliary 
dipoles  3  and  4  are  turned  down  to  produce  about  a  3-dB  SNR.  Thus,  only  about  3  dB  of  nulling 
will  occur  at  sensor  positions  3  and  4  as  the  adaptive  algorithm  reduces  the  interference  to  the  noise 
level  of  the  receiver.  The  reason  for  choosing  these  null  depths  will  become  apparent  with  the  data 
that  follow.  Figure  18  shows  the  two-dimensional  radiation  pattern  after  nulling  with  four  auxiliary 
sensors.  Two  deep  adaptive  nulls  at  x  =  ±15  cm  occur  as  expected,  and  weak  nulling  occurs  at 
z  =  ±10  cm,  also  as  expected.  The  two  deep  nulls  in  the  2  =  0  cut  are  quantified  in  Figure  19, 
where  greater  than  35  dB  of  interference  nulling  or  pattern  reduction  occurs  at  x  =  ±15  cm.  The 
peak  level  at  the  focus  is  adjusted  to  0  dB  for  both  the  quiescent  and  adaptive  patterns.  Two 
weak  adaptive  nulls  are  in  effect  in  the  x  =  0  radiation  pattern  cut  shown  in  Figure  20;  however, 
weak  nulls  are  desired  in  this  cut  due  to  temperature  considerations,  as  shown  in  the  next  section. 
The  weak  nulls  in  effect  in  the  adaptive  patterns  reduce  variation  from  the  quiescent  radiation 
pattern.  Finally,  the  transmit  array  weights  before  and  after  nulling  and  the  covariance  matrix 
eigenvalues  are  shown  in  Figures  21  and  22,  respectively.  The  adaptive  transmit  weights  exhibit 
a  5-dB  dynamic  range  in  Figure  21(a).  There  are  two  large  eigenvalues  and  two  weak  (nonzero) 
eigenvalues  shown  in  Figure  22.  These  eigenvalues  are  directly  associated  with  the  two  high-SNR 
auxiliary  sensors  and  the  two  weak-SNR  auxiliary  sensors.  Note  that  the  0-dB  level  in  Figure  22  is 
equal  to  the  noise  level.  The  probe-array  output  power  before  and  after  adaptive  nulling  is  31.4  dB 
and  0.9  dB.  respectively.  This  difference  in  power  before  and  after  nulling  means  that  the  adaptive 
cancellation  is  -30.5  dB.  The  output  file  for  the  moment-method  simulation  of  this  section  is  given 
in  Appendix  B. 


bolus.  A  comparison  of  Sullivan's  110-MHz  SAR  simulation  (Figure  5(a)  of  his  paper)  and  the 
results  (at  120  MHz)  presented  in  Figures  16  and  17  of  this  report  indicate  reasonable  agreement 
over  the  elliptical  target  region. 
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Figure  18.  Simulated  two-dimensional  adaptive  radiation  pattern  at  120  MHz  for  eight- 
element  ring  array  in  infinite  homogeneous  conducting  medium  (phantom  muscle  tissue: 
tr  =  73.5.  a  =  0  5  j.  Radiation  contour  levels  are  given  in  10-dB  steps.  Four  auxiliary 
sensors  are  used  in  forming  the  adaptive  pattern.  The  quiescent  focus  is  at  (0,0). 
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Figure  19  Simulated  one-dimensional  quiescent  and  adaptive  radiation  patterns  in  the 
2=0  cut  at  120  MHz  for  eight-element  ring  array  in  infinite  homogeneous  conducting 
medium  (phantom  muscle  tissue:  er  =  73.5,  a  =  0.5,).  Four  auxiliary  sensors  are  used  in 
forming  the  adaptive  pattern. 
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Figure  20.  Simulated  one-dimensional  quiescent  and  adaptive  radiation  patterns  in  the 
x  —  0  cut  at  120  MHz  for  eight-element  nng  array  in  infinite  homogeneous  conducting 
medium  (phantom  muscle  tissue:  er  =  73.5,  o  =  0.5 ).  Four  auxiliary  sensors  are  used  in 
forming  the  adaptive  pattern. 
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Figure  21.  Transmit  array  weights  before  and  after  adaptive  nulling.  A  dynamic  range 
of  about  5  dB  is  evident  for  the  adaptive  weights.  Four  auxiliary  sensors  are  used  in  the 
adaptive  process,  (a)  Amplitude:  (b)  Phase. 
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Figure  22.  Channel  covariance  matrix  eigenvalues  (degrees  of  freedom)  used  in  the  adap¬ 
tive  process  with  four  auxiliary  sensors. 


3.2  TEMPERATURE  DISTRIBUTION  IN  ELLIPTICAL  PHANTOM 

In  this  section,  the  transient  thermal  analysis  software  is  used  to  compute  the  temperature 
distribution  in  an  elliptical  phantom  surrounded  with  a  constant-temperature  water  bolus.  The 
41  by  41  two-dimensional  E-field  radiation  pattern  data  from  the  previous  section  are  used  as  the 
power  source  for  the  thermal  node  network.  Two  spacings  of  nodes  are  considered  here:  initially, 
the  node  spacing  At  =  Az  =  Al  =  1.905  cm  (coarse  grid)  is  used;  to  check  convergence,  the  node 
spacing  is  then  decreased  by  a  factor  of  two  to  Al  =  0.9525  cm  (fine  grid).  The  coarser  spacing  is 
shown  to  be  adequate. 

3.2.1  1.905-cm  Thermal  Node  Spacing 

Note:  the  scale  factors  used  to  convert  the  normalized  E- field  distributions  to  a  power  level 
that  induces  a  46°C  peak  temperature  at  t  =  20  minutes  are  94.1  dB  and  96.0  dB  for  the  quiescent 
and  adaptive  patterns,  respectively.  These  scale  factors  are  determined  by  trial  and  error.  From 
equations  given  in  Section  2.3.2  and  the  parameter  values  in  Table  1,  all  resistors  in  the  phantom 
muscle  tissue  had  a  value  of  96.5°C/VV  and  all  resistors  in  the  water  bolus  had  a  value  of  87.2°C/W. 
The  value  of  the  capacitors  in  the  phantom  muscle  tissue  is  23.6  J/°C;  values  for  capacitors  in  the 
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water-bolus  region  are  not  used  in  the  input  to  the  transient  thermal  analysis  software.  Instead,  a 
constant  temperature  of  10°C  is  enforced  at  each  water-bolus  node  With  a  41  by  41  grid,  a  total 
of  3280  resistors  and  1681  capacitors  are  used  in  the  thermal  simulation.  The  CPU  time  required 
to  compute  the  temperature  distribution  on  a  41  by  41  grid  is  under  four  minutes. 

Figure  23  shows  the  two-dimensional  temperature  distribution  produced  at  time  t  =  20  min¬ 
utes  in  the  elliptical  phantom  muscle  tissue  target  without  adaptive  nulling.  To  geneiate  Figure  23, 
the  power  source  used  in  the  transient  thermal  analysis  is  the  quiescent  radiation  pattern  given  in 
Figure  14.  The  initial  temperature  (at  time  t  =  0)  is  25°C.  Notice  the  occurrence  of  two  hot  spots 
on  the  left  and  right  sides  of  the  elliptical  phantom.  The  peak  temperature  on-focus  is  46°C,  which 
is  achieved  by  scaling  the  normalized  quiescent  E-field  as  described  earlier.  The  two  hot  spots  are 
quantified  in  the  r  =  0  temperature  pattern  cut  shown  in  Figure  24.  While  the  peak  temperature 
at  each  hot  spot  in  Figure  24  is  only  41°C,  the  temperature  profile  for  x  =  0  in  Figure  25  shows  no 
hot  spots.  As  any  undesired  hot  spot  is  a  potential  source  for  compromising  the  therapy  session, 
adaptive  nulling  is  used  to  reduce  the  sidelobes  corresponding  to  the  hot  spots. 

Figure  26  shows  the  simulated  two-dimensional  thermal  distribution  at  time  t  =  20  minutes, 
with  adaptive  nulling  at  four  auxiliary  sensors  in  effect.  The  focal-spot  diameter  with  adaptive 
nulling  is  equivalent  to  the  focal-spot  diameter  before  adaptive  nulling.  Hot  spots  on  the  left  and 
right  sides  of  the  target  are  eliminated.  A  comparison  of  the  temperature  distribution  before  and 
after  nulling  along  the  major  axis  [z  =  0)  of  the  target  ellipse  is  made  in  Figure  27.  Similarly,  the 
temperature  distribution  before  and  after  nulling  along  the  minor  axis  of  the  target  ellipse  is  shown 
in  Figure  28. 


3.2.2  Convergence  Check:  0.9525-cm  Thermal  Node  Spacing 

The  convergence  of  the  previous  thermal  simulations  was  verified  by  increasing  the  density 
of  E-field  observation  probe  positions  by  a  factor  of  two,  with  a  new  spacing  between  points  of 
0.9525  cm.  still  with  a  41  by  41  grid.  The  ring  array  operates  as  before  at  120  MHz,  and  there  are 
four  auxiliary  sensors  laid  out  as  shown  in  Figure  13.  As  the  auxiliary  positions  are  the  same,  the 
adaptive  weights  and  covariance  matrix  eigenvalues  in  Figures  21  and  22.  respectively,  remain  the 
same.  From  the  parameter  values  in  Table  1.  all  resistors  in  the  finer-grid  muscle-tissue  phantom 
had  a  value  of  193.0  C 'W  and  all  resistors  in  the  water  bolus  had  a  value  of  174.4CC/W.  The 
value  of  the  capacitors  in  the  phantom  muscle-tissue  is  2.95  J/rC.  Again,  a  constant  temperature 
of  10CC  is  enforced  at  each  water-bolus  node.  The  E-field  scaling  factors  to  raise  the  focal-point 
temperature  to  46CC  before  and  after  nulling  are  76.5  dB  and  78.4  dB.  respectively.  The  finer-grid 
two-dimensional  thermal  distributions  before  and  after  nulling  are  showm  in  Figures  29  and  30, 
respectively.  Although  the  temperature  contours  are  smoother,  the  general  agreement  between 
these  patterns  and  the  coarser-grid  patterns  in  Figures  23  and  26  are  evident.  Similarly,  one¬ 
dimensional  thermal  pattern  cuts  with  the  finer  grid  are  shown  in  Figures  31  (x  axis)  and  32  ( z  axis); 
good  agreement  with  the  coarse-grid  patterns  (Figures  27  and  28)  is  observed.  In  particular,  the 
finer  detail  in  Figure  31  shows  that  the  hot  spots  are  at  a  42CC  level  compared  to  the  41°C  level 
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observed  for  the  coarse  grid  (Figure  24).  Thus,  convergence  of  the  coarse-grid  thermal  patterns  is 
demonstrated. 


172660-23 


Figure  23.  Simulated  two-dimensional  thermal  pattern  at  time  t  =  20  minutes  be¬ 
fore  nulling  in  elliptical  phantom  muscle-tissue  target  surrounded  with  1CF  C  constant- 
temperature  water  bolus.  The  incident  RF  power  distribution,  from  Figure  14,  is  at 
120  MHz.  Temperature  contour  levels  are  given  in  2C  C  steps.  Hot  spots  on  the  left 
and  right  sides  of  the  target  are  observed. 
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Figure  24 •  Simulated  one-dimensional  (z  =  Oj  thermal  pattern  at  time  t  =  20  minutes 
before  nulling  in  elliptical  phantom  muscle-tissue  target  surrounded  with  1CFC  constant- 
temperature  water  bolus.  The  quiescent  incident  RF  power  distribution,  from  Figure  14 , 
is  at  120  MHz.  Hoi  spots  on  the  left  and  right  sides  of  the  target  are  observed. 
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LINEAR  CUT  OF  z  AXIS  @  x  =  0  (cm) 

Figure  25.  Simulated  one-dimensional  (x  —  0^  thermal  pattern  at  time  t  =  20  minutes 
before  nulling  in  elliptical  phantom  muscle-tissue  target  surrounded  with  1CFC  constant- 
temperature  water  bolus.  The  quiescent  incident  RF  power  distribution,  from  Figure  14, 
is  at  120  MHz.  No  undesired  hot  spots  are  present. 
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Figure  26.  Simulated  two-dimensional  thermal  pattern  at  time  t  =  20  minutes  (with 
adaptive  nulling  at  four  auxiliary  sensors  in  effect)  in  elliptical  phantom  muscle-tissue 
target  surrounded  with  l(f  C  constant-temperature  water  bolus.  The  adapted  incident  RF 
power  distribution,  from  Figure  18.  is  at  120  MHz.  Temperature  contour  levels  are  given 
in  2C  C  steps.  Hot  spots  on  the  left  and  right  sides  of  the  target  are  eliminated. 
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Figure  27.  Simulated  one -dimensional  (z  —  Q)  thermal  patterns  at  time  t  =  20  minutes 
before  and  after  nulling  in  elliptical  phantom  muscle-tissue  target  surrounded  with  1(F  C 
constant-temperature  water  bolus.  The  quiescent  incident  RF  power  distribution,  from 
Figure  14 ■  is  at  120  MHz.  Hot  spots  on  the  left  and  right  sides  of  the  target  are  clearly 
eliminated  by  the  adaptive  nulling  process.  Four  auxiliary  sensors  are  used  in  the  adaptive 
process. 
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LINEAR  CUT  OF  z  AXIS  @  x  =  0  (cm) 


Figure  28.  Simulated  one- dimensional  (x  =  0)  thermal  patterns  at  time  t  —  20  minutes 
before  and  after  nulling  in  elliptical  phantom  muscle-tissue  target  surrounded  with  1(F  C 
constant-temperature  water  bolus.  The  quiescent  incident  RF  power  distribution,  from 
Figure  14.  is  at  120  MHz.  No  undesired  hot  spots  are  present.  Four  auxiliary  sensors  are 
used  in  the  adaptive  process. 


47 


172666-29 


Figure  29.  E-field  probe-sample-spacing  convergence  check  for  simulated  two-dimensional 
thermal  pattern  at  time  t  =  20  minutes  ( before  adaptive  nulling)  in  elliptical  phantom 
muscle-tissue  target  surrounded  with  1(F  C  constant-temperature  water  bolus.  The  quies¬ 
cent  incident  RF  power  distribution  is  at  120  MHz.  Temperature  contour  levels  are  given 
in  2°  C  steps.  The  grid  spacing  is  one-half  the  spacing  of  that  in  Figure  23.  Hot  spots 
on  the  left  and  right  sides  of  the  target  are  present,  as  previously  observed  for  the  coarser 
probe-sample  spacing. 
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Figure  30.  E-field  probe-sample-spacing  convergence  check  for  simulated  two-dimensional 
thermal  pattern  at  time  t  =  20  minutes  (with  adaptive  nulling  in  effect)  in  elliptical 
phantom  muscle-tissue  target  surrounded  with  ltFC  constant -temperature  water  bolus. 
The  adapted  incident  RF  power  distribution  is  at  120  MHz.  Temperature  contour  levels 
are  given  in  2~  C  steps.  The  grid  spacing  is  one-half  the  spacing  of  that  in  Figure  26.  Hot 
spots  on  the  left  and  right  sides  of  the  target  are  eliminated,  as  previously  observed  for  the 
coarser  probe-sample  spacing. 
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Figure  31.  E-field  probe-sample-spacing  convergence  check  for  simulated  one -dimensional 
(z  =  0)  thermal  patterns  at  time  t  =  20  minutes  before  and  after  nulling  in  elliptical  phan¬ 
tom  muscle-tissue  target  surrounded  with  Iff1  C  constant-temperature  water  bolus.  The 
quiescent  incident  RF  power  distribution  is  at  120  MHz  with  0.9525-cm  sample  spacing. 
Hot  spots  on  the  left  and  right  sides  of  the  target  are  eliminated  by  the  adaptive  nulling  pro¬ 
cess.  as  previously  observed  in  Figure  27.  Four  auxiliary  sensors  are  used  in  the  adaptive 
process. 
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Figure  32.  E-field  probe- sample-spacing  convergence  check  for  simulated  one- dimensional 
(x  =  0)  thermal  patterns  at  time  t  =  20  minutes  before  and  after  nulling  in  elliptical  phan¬ 
tom  muscle-tissue  target  surrounded  with  IGF  C  constant-temperature  water  bolus.  The 
quiescent  incident  RF  power  distribution  is  at  120  MHz  with  0.9525-cm  sample  spac¬ 
ing.  A’o  undesired  hot  spots  are  present.  Four  auxiliary  sensors  are  used  in  the  adaptive 
process. 
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4.  CONCLUSION 


Focused  near-field  adaptive  nulling  has  been  applied  to  the  problem  of  generating  therapeutic 
temperature  distributions,  free  of  undesired  hot  spots,  for  deep-regional  cancer  treatment  in  a  target 
body.  A  noninvasive  adaptive  hyperthermia  system  concept  has  been  described  and  is  analyzed  in 
detail.  Auxiliary  short-dipole  nulling  sensors  are  positioned  on  the  surface  of  the  target.  Due  to 
the  finite  width  of  an  electric-field  null,  a  null  formed  at  the  target  surface  extends  into  the  interior 
region  of  the  target.  When  the  width  of  the  null  is  properly  chosen,  undesired  hot  spots  inside 
the  target  are  eliminated  while  simultaneously  focusing  on  a  deep-seated  tumor.  The  resolution 
between  a  deep  null  and  focus  is  limited  by  the  half-power  beamwidth  of  the  hyperthermia  array. 
The  resolution  can  be  somewhat  enhanced  by  using  weak  nulls  whenever  the  spacing  between  the 
null  and  focus  is  less  than  one  half-power  beamwidth.  The  depth  of  null  is  controlled  by  the  SNR 
at  the  auxiliary  probe  position. 

A  theory  for  analyzing  adaptive  nulling  with  a  phased  array  in  an  infinite  homogeneous 
conducting  medium  has  been  developed.  The  method  of  moments  is  used  to  compute  the  electric 
field  at  a  short-dipole  sensor  due  to  a  thin-wire  dipole  ring  array.  The  SMI  algorithm  is  used  to 
adaptively  control  the  transmit  array  weights  and  to  form  radiation  pattern  nulls.  The  method  used 
for  analyzing  the  transient  thermal  behavior  of  an  RF-illuminated  target  by  using  an  equivalent 
electric  circuit  network  has  been  described. 

Computer  simulations  of  the  electric  field  and  induced  temperature  distribution  for  an  eight- 
element  dipole  ring  array  at  120  MHz  have  been  presented.  Four  auxiliary  E- field  short-dipole 
sensors  are  used  in  eliminating  two  widely  separated  hot  spots  surrounding  the  array  focus  at  the 
center  of  an  elliptical  target  with  dimensions  30  cm  by  20  cm.  Calculated  data  indicate  that  si¬ 
multaneous  nulling  and  focusing  can  yield  desirable  electric-field  distributions  for  hyperthermia. 
Undesired  hot  spots  controlled  with  an  insufficient  number  of  auxiliary  nulling  probes  have  been 
shown  to  simply  redistribute  within  the  target.  Thus,  an  adequate  number  of  auxiliary  probes 
must  be  provided  to  ensure  that  no  undesired  hot  spots  are  formed  in  deep-regional  hyperther¬ 
mia  therapy.  Although  the  computer  simulations  are  implemented  for  a  phantom  target,  adaptive 
nulling  may  prove  beneficial  in  clinical  trials  of  hyperthermia  with  living  targets.  Further  adaptive- 
nulling  simulations  with  inhomogeneous  targets  are  desirable,  as  are  experimental  measurements  of 
adaptive  hyperthermia  in  phantom  targets.  A  study  of  proper  number  and  positioning  of  the  non¬ 
invasive  sensors  for  treating  arbitrary  tumor  positions  needs  to  be  performed.  The  adaptive-nulling 
technique  described  here  should  be  applicable  to  both  invasive  and  noninvasive  RF  hyperthermia 
systems,  as  well  as  ultrasound  hyperthermia  systems. 
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APPENDIX  A 

SYSTEM  DEGRADATION  DUE  TO  INSUFFICIENT  NUMBER  OF 

AUXILIARY  PROBES 


This  section  shows  simulation  results  which  indicate  that  an  insufficient  number  of  auxiliary 
nulling  probes  will  degrade  adaptive  hyperthermia  system  performance.  The  simulation  parameters 
are  the  same  as  presented  in  Section  3.  except  only  two  auxiliary  probes  ( Naux  =  2)  are  used  to 
form  nulls  at  the  surface  of  the  elliptical  muscle-tissue  target.  The  geometry  used  in  the  simulations 
is  shown  in  Figure  A-l.  A  60-cm  diameter  ring  array  of  eight  dipoles  uniformly  surround  a  fictitious 
elliptical  target  zone  with  major  axis  30  cm  and  minor  axis  20  cm.  The  length  of  each  perfectly 
conducting  center- fed  dipole  array  element,  at  120  MHz  in  the  infinite  homogeneous  muscle  tissue, 
is  A/2  or  13.25  cm.  The  array  focus  is  assumed  to  be  at  the  origin,  and  two  auxiliary  short- 
dipole  sensors  with  length  1.27  cm  are  positioned  at  (15  cm,  0)  and  (-15  cm,  0).  In  rectangular 
coordinates,  each  dipole  is  oriented  in  the  y  direction  and  the  feed  terminals  of  each  dipole  are 
located  at  y  =  0.  The  quiescent  radiation  patterns  and  thermal  distributions  are  the  same  as 
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Figure  A-l.  Geometry  for  eight-element  nng  array  and  two  E-field  auxiliary  sensors. 
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presented  in  Figures  14,  15,  16,  17,  23,  24,  and  25. 

With  only  two  auxiliary  sensors,  the  two-dimensional  adapted  electric-field  radiation  pat¬ 
tern  shown  in  Figure  A-2  indicates  that  both  auxiliary  positions  are  nulled.  However,  the  two- 
dimensional  adaptive  thermal  distribution  in  Figure  A-3  shows  that  the  quiescent  hot  spots  (see 
Figure  23)  have  only  been  redistributed.  The  one-dimensional  thermal  distribution  in  Figure  A-4 
shows  that  the  hot  spots  are  reduced  on  the  left  and  right  sides  of  the  elliptical  target,  and  Fig¬ 
ure  A-5  shows  that  the  hot  spots  are  redistributed  to  lie  along  the  z  axis.  Thus,  two  additional 
auxiliary  probes  (probes  3  and  4  in  Figure  13)  are  necessary  to  avoid  the  possibility  of  redistributed 
hot  spots.  Finally,  the  adaptive  weights  and  covariance  matrix  eigenvalues  are  shown  in  Figures  A-6 
and  A-7,  respectively.  There  are  only  small  differences  between  the  weights  and  eigenvalues  com¬ 
pared  to  the  four-auxiliary  probe  case  (see  Figures  21  and  22).  The  transmit  weights  cover  a  7-dB 
dynamic  range  and  there  are  only  two  nonzero  eigenvalues.  The  probe-array  power  outputs  before 
and  after  adaptive  nulling  are  31.4  dB  and  0.0  dB,  respectively.  Thus,  the  interference  cancellation 
is  to  the  noise  level,  or  -31.4  dB.  An  output  file  for  the  moment-method  simulation  for  this  section 
is  given  in  Appendix  B. 
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E-FIELD  AFTER  NULLING;  2  AUXILIARIES 


Figure  A-2.  Simulated  two-dimensional  adaptive  radiation  pattern  at  120  MHz  for  eight- 
elemen .  ring  array  in  infinite  homogeneous  conducting  medium  ( phantom  muscle  tissue : 
er  =  73.5,  o  =  0.5/.  Radiation  contour  levels  are  given  in  10-dB  steps.  Two  auxiliary 
sensors  are  used  in  forming  the  adaptive  pattern.  The  quiescent  focus  is  at  (0,0). 
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Figure  A-3.  Simulated  two-dimensional  thermal  pattern  at  time  t  —  20  minutes  ( with 
adaptive  nulling  at  only  two  auxiliary  sensors  in  effect)  in  elliptical  phantom  muscle-tissue 
target  surrounded  with  1CF  C  constant-temperature  water  bolus.  The  adapted  incident  RF 
power  distribution,  from  Figure  18,  is  at  120  MHz.  Temperature  contour  levels  are  given 
in  Z~  C  steps.  Initial  hot  spots  ( before  nulling)  on  the  left  and  right  sides  of  the  target  (see 
Figure  23)  are  redistributed  to  the  top  and  bottom  ( anterior  and  posterior  positions). 
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Figure  A-f.  Simulated  one -dimensional  (z  =  0J  thermal  patterns  at  time  t  =  20  minutes 
before  and  after  nulling  in  elliptical  phantom  muscle-tissue  target  surrounded  with  Iff-  C 
constant -temperature  water  bolus.  The  quiescent  incident  RF  power  distribution,  from 
Figure  14.  is  at  120  MHz.  Hot  spots  on  the  left  and  right  sides  of  the  target  are  eliminated 
by  the  adaptive  nulling  process.  Two  auxiliary  sensors  are  used  in  the  adaptive  process. 
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Figure  A- 5.  Simulated  one- dimensional  (x  =  0)  thermal  patterns  at  time  t  =  20  minutes 
before  and  after  nulling  in  elliptical  phantom  muscle-tissue  target  surrounded  with  10PC 
constant-temperature  water  bolus.  The  quiescent  incident  RF  power  distribution,  from 
Figure  14 ■  is  at  120  MHz.  Two  undesired  hot  spots  are  present  in  the  adaptive  pattern. 
Two  auxiliary  sensors  are  used  in  the  adaptive  process. 
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Figure  A-6  Transmit  array  weights  before  and  after  adaptive  nulling.  A  dynamic  range 
of  about  7  dB  is  evident  for  the  adaptive  weights.  Two  auxiliary  sensors  are  used  in  the 
adaptive  process,  (a)  Amplitude;  (b)  Phase. 
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Figure  .4-7.  Channel  covariance  matrix  eigenvalues  (degrees  of  freedom)  used  in  the 
adaptive  process  with  two  auxiliary  sensors. 
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APPENDIX  B 

SOFTWARE  DOCUMENTATION 


This  appendix  lists  sample  input  and  output  files  for  the  adaptive  hyperthermia  simulation. 
Listings  of  the  key  FORTRAN  programs  that  implement  adaptive  nulling  in  a  homogeneous  con¬ 
ducting  medium  are  also  given.  Additional  moment-method  software  not  listed  here  is  contained 
in  Richmond  [37j. 

The  first  data  file  (four  auxiliary  sensors)  is  used  to  generate  the  El- field  results  presented 
in  Section  3.1.  The  second  data  file  (two  auxiliary  sensors)  was  used  in  calculating  the  electric- 
field  distribution  presented  in  Appendix  A.  The  corresponding  output  files  give  the  values  for  the 
array  mutual  coupling,  quiescent  and  adaptive  transmit  weights,  covariance  matrix,  eigenvalues, 
and  cancellation. 


•••Input  data  file  for  adaptive  nulling  with  four  auxiliary  probes, 
•••sdipiamhyper  data  file,  filename  adipjamhyper .datacirconicrr4 
ADIPOLE  NC0LX*8 ,NR0VY*1 ,HZIN-0 .001 ,HLIN*2 . 6 , 

DXIN*6 . 888 ,  DYIN*6 . 888 .  ARADIN*0 . 0039 , 

ICIRC-1 .CRADIN-ll . 81 ,HLSIN*0 . 2S . 

IVL-0 , FCHZ-120 . 0E6 , BWFHZ*1 . OEO .NFREQ-5 , IVR-0 . 

NPDX-0 ,  NXDUM-0 ,  NYDUM-0 , ER3*73 . 5 , SIG3-0 .  5 ,  TD3—  1 . 0 , 

ZL0AD-50 . 0 , ZCHAR*0 . 0 , NGEN=8 . IGEN*0 , 

NSCANS-1, 

THSINC*5 .0, 

IHUT*1 , IBLTSL*0, 

IEN0RM-1,ICHEB-0,SLLDB=20. ,EDGTDB-0. , 

IPATRN-1 , IPRC0M=1 , IANG1P=0 , 

NPHCT-0,NTHPT-499, 

THDR-180. ,THDMIN=-90. , 

NC0LXN*121 , NR0VYN* 1 , NCQLZN* 1 , 

RJLSYIN-0 . 0 ,RLSZIN*0 . 0 , 

NC0LXN=1 17 , NR0VYN* 1 , NC0LZN= 1 , 

INEAR* 1, 

IP0l*2 , IGRNDP*0 , 

ITITPR-O , ITLT0P-O , 

NFC0LX-1 ,NFR0VY=1 , 

IANTX-1 , IANTY-0 ,NP0VER*0 , 

IPC0NN*1 , IPCONF-O , IPCU7F*0 , IPCFX-0 , IPCFY-0 , IPCFZ*0 , 

ITEK-0 , 

IQUAN*0 , IRNERR*0 , ELERDB-0 . 02 , ELERDG-0 . 2 , NBM0D* 12 , 

NRAN-1  ,NBADWT*32 .  AtfERDB-0 . 0 .  AVERDG*0 . 0 , 

NAUX*7 ,IAUXA(1)*1 ,2, 3, 4, 5, 6, 7, 

IATTEN*1 , 

N  JA."S-7 ,  ISLC-0 ,  AUXADB  (1 )  -8*0 . 0 , 

PVRJDB ( 1 ) *40 .  ,40.  ,15.  ,15.  ,3*-99.0, 

YJAMIN (1 )*7*0 -  0 , 

XF0CIN*0 . 0 , ZFOCIN*0 . 0 , 

XJAHIN(l)*5.9,-5.9,2*0.0,3*-0.0, 

ZJ AMIN(l)-0. 0,0. 0,4.0, -4. C, 3*0.0 
YNIN-0.0, 

NC0LXN*101 ,NTHPT*101 , 

XNIN--1E. ,ZNIN*-15 . ,RLSXIN*30. .RLSZIN-30. .NCOLXN-41 .NCOLZN-41 , 

A  EN  D 
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•••Output  data  file  for  adaptive  nulling  with  four  auxiliary  probes . 
•••sdipjaahypar  data  file,  filename  sdipjamhyper.datacirconierr4 
XJAMIN(l) , YJAMIN(l) .ZJAMIN(l)-  6.900  0.000  0.000 

FGHZ-  0.1200000 

DX.DY.HL.ARAD-  0.17496  0.17496  0.06604  0.00010 

IV, VTA-  1  0.6916E-02  0.2211E-02 

IV. VTA-  2  0.5916E-02  0.2211E-02 

IV, VTA-  3  0.5916E-02  0.2211E-02 

IV, VTA-  4  0 . 6916E-02  0.2211E-02 

IV, VTA-  5  0.6916E-02  0.2211E-02 

IV, VTA-  6  0.5916E-02  0.2211E-02 

IV, VTA-  7  0.S916E-02  0.2211E-02 

IV, VTA-  8  0.6916E-02  0.22UE-02 

•••••••••array  nutual  impedance  matrix  (firet  rov)»«»»»*» 

2(1,  1)»  0.91706E+02  0.12403E+02 

2(1,  2)— 0.33552E-01  0.19136E+00 

ZCl,  3)— 0.14488E-01-0.64125E-02 

Z(l,  4)-  0.29420E-02  0.13780E-02 

Z(l,  5)-  0 . 15283E-02-0. 11371E-02 

Z(l,  6)-  0 . 29420E-02  0.13780E-02 

Z(l,  7)— 0.14488E-01-0.54125E-02 

Z(l,  8)— 0.335S2E-01  0.19136E+00 

(TEL ,  NCOLX ,  NRO  VY-  8  8  1 

CURRENTS 


1 

1.000000 

0.683E-04 

13 

2 

1 . 000000 

0.683E-04 

13 

3 

1.000000 

0 . 683E-04 

13 

4 

1 . 000000 

0.683E-04 

13 

S 

1.000000 

0.683E-04 

13 

6 

1.000000 

0 . 683E-04 

13 

7 

1 . 000000 

0.683E-04 

13 

8 

1.000000 

0 . 683E-04 

13 

IC.RVTA-  1  0.35355E+00 

IC.RVTA-  2  0.35355E+00 

IC.RVTA-  3  0.35355E+00 

IC.RVTA-  4  0.35355E+00 

IC.RVTA-  5  0.35355E+00 

IC.RVTA-  6  0.35355E-00 

IC.RVTA-  7  0.35355E+00 

IC.RVTA-  8  0 . 35355E+00 

I,CVTA(I,1)»  1  0 . 34510E+00-0 . 76830E-01 

I , CVTA ( I , 1 ) «  2  0 . 34510E+00-0 . 76830E-0 1 

I,CVTA(I,1>-  3  0.34510E+00-0.76830E-01 

I , CVTA (I , 1)-  4  0 . 345 1 OE+OO-O . 7 6830E-0 1 

I.CVTAU.l)-  5  0 . 34510E+00-0 . 76830E-01 

I .CVTA (I , 1)«  6  0. 34510E+00-0. 76830E-01 

I.CVTAU.l)-  7  0 . 34510E+00-0 . 76830E-01 

I.CVTAU.l)-  8  0. 34510E*00-0 .76830E-01 

•••••ring  array  weights  before  nulling  (amp, phase) 


I- 

1 

CVTADB , CVTADG- 

0 . 00000 

-12.56098 

I- 

2 

CVTAr"  .  'UADG- 

0 . 00000 

-12.56099 

I- 

3 

CVTADB.CVTADG- 

0.00000 

-12.55098 

I- 

4 

CVTADB  .CVTADG- 

0.00000 

-12.65097 

I* 

6 

CVTADB.CVTADG- 

0 . 00000 

-12.55098 

I- 

6 

CVTADB ,  CVTADG- 

0.00000 

-12.55097 

1- 

7 

CVTADB ,  CVT  ADG- 

0 . 00000 

-12.55098 

I- 

8 

CVTADB .CVTADG - 

0 . 00000 

-12.55097 

Z(l, 

D-  0.91705E+02  0 

12403E*02 

ZU, 

2)— 0.33552E-01  0 

19136E+00 
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2(1,  3)— 0.14488E-01-0.64125E-02 
Z(l,  4)-  0.29420E-02  0.13780E-02 
Z(l,  5)-  0. 15283E-02-0 . 11371E-02 
Z(l,  6)-  0.29420E-02  0.13780E-02 
Z(l,  7)— 0.14488E-01-0.64125E-02 
Z(l,  8)— 0.33552E-01  0.19136E+00 


NEL.NCOLX.NROWY- 

8 

8  1 

II.PWKJ  IN  POWER* 

1 

0 . 10000E+05 

II.PWRJ  IN  POWER- 

2 

0 . 10000E+05 

II.PWRJ  IN  POWER* 

3 

0.31623E+02 

II.PWRJ  IN  POWER- 

4 

0 . 31623E-02 

II.PWRJ  IN  POWER- 

5 

0.12S89E- 

•09 

II.PWRJ  IN  POWER- 

6 

0. 12589E- 

■09 

II.PWRJ  IN  POWER- 

7 

0.12S89E- 

-09 

I, J. CNDB. PHASEN- 

1 

1 

36.35 

0.00 

I, J. CNDB. PHASEN- 

1 

2 

31 .68 

91.33 

I , J , CNDB , PHASEN* 

1 

3 

24.99 

-117.95 

I . J , CNDB , PHASEN* 

1 

4 

20.88 

6.26 

I, J. CNDB, PHASEN- 

1 

5 

20.72 

-0.14 

I, J. CNDB, PHASEN- 

1 

6 

20.88 

6.26 

I, J, CNDB, PHASEN- 

1 

7 

24.99 

-117.95 

I, J, CNDB, PHASEN- 

1 

8 

31.68 

91.33 

I, J, CNDB, PHASEN- 

2 

1 

31.68 

-91.33 

I, J, CNDB, PHASEN- 

2 

2 

27.05 

0.00 

I, J, CNDB, PHASEN- 

2 

3 

20.23 

155.55 

I, J, CNDB, PHASEN- 

2 

4 

8.12 

177.91 

I, J, CNDB, PHASEN- 

2 

S 

20.90 

-6.14 

I, J, CNDB, PHASEN- 

2 

6 

8.28 

177.98 

I, J, CNDB, PHASEN- 

2 

7 

20.24 

155.46 

I. J, CNDB, PHASEN- 

2 

8 

27.04 

0.00 

I, J, CNDB, PHASEN- 

3 

1 

24.99 

117.95 

I, J, CNDB, PHASEN- 

3 

2 

20.23 

-155.55 

I, J, CNDB. PHASEN- 

3 

3 

16.67 

0.00 

I, J. CNDB, PHASEN- 

3 

4 

20.26 

-155.49 

1. J, CNDB, PHASEN- 

3 

5 

25.01 

117.95 

I, J, CNDB, PHASEN- 

3 

6 

20.27 

-156.41 

I, J, CNDB, PHASEN- 

3 

7 

16.52 

o.oc 

I, J, CNDB. PHASEN- 

3 

8 

20.24 

-155.46 

I, J, CNDB, PHASEN- 

4 

1 

20.88 

-6.26 

I, J, CNDB. PHASEN- 

4 

2 

8.12 

-177.91 

I, J, CNDB, PHASEN- 

4 

3 

20.26 

155.49 

I, J, CNDB, PHASEN- 

4 

4 

27.08 

0.00 

I, J, CNDB, PHASEN- 

4 

5 

31.71 

-91.33 

I, J. CNDB, PHASEN- 

4 

6 

27.07 

0.00 

I, J. CNDB. PHASEN- 

4 

7 

20.27 

155.41 

I, J, CNDB, PHASEN- 

4 

e 

8.28 

-177.98 

I, J, CNDB. PHASEN- 

5 

1 

20.72 

0.14 

I, J. CNDB. PHASER- 

5 

2 

20.90 

6.14 

I. J. CNDB, PHASEN- 

5 

3 

25.01 

-117.95 

I, J, CNDB, PHASEN- 

5 

4 

31.71 

91.33 

I, J, CNDB, PHASEN- 

5 

5 

36.38 

0.00 

I, J, CNDB. PHASEN- 

5 

6 

31.71 

91.33 

I, J, CNDB. PHASEN- 

5 

7 

25.01 

-117.95 

I, J, CNDB, PHASEN- 

S 

8 

20.90 

6.14 

I, J, CNDB, PHASEN- 

6 

1 

20.88 

-6.26 

I, J, CNDB, PHASEN- 

6 

2 

8.28 

-177.98 

I. J. CNDB, PHASEN- 

6 

3 

20.27 

155.41 

I.J.CWDP  PM4CFW- 

* 

4 

77  A7 

0.00 

I.J.CKDB.PHASEN- 

6 

S 

31.71 

-91.33 

I , J , CHUB , PHASEN* 

6 

6 

27.08 

0.00 

I . J , CHDB , PHASEN* 

6 

7 

20.26 

165.49 

I, J.CNDB, PHASER* 

6 

8 

8.12 

-177.91 

I, J, CHDB, PHASEN* 

7 

1 

24.99 

117.96 

I, J, CHDB, PHASEN* 

7 

2 

20.24 

-166.46 

I, J. CHDB, PHASEN* 

7 

3 

16.  S2 

0.00 

I, J, CHDB, PHASEN* 

7 

4 

20.27 

-165.41 

I , J , CHDB . PHASEN* 

7 

5 

26.01 

117.95 

I, J.CNDB, PHASEN* 

7 

6 

20.26 

-165.49 

I, J. CHDB, PHASEN* 

7 

7 

16.67 

0.00 

I , J , CNDB , PHASEN* 

7 

8 

20.23 

-155.65 

I,J, CHDB, PHASEN* 

8 

1 

31.68 

-91.33 

I , J , CNDB , PHASEN* 

8 

2 

27.04 

0.00 

I , J , CNDB , PHASEN* 

8 

3 

20.24 

155.46 

I, J, CHDB. PHASEN* 

8 

4 

8.28 

177.98 

I, J.CNDB, PHASEN* 

8 

5 

20.90 

-6.14 

I, J.CNDB, PHASEN* 

8 

6 

8.12 

177.91 

I, J.CNDB, PHASEN* 

8 

7 

20.23 

165.55 

I, J.CNDB, PHASEN* 

8 

8 

27.05 

0.00 

••••••••*«igenvaluas ****** ••••****•**•*••* 


I.EVLDBN* 

1 

37 . 437 

I.EVLDBN* 

2 

37.195 

I.EVLDBN* 

3 

2.771 

I , EVLDBN* 

4 

3.103 

I.EVLDBN* 

5 

0.000 

I.EVLDBN* 

6 

0.000 

I.EVLDBN- 

7 

0.000 

I, EVLDBN- 

8 

0.000 

CPR0D1 .CPR0D2*  0. 

13839E+04-0 

85265E-13  0.10000E+01  O.OOOOOE+OO 

INR  NORMALIZATION 

PARAMETER. 

INRNOR*  1 

I-  1  WANDB* 

-12.318 

I-  2  WANDB* 

-8.734 

I*  3  WANDB* 

-7.602 

I-  4  WANDB* 

-8.734 

I*  5  WANDB* 

-12.318 

I-  6  WANDB* 

-8.734 

I*  7  WANDB* 

-7.602 

I*  8  WANDB* 

-8.734 

I.WANU.l)- 

1 

0 . 1 1307E+00-0 . 21414E+00 

I .WAN (I , 1)* 

2 

0 . 34937E+00  0.10851E+00 

I .WAN  Cl  ,1)* 

3 

0 . 32832E+00-0 . 2567 1E+00 

I , WAN (1,1)* 

4 

0.34937E+00  0.10851E+00 

I .WAN (1,1)- 

6 

0 . 1 1307E+00-0 . 21414E+00 

I.WANCI.D- 

6 

0.34937E+00  0.10851E+00 

I , WAN  C 1 ,  1  )* 

7 

0 . 32832E+00-0 . 25671E+00 

I , WAN (1,1)* 

8 

ve 

0 . 34937E+00  0 . 1085 1E+00 
array  weights  (amp. .phase)1 

I* 

1 

WANDB. WANDG- 

-4.71576 

-62.16515 

I- 

2 

WANDB.WANDG* 

-1.13211 

17.26469 

I- 

3 

WANDB, WANDG- 

0.00000 

-38.02192 

I- 

4 

WANDB.WANDG* 

-1.13210 

17.25474 

I- 

5 

WANDB, WANDG- 

-4.71576 

-62.16529 

I- 

6 

WANDB.WANDG* 

-1.13210 

17.25474 

I* 

7 

WANDB, WANDG* 

0 . 00000 

-38.02193 

I- 

8 

WANDB, WANDG- 

-1.13210 

17.25469 

••••••••cancellation*** •*••*• •*•••**•** 

CPR0D1.CPR0D2-  0.12271E+01  0.1218SE-13  0 . 10000E+01  0 . OOOOOE+OO 
I HR  NORMALIZATION  PARAMETER,  INRNOR*  1 
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INR-  QUI.ADAP  .CANCEL-  31.411  0.889  -30.622  DB 

CANCDB ,NRAN , AVECAN-  -30.52207  1  -30.62207 

—Input  data  file  lor  adaptive  nulling  with  two  auxiliary  probes, 
•••sdipjamhyper  data  file,  filename  edipjaahyper.datacirconicrr3 
ADIPOLE  NCOLX-8 .NROWY-l .HZIN-O.OOl .HLIN-2 .6 , 

DIIN-6 . 888 .  DYIN-6 . 888 ,  ARADIN-0 . 0039  , 

ICIRC-1  .CRADIN- 1 1 . 81  .HLSIN-0 . 25 , 

ItfL-0 , FCHZ- 120 . 0E6 , BVFHZ-1 . 0E0 , NFREQ-5 , IWR-0 , 

NPDI-0 ,  NXDUH-0 ,  NYDUN-0 ,  ER3-73 . 5 ,  SIG3-0 . E.TD3— 1 . 0 , 

ZLOAD-50 . 0 , Z CHAR-0 . 0 , NGEN-8 , IGEN-0 . 

BSCANS-1, 

THSINC-S.O, 

IMUT-1 .IBLTSL-0 , 

IENORM-1 , ICHEB-0 , SLLDB-20 . , EDGTDB-0 . . 

IPATRN-1 . IPRCOH-l . IANGLP-0 , 

NPHCT-0 ,  NTHPT-499 , 

THDR-180 .  .THDMIN— 90.  , 

RLSYIN-0 . 0 .RLSZIN-0 . 0 . 

NC0LXN-117 .NROWYN-l .NCOLZN-l , 

INEAR-1 . 

IPOL-2 , IGRNDP-0 , 

ITLTFR-0 , ITLTDP-0 , 

NFCOLI-1 .NFROWY-l , 

IANTX-1 , IANTY-0 . NPOWER-O , 

IPCONN-l , IPCONF-O , IPCUTF-0 . IPCFX-0 , IPCFY-0 , IPCFZ-0 , 

ITEK-0 , 

IQUAN-0 , IRNERR-0 . ELERDB-0 . 02 . ELERDG-0 . 2 . NBMQD- 1 2 , 

NRAN- 1 . NBADVT-32 , AWERDB-0 . 0 . AVERDG-0 . 0 . 

NAUX«7,IAUXA(1)«1 .2,3,4 ,5,6,7. 

IATTEN-1 , 

INRNOR-1, 

N JAKS-7 , ISLC-0 , AUXADB ( 1 ) -8-0 . 0 , 

PtfRJDB(l)-40. .40. , 5--99 . 0, 

YJAMIN(l)-7»0.0, 

XFOCIN-O . 0 , ZFOCIN-O . 0 , 

XJAHIN(l)-5.9,-5.9,2-0.0,2— 4.0, 1-0.0, 

ZJAMIN (1) -0.0, 0.0, 4. 0,-4. 0,3. 0,-3. 0,1 -0.0 
YNIN-O.O, 

NC0LXN-101 .NTHPT-101 , 

ININ— 16.  ,ZNIN— 15.  .RLSXIN-30.  .RLSZIN-30.  .NC0LXN-41  .NCOLZN-41 , 
AEND 

•••Output  data  file  for  adadptive  nulling  with  two  auxiliary  probes, 
•••sdipjamhyper  data  file,  filename  sdipjamhyper .datacirconicrr3 
IJAMINC1)  ,YJAMN(1)  ,ZJAMIN(1)-  5.900  0.000  0.000 


FGHZ-  0.1200000 

DX.DY.HL.ARAD-  0.17496 


IV, VTA-  1 
IV. VTA-  2 
IV. VTA-  3 
IV, VTA-  4 
IV, VTA-  5 
IV, VTA-  6 
IV, VTA-  7 
IV, VTA-  8 


0 . 5916E-02 
0 . 5916E-02 
0.5916E-02 
0 . 5916E-02 
0 . 6916E-02 
0.6916E-02 
0.5916E-02 
0 . 5916E-02 


0.2211E-02 
0.2211E-02 
0 . 221  IE-02 
0.2211E-02 
0 . 221  IE-02 
0.22' tE-02 
0. 2211E-02 
0.2211E-02 


0.17496 


0.06604 


array  mutual  impedance  matrix  (first  row)«»**«*» 


Z(l,  D-  0.91705E+02  0.12403E-02 
Z(l,  2)— 0.33552E-01  0.19136E+00 
Z(l,  3)--0 . 14488E-01-0 . 54125E-02 


0.00010 
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2(1. 

4)-  0.29420E-02 

0 . 13780E-02 

2(1, 

5)-  0.15283E-02- 

-0. 11371E-02 

2(1. 

6)»  0.29420E-02 

0.13780E-02 

2(1, 

7)— 0.14488E-01- 

-0.5412SE-02 

2(1. 

8)— 0.33582E-01 

0.19136E+00 

NEL.NCOLX.NROWY-  8 

CURRENTS 

8  1 

1 

1.000000 

0.683E-04 

13. 

2 

1.000000 

0.683E-04 

13. 

3 

1.000000 

0.683E-04 

13. 

4 

1.000000 

0.683E-04 

13. 

5 

1.000000 

0.683E-04 

13. 

6 

1.000000 

0 . 683E-04 

13. 

7 

1.000000 

0.683E-04 

13. 

8 

1.000000 

0.683E-04 

13. 

IC.RWTA-  1 

0.35355E+00 

IC.RWTA-  2 

0 . 35355E+00 

IC.RWTA-  3 

0 . 35355E+00 

IC.RWTA-  4 

0.3S355E+00 

IC.RWTA-  5 

0.3S3S5E+00 

IC.RWTA-  6 

0 . 35355E+00 

IC.RWTA-  7 

0 . 353S5E+00 

IC.RWTA-  8 

0  - 35355E+00 

I.CWTA(I.l)- 

1  0 . 34510E+00-0 . 76830E-01 

I.CWTA(I.l)- 

2  0 . 34510E+00-0 . 76830E-01 

I.CWTAd.l}- 

3  0 . 345 10E+00-0 . 76830E-01 

I.CWTAd.l)- 

4  0 . 34510E+00-0 . 76830E-01 

I.CWTA(I.l)- 

6  0 . 34510E+00-0 . 76830E-01 

I.CWTA(I.l)- 

6  0 . 345 10E+00-0 . 76830E-01 

I ,CWTA(I , 1)« 

7  0 . 34510E+00-0 . 76830E-01 

I ,CWTA (I , 1)» 

8  0 . 34510E+00-0 . 76830E-01 

•ring  array  weights 

;  before  nulling 

(amp.phase)** 

I* 

1 

CWTADB .  CWTADG- 

0.00000 

-12  55098 

I- 

2 

CWTADB ,  CWTADG- 

0.00000 

-12.65099 

I- 

3 

CWTADB.  CWTADG- 

0 . 00000 

-12.55098 

I- 

4 

CWTADB  .CWTADG- 

0 . 00000 

-12.55097 

I- 

5 

CWTADB .CWTADG- 

0 . 00000 

-12.55098 

I- 

6 

CWTADB, CWTADG- 

0 . 00000 

-12.55097 

I- 

7 

CWTADB. CWTADG- 

0 . 00000 

-12.55098 

I- 

8 

CWTADB.  CWTADG- 

0 . 00000 

-12.55097 

2(1. 

1)-  0.91705E+02 

0 . 12403E-02 

2(1. 

2)— 0.33552E-01 

0 . 19136E+00 

2(1, 

3)— 0. 14488E-01- 

■0 . 54125E-02 

2(1, 

4)-  0 . 29420E-02 

0. 13780E-02 

2(1 . 

5)-  0 . 15283E-02- 

■0.U371E-02 

2(1. 

6)-  0.29420E-02 

0.13780E-02 

2(1. 

7)— 0. 14488E-01- 

-0.5412SE-02 

2(1. 

8)— 0.33552E-01 

0. 19136E+00 

NEL. 

NCOLX.NROWY-  8 

8  1 

II.PWRJ 

IN  POWER-  1 

0. 10000E-05 

II.PWRJ 

IN  POWER-  2 

0.10000E+0S 

II.PWRJ 

IN  POWER-  3 

0.12S89E-09 

II.PWRJ 

IN  POWER-  4 

0.12589E-09 

II.PWRJ 

IN  POWER-  5 

0.12589E-09 

II.PWRJ 

IN  POWER-  6 

0. 12S89E-09 

II.PWRJ 

1  IN  POWER-  7 

0.12589E-09 

•covariance  matrix' 

I.J. 

CND8.PHASEN-  1 

1  36.35 

0.00 

I.J. 

CNDB.PHASEN-  1 

2  31.68 

91.33 
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I .J.CNDB, PHASES-  1  3 

I, J.CNDB, PHASES -  1  4 

I , J .CNDB.PHASEN-  1  5 

I, J.CNDB, PHASEN-  1  6 

I, J.CNDB.PHASEN-  1  7 

I, J.CKDB, PHASED-  1  8 

I, J.CKDB, PHASED-  2  1 

I, J.CKDB, PHASEN-  2  2 

1, J.CKDB. PHASED-  2  3 

I, J.CKDB, PHASER  2  4 

I, J.CKDB, PHASED-  2  S 

I, J.CKDB, PHASER-  2  6 

X. J.CKDB. PHASEN-  2  7 

I, J.CKDB. PHASEN-  2  8 

I, J.CNDB.PHASEN-  3  1 

I, J.CKDB. PHASEN-  3  2 

I, J.CNDB, PHASEN-  3  3 

I, J.CKDB, PHASEN-  3  4 

I, J.CKDB. PHASEN-  3  5 

I, J.CNDB, PHASEN-  3  6 

I, J.CNDB, PHASEN-  3  7 

1, J.CNDB, PHASEN-  3  8 

I, J.CNDB. PHASEN-  4  1 

I, J.CNDB, PHASEN-  4  2 

I. J.CKDB. PHASEN-  4  3 

I,  J.CKDB,  PHASEN-  4  4 

I. J.CNDB, PHASEN-  4  S 

I,  J.CNDB.  PHASEN-  4  6 

I,  J.CNDB,  PHASEN-  4  7 

I, J.CNDB, PHASEN-  4  8 

I, J.CNDB, PHASEN-  S  1 

I, J.CNDB. PHASED-  5  2 

I,  J.CNDB,  PHASEN-  S  3 

I, J.CNDB. PHASEN-  5  4 

I, J.CNDB. PHASEN-  5  5 

I. J.CNDB. PHASEN-  5  6 

I, J.CNDB. PHASEN-  S  7 

I, J.CNDB. PHASEN-  5  8 

I,  J.CNDB. PHASEN-  6  1 

1, J.CNDB. PHASEN-  6  2 

I, J.CKDB. PHASED-  6  3 

I, J.CKDB. PHASEN-  6  4 

I.  J.CNDB.PHASEN-  6  5 

I. J.CNDB. PHASEN-  6  6 

I. J.CNDB, PHASEN-  6  7 

I. J.CNDB, PHASEN-  6  8 

I, J.CNDB. PHASEN-  7  1 

I.  J.CNDB.  PHASEN-  7  2 

I, J.CNDB. PHASED-  7  3 

I, J.CKDB. PHASEN-  7  4 

I. J.CKDB. PHASEN-  7  5 

I, J.CNDB. PHASEN-  7  6 

I, J.CNDB. PHASEN-  7  7 

I,  J.CKDB,  PHASEN-  7  .! 

I, J.CNDB. PHASEN-  8  1 

I, J.CKDB. PHaSEN-  8  2 

I, J.CNDB. PHASEN-  8  3 

I. J.CNDB. PHASEN-  8  4 

I .J.CNDB.PHASEN-  8  5 
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.90 

6 

.17 

25 

.01 

-117 
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I , J ,CNDB .PHASEN* 

8 

6 

8.24 

177.97 

I , J , CNDB , PHASEN* 

8 

7 

20.24 

155.47 

I. J.CHDB, PHASEN* 

8 

8 

27.05 

0.00 

8*« 

I , EVLDBN* 

1 

37 . 195 

I.EVLDBN* 

2 

37.437 

I. EVLDBN- 

3 

0.000 

I.EVLDBN- 

4 

0.000 

I.EVLDBN- 

5 

0.000 

I.EVLDBN- 

6 

0.000 

I.EVLDBN- 

7 

0.000 

I.EVLDBN- 

8 

0.000 

CPR0D1 , CPR0D2-  0 . 13835E+04-0 . 66843E- 13  O.IOOOOE+Ol  O.OOOOOE+OO 
I HR  NORMALIZATION  PARAMETER,  INRNOR-  1 


I-  1  WANDB- 

-14.003 

I-  2  WANDB- 

-8.939 

I-  3  WANDB- 

-6.886 

I-  4  WANDB- 

-8.939 

I-  5  WANDB- 

-14.003 

1-  6  WANDB- 

-8.939 

I-  7  WANDB- 

-6.886 

1-  8  HANDB- 

-8.939 

I.VAN(I.l)- 

1 

0 . 48427E-01 -0 . 19349E+00 

I.MAN(I.l)- 

2 

0.3S635E+00  0.26440E-01 

I .WAN (I , 1 ) — 

3 

0 .43481E*00-0 . 12564E-00 

I .VAN (1 , 1 )  — 

4 

0.35635E+00  0.26441E-01 

I .WAN (I , 1)- 

5 

0 . 48427E-01-0 . 19349E+00 

I.VANU.l)- 

6 

0.35635E+00  0.26441E-01 

I ,UAN(I , 1)- 

7 

0 . 43481E+00-0 . 12S64E*00 

I , UANCI , 1 )* 

8 

0.35635E+00  0.26441E-01 

)«**• 

•  •4 

■•adaptive  array 

»* 

I* 

1 

WANDB.VANDG- 

-7.11720  -75.94859 

I- 

2 

WANDB.NANDC- 

-2.05298  4.24348 

I- 

3 

WANEB.WANDG- 

0.00000  -16.11658 

I- 

4 

WANDB.WANDG- 

-2.05297  4.24352 

I- 

5 

VANDB.VANDG- 

-7.11719  -75.94872 

I- 

6 

WANDB.WANDG- 

-2.05298  4.24352 

I- 

7 

WANDB.WANDG- 

0.00000  -16.11658 

I- 

8 

WANDB.WANDG- 

-2.05297  4.24348 

CPR0D1, 

lCPR0D2*  0. 10001E+01-0.63768E-14  0.10000E+01 

0 . OOOOOE+OO 

I  NR 

NORMALIZATION  PARAMETER.  INRNOR-  1 

INR* 

■  QUI.ADAP  .CANCEL- 

31.410  0.000  -31 

.410  DB 

CANCDB i 

.  NRAN , AVEC AN- 

-31.40966  1  -31.40966 

••• •••••••••••••makef ile  for  moment  method  software******* 

••••file  adipjamhyperMake***** 

EXECUTABLE  *  sdipjasihyper .  out 
OBIS  *  sdipjamhyper . o  dzabgnloss . o  \ 

dpack2o  zabgenloss.o  pack2.o  fwgh.o  \ 
chebafo  wquan.o  taylor.o  \ 
eigenv . o  reordr . o  ydipsubloss . o  \ 
plothyper.o  contek.o  plabel.o  \ 
conturek.o  circsubloss . o 
*  use  tabs  for  continues 
IL  *  /usr/lib/f 68881 /libm . il 
FLAGS  *  -01  -v  -f 68881 
*( EXECUTABLE)  :  t(OBJS) 

dislink  -1F77  t  (FLAGS)  -o  I  (EXECUTABLE)  KOBJS)  *(IL) 
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sdipjamhyper . o :  /home/ajf /hyperthermia/sdipjamhyper . f 

til  KFLAGS)  -c  /home/ajf /hyperthermia/sdipjamhyper . f  KIL) 

plothypar . o :  /home/ajf /hyperthermia/plothyper . f 

ill  KFLAGS)  -c  /home/ajf /hyperthermia/plothyper  .t  KIL) 

dzabgnloss .0:  /home/ajf /hyperthermia/dzabgnloss .f 

ill  KFLAGS)  -c  /home/ajf  /hyperthermia/dzabgnloss  .f  KIL) 

dpack2 . 0 :  /home /a j  f /mon j  am/dpack2 . f 

til  KFLAGS)  -c  /home/ajf /monjam/dpack2.f  KIL) 

contek . o :  /home/ajf /monjam/contek . f 

f 77  KFLAGS)  -c  /home/ajf /mon jam/contek . f  KIL) 

plabel . o :  /home/ a j  f /mon j  tr/plabel . f 

til  KFLAGS)  -c  /home/ajf/monjtr/plabel.f  KIL) 

conturek.o :  /home/ajf  /hyperthermia/conturek  .f 

ill  KFLAGS)  -c  /home/ajf/hvperthermia/conturek.f  KIL) 
zabgenloss . 0 :  /home/ajf /hyperthermia/zabgenloas . f 

ill  KFLAGS)  -c  /home/ajf /hyperthermia/zabgenloss  .f  $(IL) 
pack2.o:  /home/ajf /monjam/pack2 . f 

ill  KFLAGS)  -c  /home/ajf /mon;am/pack2.f  KIL) 
fvgh.o:  /home/ajf /monjam 'fwgh . f 

ill  KFLAGS)  -c  /home/ajf /monjam/fwgh . i  KIL) 

chebaf.o:  /home/ajf /mon jam/chebaf . f 

ill  KFLACS)  -c  /home/ajf /monjam/chebaf  .  f  KIL) 

wquan.o:  /home/ajf /monjam/wquan .  f 

ill  KFLAGS)  -c  /home/ajf /monjam/irvquan . f  SCIL) 

taylor . 0 :  /home/ajf /mon jam/taylor . f 

ill  KFLAGS)  -c  /home/ajf /monjam/taylor . f  KIL) 

eigenv . 0 :  /home/ajf /monjam/eigenv . f 

ill  KFLAGS)  -c  /home/ajf /monjam/eigenv.f  KIL) 
reordr.o:  /home/ajf /monjam/reordr . f 

ill  S(FLAGS)  -c  /home/ajf /monjam/reordr  .  f  KIL) 

ydipsubloss . o :  /home/aj f /hyperthermia/ydipsubloss . f 

ill  KFLAGS)  -c  /home/ajf /hyperthermia/ydipsubloss . f  KIL) 

circ subloss . o :  /home/ajf /hyperthermia/circsubloss . f 

ill  KFLAGS)  -c  /home/ajf /hyperthermia/circsubloss  . f  SCIL) 

©  M.I.T.  LINCOLN  LABORATORY  1991,  ALL  RIGHTS  RESERVED 

•».«.«file  sdipjamhyper . f.. 

C«**PR0GRAM  SDIPJAMHYPER  F  ---  ANALYZES  FINITE  ARRAYS  OF  DIPOLES 
C***IN  LOSSY  DIELECTRIC  OR  FREE  SPACE. 

C««*THE  DIPOLES  ARE  ASSUMED  TO  BE  ORIENTED  PARALLEL 

C»»*T0  THE  PLANE  OF  THE  GRID  FOR  A  PLANAR  ARRAY,  OR  THEY  CAN  BE 

C«**ARRANGED  IN  AN  ANNULAR  (RING)  ARRAY  CONFIGURATION. 

C.. .RECEIVING  CONDITIONS  ARE  ASSUMED. 

C... DOUBLE-PRECISION  VERSION 
PARAMETER  (NUNCHN=8) 

PARAMETER  (NUMAUT*8) 
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PARAMETER  (NUMJAM-7) 

PARAMETER  (NUMELM-8) 

PARAMETER  (NUMFRQ-5) 

PARAMETER  (NUMNPT-1681) 

COMPLEX  PSC1826)  ,CZ(900)  ,  VA(NUMELM)  ,  ZOfUMELM  .NUMELM) 

COMPLEX  VTA(NUMELM)  , VREF A, VRECVA( NUMELM) 

COMPLEX  VCV(NUMELM)  .VRECVX  (NUMELM) . C J , CSUMA ,  VIA 
CC**»THE  ABOVE  MATRICES  ARE  DIMENSIONED  BY  THE  NUMBER  OF  ELEMENTS 
COMPLEX  »16  COVNF  (NUMCHN .NUMCHN)  ,  COVNFI  (NUMCHN .NUMCHN) 
COMPLEX  *16  COVAAJ (NUMCHN, NUMCHN)  ,CINVCN(NUMCHN, NUMCHN) 
COMPLEX  *16  VCHA (NUMCHN, NUMFRQ) 

COMP  LEI  *16  VMAJMA(NUM JAM, NUMFRQ) 

COMPLEX  *16  VAUIJA(NUMAUI,NUMJAM, NUMFRQ) 

COMPLEX  *8  VAXCVA(NUMAUX  .NUMNPT) 

C***N0TE  DIMENS.  VCHA(NCHAN  .NFREQ)  ,  VMAJMA(NJAM,  NFREQ) 

C  VAUXJA(NAUX,NJAM, NFREQ) , VAXCWA(NAUI ,NPTS) 

COMPLEX  *8  VXMANA (NUMNPT) .ETASP.GAMSP 

COMPLEX  *16  EIGV  AN  (NUMCHN)  .EIGVEN  (NUMCHN  ,  NUMCHN) 

COMPLEX  *16  NANA (NUMCHN , 1) 

COMPLEX  *16  tfTCTR(l, NUMCHN) , WQSLC (NUMCHN , 1 ) 

COMPLEX  *16  WQNA(NUMELM.l)  ,CMPROD(l, NUMCHN )  ,UAN (NUMCHN,  1) 
COMPLEX  *16  WANDMA(NUMELM.l) 

COMPLEX  *16  CWTA(NUMELM.l)  .ETADP.GAMDP 
DIMENSION  ACALPH(NUMELM)  .RWTA(NUMELM)  .RUTB(NUMELM) 

DIMENSION  VI ADB ( NUMNPT ),VXAPH( NUMNPT) 

DIMENSION  XC(NUMNPT)  .YC(NUMNPT)  .ZC(NUMNPT) 

DIMENSION  THSD(l) , PHSD (1 ) ,THCT(1) ,PHCT(1) 

DIMENSION  X J AMIN (NUM JAM) , YJAMIN(NUMJAM) .ZJAMIN(NUMJAM) 
DIMENSION  XJAM(NUMJAM) , YJAM(NUMJAM) .ZJAM(NUMJAM) 

DIMENSION  PVRJDB(NUMJAM) 

DIMENSION  PCHADB( NUMNPT) 

REAL  *8  CNDB ,PHASEN,DPDCR, WKE(144) .EVLDBN 
REAL  *8  PREALN .PIMAGN 
REAL  *8  QINRDB , AINRDB , CANCDB .EXCPCB 
REAL  *8  QINRAA , AINRAA 

REAL  *8  CNCLNA .ELSGDB , ELSGDG , AUSGDB , AVSGDG 

REAL  *8  SUMC , SUMAA , AVECAN , AVCAA , SAVE , SA VESR 

REAL  *4  RX(IOO) ,RY(100) ,RZ(100) ,FM(41,41) ,DMF (41 ,41 ) 

REAL  *4  AUIADB(NUMAUX) , VATENA(NUMAUX) 

INTEGER  IAUXA(NUMAUX)  , LTMP ( NUMCHN )  .MTMP(NUMCHN) 

CHARACTER  DATNAM*3S ,  0UTNAM*35 

COMMON/A/DX ,DY .NCOLX .NROVY ,NEL,HZ ,HL , ARAD , ZLO AD , ZCHAR 

COMMON  /B/  NGEN.IGEN, THETAS, PHIS, IMUT.IBLTSL, IPATRN 

COMMON  /D/  FGHZ , RLAMDA , IWL , IS , NSCANS 

COMMON  /DEG/  TKDR.THDMIN 

COMMON  /PLT/  INEAR, IPLOTM 

COMMON/XL YL/XL.YL 

COMMON/NEAR/I ANTX , I ANTY .NPOWER , IPCONN , IPCONF , IPCUTF 
C0MM0N/NEAR2/IPCFX , IPCFY , IPCFZ 

COMMON  /P/IN , YN , ZN . RLSX , RLSY , RLSZ , NCOLIN , NROVYN , NCOLZN 

COMMON  /P2/EDGET , ICOMB .PUNFLX .PUNFLY 

COMMON  /CHEBY/  ICHEB ,SLLDB 

COMMON  /GROUND/  IGRNDP 

COMMON  /NORM/  IENORM.BIGNDB 

COMMON  /NORMAL/  INRNOR 

COMMON /PCENTR/PCDX IN 

COMMON  /WRITE/  IVR 

COMMON  /SPLOSS/  ETASP.GAMSP 

COMMON  /DPLOSS/  ETADP.GAMDP 


DIP00200 

DIP00210 


DIP00320 


DIP00340 

DIP00370 

D1P00380 

DIP00390 

DIP00400 

DIP00410 


DIP00460 

DIP00470 

DIP00480 

DIP00490 

DIP00500 

DIPOOSIO 

DIP00S20 

DIP00530 

DIP0054C 

DIP005S0 

DIP00560 

D„P00B70 

DIPC0580 

DIP00590 

DIP00600 


72 


COMMON  /F/  FHZ ,ER3 ,SIG3 ,TD3 
COMMON  /CIRCLE/  ICIRC .RADIUS ,HLS 

NAMELIST/DIPOLE/NCOLX , NROVY ,HZIN ,HLIN , ARADIN , ZLOAD ,NGEN , I GEN , 
1DXIN , DYIN , IMUT , NSCANS , THSD , PHSD , THSINC , IBLTSL , IPATRN , ER3 , SIG3 ,TD3 
2NPHCT , PHCT , NTHPT , ZCHAR , I WL , FCH2 , BVFHZ , NFREQ ,THDR, THDMIN .ICIRC, 
3XNIN . YNIN , ZNIN .NCOLXN , NROtfYN , NCOLZN , INEAR . IPRCOM , IANGLP . CRADIN . 

41 ANTX ,  IANTY , NPOVER . NFCOU , NFROVY ,  IPCONN ,  IPCONF .  IPCUTF ,  IPOL.HLSIN , 
BIPCFX , IPCFY , IPCFZ , ITEK , IGRNDP , IGAIN , ICHEB ,  IWR , IENORM , BIGNDB . 
6NXDUM . NYD'JH , NPDX  ,  EDGTDB ,  RLSXIN , RLS YIN . RLSZIN .  ISLC .  INUNIF . 

7XF0CIN , ZFOCIN . XJAKIN , YJAMIN . ZJAMIN , IATTEN , AUIADB , INRNOR , 

8N JAMS , PWRCDB , PtfRJDB . NAUI , IAUXA , IAUXB . ITLTPR , ITLTDP , 

9IQUAN .  IRNERR .  ELERDB .  ELEROG  ,  AVEROB  ,  AVERDG ,  NBMOD ,  NBADVT ,  NRAN ,  SLLDB 
C*«NOTE;  IF  ItfL-0  (INCHES).  IWL-2  (METERS) 

VRITE(6.2959) 

2959  FORMAT (IX, ’ENTER  INPUT  DATA  FILE  NAME  (typ.  idipjamhyper . data) ’ ) 
READ(5,«)  DATNAM 

OPEN (4 , FILE*DATNAM , FORM* ’ FORMATTED ’ ) 

HRITE(6,3959) 

3959  FORMAT (IX, 'ENTER  OUTPUT  DATA  FILE  NAME  (typ  sdipjamhypar .output) ’ 
R£AD(5 , »)  OUTNAM 

OPEN ( 8 , FILE*OUTNAM , FORM* ’ FORMATTED  > ) 

CALL  GETCP2(CPU1) 

PI*3. 141592654 
DCR-PI/180. 

DPDCR-DCR 
CJ*(0. ,1.) 

CINMTR*0 . 0254 

IGRNDP'l 

IGAIN-0 

ZLOAD-O.O 

ICHEB*0 

IWR*0 

IEN0RM=1 

EXCPCB*-1 .  ODC 

IQUAN*0 

NRAN*1 

TILTPR*0.0 

TILTDP=0 . C 

IATTEN*0 

IPOL* 1 

ER3*1 .0 

SIG3*0.0 

TD3*-1 .0 

ICIRC*0 

CRADIN*0.0 

KLSIN=0 . 0 

ISLC* 1 

INUNIF*0 

INRNOR-O 

READ(4, DIPOLE) 

IFtNFREQ . EC  1)  BVFHZ*0.0 
IF(ISLC.EQ.O)  WRITE(6, 81100) 

81100  FORMAT (IX,  ’FULLY  ADAPTIVE  ARRAY’) 

IFdSLC.EO.O.  AND.  ICHEB.  EQ.O)  INUNIF-1 
IF ( ICIRC . EQ . 1 )  WRITE(6,7000) 

7000  FORMAKlX,  ’RING  ARRAY  GEOMETRY’) 

fhz-fch: 


DIP00610 


DIP00660 


DIP00700 


DIP00710 

DIP00720 

DIP00730 

DIP00740 

DIP00750 

DIP00760 

DIP00770 

DIP00780 

DIP00790 

DIP00800 

DIP00810 

DIP00820 

DIP00830 

DIP00840 

DIP00850 

DIP00860 

DIP00870 


DIP00880 


IXY-0 

IF(NROWYN.EQ.l)  IXZ-1 
IF(NCOLZN.EQ.l)  IXY-1 
IFdPOL.  EQ  2)  IPECGM-0 
IFdPRCOM.EQ.O)  VRITE(6 ,7898) 

7898  FORMATdX, ’NO  PROBE  COMPENSATION’) 

WRITE(6,9276)IGRNDP 

9276  FORMATdX’ GROUND  PLANE  PARAMETER,  IGRNDP- ’ , 14) 

CC  DO  1616  IX-l.NAUX 

CC  IF(IATTEN.EQ.O)  AOXADB(IX)-0.0 

CC  VRITE(6, 2318)11, AUXADBdX) 

C2318  FORMATdX,  ’IX, AUXADB(IX)-’  ,14, 2X.F12. 2) 

CC  VATENA(IX)-10.-»(AUXADBdX)/20.) 

C1615  CONTINUE 

IF(NJAMS.GT.l)  WRITE(8,2009)XJAMIN(1),YJAMIN(1) .ZJAMIN(l) 
2009  FORMATdX,  ’XJAMIN(l)  ,YJAMIN(1)  ,ZJAMIN(1)-’ .3F12.3) 
WRITE(6,6987)  BWFHZ  ,NFR£Q 
6987  FORMATdX,  ’BWFHZ,  NFREQ=  '  .£12.5, 2X  .  15) 

IF(ITEK . EQ - 0)CALL  COMPRS 

IF (ITEK . EO . 1 ) C ALL  TEKALL (4014,480,0,1,0) 

CC  IF(INEAR  EQ.O)  CALL  PRNTDA 
IF(IMUT.EQ.O)  ZL0AD-1.0 
IFdTLTPR.EQ.  1)  TILTPR-45 . 

IFdTLTDP.EQ.  1)  TILTDP-45 . 

CTPR«COS(TILTPR*DCR) 

STPR“SIN(TILTPR*DCR) 

CTDP-COS (TILTDP-DCR) 

STDP-SIN (TILTDP*DCR) 

FGHZ-FCHZ/1.0E9 

NEL-NCOLX-NROUY 

RNEL-NEL 

XLIN-DXIN*(NC0LX-1) 

YLIN-DYIN* ! NROVi  -i) 

NAC0LX-NC0LT-;<P0X-2*NIDUM 

NAROUY-NROUY-2-NYDUM 

nael*nacolx*narovy 

PCDXIN*NPDX*DXIN 
NAUXP1-NAUX-*1 
NAUXP2-NAUX+2 
IF(ISLC.EQ.l)  NMAX-NAUXP1 
IF(ISLC.EQ.O)  NMAX-NEL 
NMAXP1-NMAX+1 
ELSGDB*ELERDB*SQRT(3 . ) 

ELSGDG*ELERDG*SQRT  ( 3  ) 

A9SGDB-AVERDB*SQRT(3. > 

AVSGDG-AtfERDG*SQRT ( 3 . ) 

:nitrn=i 

IF(ISLC.EQ.O)  THEN 
NAUX-NEL 
NAUXP1-NAUX-M 
NAUXP2-NAUX-2 
DO  79130  1*1, NEL 
IAUXA(I)*I 

VRITE(6, 76767)1 ,IAUXA(I) 

76767  FORMATdX,  ’I, IAUXA(l)*’  .215) 

79130  CONTINUE 
ENDIF 

DO  1615  IX-l.NAUX 
IFdATTEN.EO  O,  AUXADB(IX)-0.0 


DiP00890 

DIP00900 

DIP00910 

DIP00920 

DIP00930 

DIP00940 

DIP00950 

DIP00960 

DIP00970 

DIP00980 

DIP00990 

DIP01000 

DIP01010 

DIP01020 

DIP01030 

DIP01040 

DIP01060 

DIP01060 

DIP01070 

DIP01080 

DIP01090 

DIPOllOO 

DIP01110 

DIP01120 

DIP01130 

DIP01140 

DIP01160 

DIP01160 

DIP01170 

DIP0H80 

DIP01190 

DIP01200 

DIP01220 

DIP01230 


DIP012S0 
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WRITE ( 6 , 23 1 8 ) I X , AUXiDB (I X ) 

2318  F0RMAT(1X, • IX , AUXADB(IX)*  \  14 .21 ,F12 .2) 

VATENA(IX)*lO . — (AUXADBCIX) /20 . ) 

1615  CONTINUE 
30  CONTINUE 
NR-1 

IF(NGEN.EQ.l)  NSCANS-0 

ICC-NEL 

ICC1-ICC 

IF(IBLTSL.EQ.l)  ICC1-NR0WY 
CC  IFCIWL.EQ.  1)G0  TO  50 

WRITE(8,40)FGHZ 

40  FORMATS, IX. ’FGHZ=' ,F16.7) 

. . — • . . . . . . . — — — 

c».». -COMPUTE  FREE  SPACE  LAMBDA— . 

C  —  —  .ALL  UNITS  HAVE  TO  BE  IN  METERS 
RLAMDA-2 . 997925E10/FCHZ/2 . 54 
C— -PARAMETER  CONVERSIONS  TO  PROPER  UNITS 
55  CONTINUE 

IF(IWL . EQ . 2)  GO  TO  50 
DX-DXIN-CINMTR 
DY-DYIN-CINKTR 
HL-HLIN-CINMTP. 

HLS-HLSIN-CINMTR 
ARAD-ARADIN-CINMTR 
HZ-HZIN-CINMTR 
RLSX-RLSXIN-CINMTR 
RLSY-RLSYIN-CINMTR 
RLSZ-RLSZIN-CINMTR 
XN-XNIN-CINMTR 
YN-YN'IN-CINMTR 
ZN=ZNIN»CINKTR 
XFOC-XFOCIN-CINMTR 
ZFOC-ZFOCIN-CINMTR 
XL-XLIN-CINMTR 
YL-YLIN-CINMTR 
PCDX-PCDXIN-CINMTR 
RADIUS-CRADIN-CINMTR 
50  CONTINUE 

WRITE C 8 . 60 ) DX , 0 Y , HL , ARAD 
60  FORMAT ( IX , ’ DX . DY , HL , ARAD= ’ ,2X ,4F14 ■ 5) 

C 

C— •  COMPUTE  CALIBRATION  CONSTANTS  (PHASE  ONLY)  TO 
C— MAXIMIZE  GAIN  (FOCUS  ANTENNA)  TO  NEAR  FIELD  RANGE 
C 

C  —  -PHASE  CENTER  'A'  VOLTAGE  EXCITATION 

X0A--XL/2 . -NXDUM-DX- (NACOLX- 1 ) /2 . *DX 
YP-0.0 

IF(ICIRC.EQ.O)  WRITE(6, 24690) 

24690  F0RMAT(1X, ’CALLING  NFDPX2 ’ ) 

IF(ICIRC.EQ.O) 

2CALL  NFDPX2 ( CTPR , STPR , CTDP , STDP , XFOC , YP , ZFOC , XOA , 0 , VA , VREFA) 
IFdCIRC.EQ.  1) 

2CALL  NFDPC2 ( XFOC , YP , ZFOC , XOA , 0 , VA , VREFA ) 

C-  — SAVE  INCIDENT  VOLTAGES 
DO  65  IV-l.NEL 
VTAdV)-VA(IV) 

WRITE(6,3757) IV, VTA (IV) 

WRITE (8 ,3757)IV,VTA(IV) 
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3767 

FORMAT ( IX , > I V , VTA- ’,14 , 2X . 2E12 . 4) 

DIP01760 

66 

CONTINUE 

DIP01770 

WRITE(6,9876)VTA(2) 

DIP01780 

9876 

F0RMAT(1X, ’VTA(2)-’ ,2E12.4) 

DIP01790 

NB-NCOLX 

DIP01800 

IDMB-NROWY 

DIP01810 

IDM-NEL 

DIPC1820 

ICC-NEL 

DIP01830 

IBLT-IBLTSL 

DIP01840 

IDM1-IDM3 

DIP01850 

IF(IBLT.EQ.O)  IDM1-NEL 

DIP01860 

12-1 

DIP01890 

8889 

CONTINUE 

DIP02350 

CALL  ZMATRI (CTDP , STDP .  IBLTSL ,  ICC1 ,  ICC . CZ , Z) 

C— SOLVE  SYSTEM  OF  EQUATIONS  FOR  THE  UNKNOWN  CURRENTS 

IF(IBLT.EQ.l)  GO  TO  240 

DIP02370 

ISYM-0 

DIP02380 

112-1 

DIP02390 

WRITE (6 ,9876) VTA (2) 

DIP02400 

WRITE(6 ,6110) 

DIP02410 

6110 

FORMAT (IX, ’CALL  CROUT’) 

DIP02420 

CAT  I  CROUT ( 2, V A, ICC 1 ,ICC,ISYM,IWR,I12,NEL) 

DIP02430 

112-2 

DIP02440 

WRITE(6,9876)VTA(2) 

DIP02450 

GO  TO  255 

DIP02460 

240 

IENTRY-4 

DIP02470 

CALL  BLTSOL (CZ . VA , PS , NCOLX , IDMB , I ENTRY ) 

DIP02480 

I  ENTRY -3 

DIP02490 

256 

WRITE(8,270) 

DIP02500 

270 

FORMAT ( IX , ’CURRENTS’) 

DIP02510 

IF (NEL . LT . 40)  CALL  CNORM(VA.NEL) 

DIP02520 

WRITE(6 , 280) 

DIP02530 

280 

FORMAT (IX. ’AFTER  CURRENTS  SOLUTION’) 

DIP02540 

C— VA  ARE  CURRENTS  (AMPERES)  NOW 

DIP02550 

C- 

DIP02560 

C... COMPUTE  RECEIVED  VOLTAGES 

DIP02570 

DO  285  IC-l.NEL 

DIP02580 

WRITE  (6 , 7531 )  ZLOAD 

DIP02590 

7631 

FORMAT (IX. ’ZLOAD- ’ ,F12 . 5) 

DIP02600 

WRITE (6 ,1134)IC,VTA(IC) 

DIP02610 

VRECVA(IC)-VA(IC) -ZLOAD 

DIP02620 

WRITE (6 ,1174) IC, VRECVA (IC) 

DIP02630 

1134 

FORMAT ( IX , ’ IC , VTA- ’ , 14 ,2X ,2E12 . 4) 

DIP02640 

1174 

FORMAT ( 11, ’IC, VRECVA-’ , 14 ,2X ,2E12 .4) 

DIP02650 

WRITE (6, 1135)IC,VA(IC) 

DIP02660 

1135 

FORMAT ( 11 , ’ IC, VA- ’ , 14 ,2X ,2E12 .4) 

DIP02670 

VRADB-20 . -AL0G10 (CABS (VRECVA ( IC) ) ) 

DIP02680 

VRAPH-ATAN2(AIMAG(VRECVA (IC) ) , REAL ( VRECVA ( IC) ) ) /DCR 

DIP02690 

C»* -COMPUTE  CALIBRATION  CONSTANTS  (PHASE  ONLY) 

DIP02700 

ACALPH(IC)— VRAPH 

DIP02710 

286 

CONTINUE 

DIP02720 

DO  7777  IC-l.NEL 

DIP02730 

WRITE  (  6 . 6667 )  I C ,  AC  ALPH  ( I C ) 

DIP02740 

CC 

WRITE (8 ,6667) IC , ACALPH( IC) 

DIP02750 

6667 

FORMAT ( IX , ’ IC , ACALPH- ’ , 14 , 2X ,F10 . 2 , ’  DECS’) 

DIP02760 

7777 

CONTINUE 

DIP02770 

C---CQMPUTE  NEAR  FIELD  PATTERN  OF  FOCUSED  ARRAY 

DIP02780 

c 

DIP02790 

TOMPUTE  BEAMFORMER  WEIGHTS  (I.E.  TAPER) 

DIP02800 

76 


IF(INUNIF.EQ.O) 

2CALL  VRCVWT ( NACOLX , NAROWY , NXDUM , NYDUM , EDGTDB , RWTA , RUTB ) 

DO  1199  IC-l.NEL 

IFCINUNIF . EQ.  1 .  AND .  ISLC .EQ . 0)  RWTA(IC)-1 . /SQRT(RNEL) 
URITE(6,S111)IC,RWTA(IC) 

WRITEC8 ,511 1 ) IC , RWTA (IC) 

5111  FORMAT C 11, ’ IC.RWTA- ’ ,14,21 ,E12 . 5) 

1199  CONTINUE 

DO  5333  KC-l.NEL 

CWTACKC , 1)*RWTA(KC)  *CEXP(CJ»ACAXJ>H(KC)»DCR) 

6333  CONTINUE 

IF(IQUAN.EQ.l) 

2CALL  ADQUAN (NEL.CVTA, NBMOD , IRNERR , INITRN , ELSGDB , ELSGDG ) 

BIGWDB— 299.0 
DO  91020  I-l.NEL 
WRITE(6, 61910) I ,CWTA(I , 1) 

WRITEC8 ,61910) I .CWTACI , 1) 

61910  FORMAT (11, ’ I  ,CWTA  (I ,  D*  ’ , I4.2X.2E12.5) 

CWTADB-20 . »DLOG 10 (CDABSC CWTACI , 1) ) ) 

IFCCWTADB.GT. BIGWDB)  BIGWDB-CVTADB 
91020  CONTINUE 

DO  91120  1*1, NEL 

CWTADB-20 . •DL0G10CCDABS (CWTA (1,1))) -BIGWDB 
CVT  ADG-DATAN2 (DIMAGCCWTACI.l)), DREAL  C  CWTA ( 1 , 1 ) ) ) /DCR 
WRITE ( 6 . 35990 ) I , CWTADB , CWT ADG 
WRITE ( 8 , 35990 ) I , CWT ADB . CWT ADG 
35990  FORMAT (IX , ’I« ’ , 14 , 2X , 'CWTADB , CWT ADG* ' , IX ,2F14 . 5) 

91120  CONTINUE 
INITRN-2 

C***PERFORM  NEAR  FIELD  SCAN  WITH  CW  RADIATING  DIPOLE 
IF( INEAR. EQ.O)  GO  TO  390 
WRITE(6,35C) 

WRITE ( 6 , 360)XNIN , YNIN , ZNIN , RLSXIN ,RLSYIH , RLSZIN , NCOLXN , NROWYN , 
2NC0LZN 

WRITE(6 ,320) IWL 
320  FORMATC IX , ’ IWL* ’ , 14) 

WRITE(6,330) 

330  FORMAT (IX, ’CHANGE  NEAR  FIELD  SCAN  PARAMETERS7,  ICHANG*1>) 

READ  ( 5 ,  •  ,  END-370 )  I  CHANG 
imCHANG.EQ.O)  GO  TO  370 
340  WRITE(6,350) 

350  FORMAT ( IX  ,  ’XN ,  TN ,  ZN ,  RLSXIN ,  RLSTIN ,  TI.'ZIN ,  NCOLXN , NROWYN , NCOLZN-  * ) 
READ (5 , * )XNIN ,  YNIN , ZNIN , RLSXIN, RLSY1N , RLSZ IN , NCOLXN .NROWYN ,NCOLZN 
WRITE (6 , 360) XNIN,YNIN,ZNIN .RLSXIN ,RLSYIN , RLSZIN , NCOLXN .  DWYN. 
2NC0LZN 

360  FORMAT (1X,6F10.3,2X,3I5) 

IF (IWL . EQ . 2)  GO  TO  370 
XN-XNIN-CINMTR 
YN-YNIN*CINMTR 
ZN*ZNIN*CINKTR 
RLSX-RLSXIN*CINMTR 
RLSY-RLSYIN*CINMTR 
RLSZ-RLSZIN*CINMTR 
370  CONTINUE 
C**ALL  DIMENSIONS  IN  METERS 
NELN-NCOLXN-NROWYN 
NPTSN*NELN»NCOLZN 

C*«*SET  DEFAULT  VALUES  FOR  DXN.DYN.DZN 
DXN-0.0 
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DYN-0.1 

DZN-0.0 

IF (NCOLXN . GT . 1 ) DXN-RLSX/ (NCOLXN- 1 ) 
IF(NROWYN.GT.l)DYN-RLSY/(NROWYN-l) 

IF (NCOLZN . GT . 1 ) DZN-RLSZ/ (NCOLZN - 1 ) 

IC-0 

BI6HDB— 299.0 

DO  3000  ICOLZN-1, NCOLZN 

ZPOS-ZN+DZN- ( ICOLZN- 1) 

DO  3000  IROWYN-l.NROWYN 
Y-YN+DYN* ( IROWYN- 1 ) 

DO  3000  ICOLXN-1 .NCOLXN 
CALL  GETCP2(CPU2) 

CPUSUB-CPU2-CPU1 
WRITE (6, 7319) CPUSUB 
7319  FORMAT ( II , ’ CPU  SUBTOTAL- ’ ,F14. 2) 

IC-IC+1 

X-IN+DXN • ( ICOLXN- 1 ) 

XC(IC)-X 

YC(IC)-Y 

ZC(IC)-ZPOS 

WRITE(6,6969)IC,XC(IC) .YC(IC) ,ZC(IC) 

6969  FORMATUX, ’ IC.XC , YC.ZC- ’,I4,2X ,3F12 . 3) 

IF(ICIRC.EQ.O) 

2CALL  NFDPX2(CTPR,STPR,CTDP,STDP,X>Y>ZP0S,0.0,O,VCV,VREFCW) 
IF(ICIRC . EQ . 1) 

2CALL  NFDPC2 ( X , Y , ZPOS , 0 . 0 , 0 , VCW , VREFCV ) 

IF(IMUT.EQ.O)  GO  TO  4255 

C— SOLVE  EACH  SYSTEM  OF  EQUATIONS  FOR  THE  UNKNOWN  CURRENTS 
IF(IBLT.Eq.l)  GO  TO  4040 
112-2 

WRITE(6 ,6110) 

CALL  CROUT (Z, VCV. ICC1 , ICC, ISYM, IWR, 112, VEL) 

GO  TO  4255 
4040  I ENTRY-3 

CALL  BLTSOL (CZ . VCW , PS , NCOLX , IDHB , I ENTRY) 

4255  CONTINUE 

c... COMPUTE  RECEIVED  VOLTAGES  FOR  PRESENT  SCAN 
DO  3285  IIIC-l.NEL 
VRECVX(IIIC)-VCW(IIIC)«ZLOAD 
3285  CONTINUE 

C—STORE  AUX.  CHANNEL  VOLTAGES 
DO  1681  IAX-l.NAUX 

C--MODIFICATION  TO  INCLUDE  AUX.  ATTEN. 

VAXCWA ( I AX , IC ) -VRECVX (I AUX A ( I AX ) ) »VATENA ( I AX ) 

CC  WRITEC6 ,3231) IAX  , IC ,XC(IC) , VAXCWACIAX , IC) 

3231  FORMATC IX, ’IAX .IC.XC, VAXCWACIAX, IC)-’ ,214 ,2X ,F12 . 3 ,21 ,2E12 . 5) 
1681  CONTINUE 
C—PERFORM  BEAM  FORMATION 
CSUMA-(0. ,0.) 

DO  6444  KC-l.NEL 

CSUMA-CSUMA+VRECVX (KCl -CWTA Uc , 1) 

6444  CONTINUE 
VXA-CSUMA 

IF(CABS(VXA) . EQ . 0 . ) VIA- ( 1 . E- 10 , 0 . ) 

VXADB(  10-20 .  -AL0G10  (CABS  (VIA)  ) 

IF(VXADB(IC)  .GT.BIC-UB)  BIGNDB-VXADBdC) 

VXAPH (IC ) - ATAN2 (AIMAG(VXA) , REAL (VIA) )/DCR 

VIMANA (10-10 .  ••  (VXADBCIC) /20 .  )-CEXP(CJ-VIAPH(IC)-DCR) 

CC  WRITE (6, 4457 )IC,VXMANA(IC) 
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CC  WRITE(6,4457)IC,VIMANA(IC) 

4457  FORMATdX,  ’IC.VXMANA-’  ,14, 2X.2E12. 6) 

WRITE(6,6429) IC ,XC(IC) , VIADB (IC) 

6429  FORMATdX,  ’ IC,X , VIADB- ’ ,I4,2X ,F10 . 2 ,2X ,F12 .2) 

3000  CONTINUE 

IF(IENORM.EQ.O)  GO  TO  2500 
WRITE(6 ,3765)BIGNDB 
3766  FORMATdX,  ’BIGNDB- ’ ,F12 . 2) 

DO  3020  IC-1, NPTSN 
CC  WRITE (6, 2343) IC, VIADB (IC) 

VXADB (IC) -VIADB (IC) -BIGNDB 
CC  WRITE (6 ,3343) IC .VXADB (IC) 

3020  CONTINUE 
2500  WRITE(6,3030) 

3030  FORMATdX, 'WANT  TO  PLOT  NEAR  FIELD  CUTS,  IPLOTN- 1>) 

READ ( 5 , * ) IPLOTN 

C---NEXT  LINE  ADDED  TO  AVOID  RUN  TIME  ERROR 
IF(IPCONN.EQ.l)  IPLOTN-O 
IF ( IPCONN .EQ.O)  THEN 
WRITE(30, 18889) 

WRITE ( 30 , 18888 ) XN , NCOLXN , DXN , ZN , NCOLZN , DZN 
WRITE  (30, 4547) 

4547  FORMATdX,  ’IZ, IX. VXADB(IC)  ’) 

ENDIF 

IPP-0 

DO  7788  IZ-1, NCOLZN 
DO  7788  IY-l.NROWYN 
DO  7788  IX- 1,  NCOLXN 
IPP-IPP+1 

CC  DO  7788  IPP- 1, NPTSN 

CC  WRITE (6 ,3343) IPP .VXADB (IPP) 

3343  FORMAT (IX , ' IPP , VXADB- 1 , 14 .2X.F12 . 2) 

IF ( IPCONN. EQ.O)  WRITE ( 30, »)  IZ, IX, VXADB (IPP) 

7788  CONTINUE 

3040  IF ( (NCOLXN . GT . 1 . OR . NROWYN . GT . 1 ) . AND . IPLOTN . EQ . 1 ) 

1 CALL  PLOTTR ( NCOLXN , NROWYN , NCOLZN , XC , YC , ZC , NPTSN , VXADB , VX APH . 
2VIADB.VXAPH) 

IF(NELN.EQ.l. AND. IPLOTN. EQ  1) 

1 CALL  PLOTAX (NPTSN , ZC , VXADB , VXAPH , VXADB , VX APH ) 

IF  (IPCONN.  EQ.O)  GO  TO  3939 
C---THIS  SECTION  FOR  CONTOUR  PLOTS 
NCLIN2-NC0LXN+2 
NRWYN2-NR0WYN+2 
NCLZN2-NC0LZN+2 
DO  3777  IX- 1, NCOLXN 
RX ( IX ) - ( XN+DXN- ( IX- 1 ) ) /CINMTR 
3777  CONTINUE 

DO  3008  IY-l.NROWYN 
RY(IY)-(YN+DYN«(IY'l)) /CINMTR 

3008  CONTINUE 

DO  3009  IZ-1, NCOLZN 
RZ(IZ)-(ZN+DZN-(IZ-1)) /CINMTR 

3009  CONTINUE 
IC-0 

IF(IXZ.EQ.l)  WRITE(30, 18889) 

18889  FORMATdX,  ’XN, NCOLXN, DXN, ZN, NCOLZN, DZN-’) 

IF ( IXZ . EQ . 1 )  WRITE! 30 , 18888 ) IN , NCOLXN , DXN , ZN .NCOLZN , DZN 
18888  FORMAT (IX, El^ .5,I5,2X,E14.5,2X,E14.5,I5,2X,E14.5) 
IF(IXZ.EQ.l)  WRITE (30, 4546) 
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4646 

F0RMATC1X, ’IZ,IX,FM(IZ,IX) ’) 

DO  3022  IZ-l.NCOLZN 

DO  3022  IY-l.NROVYN 

DIP04150 

DO  3022  IX-l.NCOLXN 

DIP04160 

IC-IC+1 

IF(IXY.EQ.l)  FM(IY , II)-VXADB(IC) 

IF(IIZ.EQ.l)  FM(IZ,IX)-VIADB(IC) 

IF ( IXZ . EQ . 1 . AND . IPCONN . EQ . 1 )  WRITE(30,»)  IZ, IX ,FM(IZ , IX) 

DIP04170 

3022 

CONTINUE 

IF(IXY.EQ.l) 

DIP04190 

2CALL  PLCONT (NCOLXN  .NROVYN ,NCLIN2 .NRWYN2 ,RX ,RY,FM ,DMF , -60 . , 

DIP04200 

310. . 6.1) 

IF(IIZ.EQ.l) 

2CALL  PLCONT (NCOLXN , NCOLZN , NCLIN2 , NCLZN2 .RX.RZ.FM.DMF.-BO., 

310. . 5.1) 

DIP04210 

3939 

CONTINUE 

DIP04220 

CC 

WRITE (26,*) NCOLXN 

DIP04470 

DO  1767  IDDD-1, NCOLXN 

DIP04480 

CC 

WRITE ( 25 . • ) VXADB (IDDD) 

1767 

CONTINUE 

DIP04600 

WRITE(6,3050) 

DIP04S10 

3050 

FORMAT (IX  ’PLOT  NEAR  FIELD  AGAIN?,  IPLA-1’) 

DIP04520 

READ(S , *)IPLA 

DIP04530 

IF(IPU.Eq.l)  GO  TO  3040 

DIP04540 

390 

CONTINUE 

DIP04550 

ISTOP-O 

DIP04560 

IF(ISTOP.EQ.l)  GO  TO  9999 

DIP04570 

C—CALL  PRNTDA  (PRINT  PARAMETERS) 

DIP04580 

C370 

CALL  PRNTDA 

DIP04590 

IF(IPATRN.EQ.O)  GO  TO  440 

DIP04630 

410 

WRITE(6,420) 

DIP04640 

420 

FORMATOX ,  ’  ISYLMBL  FOR  PLOTTING,  LT.  0  THEN  NOT  USED’) 

DIP04650 

READ(5,»)ISYMBL 

DIP04660 

WRITE (6 .430) 

DIP04700 

430 

FORMAT (IX, ’WANT  TO  PLOT  PATTERNS  AGAIN?,  IPFNA-1’) 

DIP04710 

READ(5,*)IPFNA 

DIP04720 

IF(IPFNA.EQ.l)  GO  TO  410 

DIP04730 

440 

CONTINUE 

DIP04740 

C---THIS  SECTION  FOR  COVARIANCE  MATRIX  COMPUTATION 

IF(NJAMS.EQ.O)  GO  TO  9999 

DIP04750 

DO  7999  ICH-l.NMAX 

DIP04770 

DO  7999  JCH-l.NMAX 

DIP04780 

COVNF (ICH, JCH)*DCKPLX(0 -ODO ,0 -ODO) 

DIP04790 

7999 

CONTINUE 

DIP04800 

8888 

IF(NJAMS.EQ.O)  GO  TO  4444 

DIP04810 

C--THIS  SECTION  FOR  JAMMER  COVARIANCE  MATRIX 

DIP04820 

FMINHZ-FCHZ-BWFHZ/2 . 

DELFHZ-0 . 0 

IFCNFREQ.GT. 1)  DELFHZ-BWFHZ/(NFREQ-1) 

DIP04830 

DO  600  IFR-1 ,  NFREQ 

DIP04850 

FHZ-FMINHZ+DELFHZ- ( IFR- 1 ) 

DIP04860 

FGHZ-FHZ/1 . 0E9 

DIP04870 

C---COMPUTE  FREE  SPACE  WAVELENGTH  AT  EACH  FREQUENCY 

RLAMDA-2  997925E10/FHZ/2 . 54 

C---NOTE :  GAMMA-  ALPHA  +J  BETA 

c*** 

AND  RLAMDA-2  PI/  BETA  (REF.  HAYT  PG.  334) 

C—.THUS  NEED  TO  COMPUTE  GAMMA,  AND  ETA  FOR  EACH  FREQ. 

IF(IWL.EQ.2)  CO  TO  8789 
DX«DXIN»CINMTR 
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DY-DYIN-CINMTR 
HL-HLIN-CINMTR 
ARAD-ARADIN-CINMTR 
HZ-HZIN-CINMTR 
8789  CONTINUE 

CALL  ZMATRX  ( CTDP ,  STOP ,  IBLTSL ,  ICC  1 ,  ICC ,  CZ ,  Z ) 

C---COMPUTE  ELEMENT  INDUCED  VOLTAGES  DUE  TO  JAMMER  SOURCES 
CALL  V JAMMRCN JAMS , NEL . PWRJDB , ICC1 , ICC , PS . CZ , Z , CVTA 
1 ,NB , IDMB , NAUX , IAUXA , ZLOAD , I J AMIN , Y J AMIN , Z J AMIN , 

2IFR . NFREQ , CTPR , STPR , CTDP , STDP , VMAJHA , VAUX J A) 
WRITE(6,1234)IFR,VMAJHA(1,IFR: . VAUXJACl , 1 , IFR) 

WRITE (8 , 1234) IFR, VMAJMA(1 , IFR) , VAUX JA(1 , 1 , IFR) 

1234  FORMAT (IX, ’  IFR .  VMA  JMA ,  VAUXA- ’ .I4.2X.4E11 .4) 

600  CONTINUE 

C---FORM  RECEIVED  VOLTAGE  MATRIX 
C—VRECVMCMAIN  A,  AUX  A1 ,  AUX  A2, . . .  AUX  AN  : ) 

DO  9001  I JAM-1, NJAMS 
DO  9002  IFR- 1,  NFREQ 

IF (ISLC . EQ . 1 )  VCHA ( 1 , IFR) -VMA JMA ( I JAM , IFR) 

CC  WRITE(6 ,6789) I JAM , IFR , VMA JMA (I JAM , IFR) 

6789  FORMAT  ( IX,  ’  I  JAM ,  IFR,  VMAJMAd  JAM,  IFR)« 1  , 214 , 2X  ,  2E12 . 5) 

DO  9003  IA-1.NAUX 
IAP1-IA+1 

IF(ISLC.EQ.l)  VCHA  (I API , IFR)-VAUXJA(IA ,IJAM, IFR) -VATENA(IA) 
IF(ISLC.EQ.O)  VCHA ( IA , IFR)-VAUXJA(IA , IJAM, IFR)*VATENA(IA) 
WRITE (6 , 8876) I J AM , IFR , IA , VAUX JA ( I A , IJAM , IFR) 

WRITE (8 , 8876) I JAM , IFR , I A , VAUX JA ( IA . IJAM , IFR) 

8876  F0RMAT(1X, ’IJAM, IFR, 1A.VAUXJA-’ ,3I4,2X,2E12.4) 

9003  CONTINUE 
9002  CONTINUE 

DO  5555  KKK-1,4 
DO  5555  LLL-1, NFREQ 

CCC  WRITE (6 , 6655 ) KKK , LLL , VCHA (KKK , LLL ) 

5565  CONTINUE 

WRITEC6 , 5533) 

5533  FORMAT (IX, ’NOW  COMPUTE  JAMMER  COVARIANCE  MATRIX’) 

C---COMPUTE  COVARIANCE  MATRIX  FOR  JTH  JAMMER  SOURCE 

IF(NFREQ.CT.l)  CALL  CO VSWCC VCHA, VCHA, NMAX, NFREQ, BWFHZ , COVAA J ) 
DO  9006  ICH-l.NKAX 
DO  9005  JCH-l.NMAX 

COVNF ( ICH , JCH) -COVNF ( ICH , JCH) ♦COVAA J (ICH , JCH) 

CC  WRITE ( 6 . 2299 ) ICH , JCH , COVNF ( ICH , JCH) 

9005  CONTINUE 
93: 1  CONTINUE 

C---ADD  RECEIVER  NOISE  TO  DIAGONAL  ELEMENTS 
4444  DO  8006  ICH- 1, NMAX 

COVNF (ICH , ICH) -COVNF (ICH , ICH) +1 . ODO 
8006  CONTINUE 

DO  2727  I-l.NMAX 
DO  2727  J-l.NMAX 

IF(CDABS (COVNF (I,J)) .EQ.O.ODO)  GO  TO  2727 
CNDB-10 . -DLOG 10 (CDABS (COVNF ( I , J) ) ) 

PREALN-DREAL ( CO VNF ( I . J ) ) 

PIMAGN-DIMAG (COVNF (I , J) ) 

PHASEN-DATAN2 (PIMAGN . PREALN) /DCR 
WRITE (8 ,4411)1 , J ,CNDB .PHASEN 
IF(I . EQ . 1 ) WRITE (6 ,441 1 ) I , J ,CNDB ,PHASEN 
4411  FORMAT ( IX , ’ I , J ,CNDB .PHASEN- ’, 214, 2X,2F 12-2) 

2727  CONTINUE 
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C« •COMPUTE  COVARIANCE  MATRIX  INVERSE 
DO  8007  ICH-l.NMAX 
DO  8007  JCH- 1. NMAX 
COVNFI ( ICH , JCH ) -CO VNF  C I CH , JCH ) 

8007  CONTINUE 

CALL  DCMINV  (COVNFI  ,LTMP  ,  MTMP ,  NMAX ,  NMAX ) 

C-- -CHECK  MATRIX  INVERSION  ACCURACY  (CINVERSE-C-I) 

CALL  CMMULT (COVNFI ,COVNE, NMAX , NMAX .NMAX .CINVCN) 

DO  8008  ICH-l.NMAX 
)0  8008  JCH- 1, NMAX 

IF ( ICH . EQ . 1 ) WRITE (6 . 8009 ) ICH , JCH , CINVCN ( ICH , JCH ) 

8009  FORMATC IX . * ICH , JCH , CINVCN- > . 214 ,2X , 2E12 . 5) 

8008  CONTINUE 

C— COMPUTE  EIGENVALUES  (ALSO  EIGENVECTORS  AND  PERFORMANCE  INDEX) 

I JOB-2 

CALL  EIGCC (COVNF .NMAX , NMAX , I JOB .EIGVAN .EIGVEN . NMAX .WKE.IER) 
WRITE(6, 155) IER 

1S5  FORMAT (IX, ’AFTER  COMPUTE  EIGENVALUES,  IER-’,I5) 

DO  200  I-l.NMAX 
WRITE (6 ,300)1, EIGVAN ( I ) 

300  FORMAT ( IX , ’ I .EIGVAN- ’ ,2X , 14 ,2E10 . 3) 

200  CONTINUE 

DO  206  I-l.NMAX 

EVLDBN-10 . *DL0G10(CDABS (EIGVAN (I ) ) ) 

WRITE ( 8 , 207 ) I . EVLDBN 
WRITE  ( 6 , 207 )  I ,  EVLDBN 

207  FORMAT ( IX , ’ I .EVLDBN- ’ , 14 . 2X ,2F12 . 3) 

205  CONTINUE 
C---  CALL  PRNTDA 

C— -FILL-IN  SIDELOBE  CANCELLER  QUIESCENT  WEIGHTS 
DO  8985  I-l.NMAX 
WQSLC ( I , 1 ) -DCMPLX ( 0 . ODO . 0 . ODO ) 

8985  CONTINUE 

WQSLC ( 1 . 1 ) -DCMPLX ( 1 . ODO , 0 . ODO ) 

DO  7769  I-l.NMAX 

WRITE (6 , 1212)1 .WQSLC (I , 1) 

WRITE (8 ,1212)1 .WQSLC (I , 1 ) 

1212  FORMAT ( IX , ’ I .WQSLC- ’ ,14 ,2X ,2E12 . 5) 

7769  CONTINUE 
C— COMPUTE  QUIESCENT  I  NR 

IFOSLC.EQ.  1)  CALL  INRTI  0(  WQSLC  .COVNF ,  NMAX  .WTCTR.CMPROD  ,  QINRDB) 
IF(ISLC.EQ.O)  CALL  INRTIO (CWTA, COVNF, NMAX .WTCTR .CMPROD , QINRDB) 

Q INRAA- 1 0 . -DLOG 1 0 ( CD ABS ( CO VNF ( 1 , 1 ) ) ) 

C— COMPUTE  AVERAGE  CANCELLATION 
SUMC-0 . 0 
SUMAA-0 . 0 
DO  1829  IR-l.NRAN 

C— -ZERO-OUT  ADAPTIVE  WEIGHTS  INITIALLY 
DO  57  I-l.NMAX 
WAN ( I , 1 ) -DCMPLX ( 0 . ODO , 0 . ODO) 

57  CONTINUE 

C— COMPUTE  ADAPTIVE  ARRAY  WEIGHTS 

IF(ISLC.EQ.l)  CALL  CMMULT  (CO  VNF  I  .WQSLC, NMAX  .NMAX  ,  1  .WAN) 
IF(ISLC.EQ.O)  CALL  CMMULKCOVNFI  .CWTA  .NMAX  ,  NMAX  ,  1  .WAN) 

C— QUANTIZE  AND  RANDOMIZE  ADAPTIVE  WEIGHT  SETTINGS 
IF( IQUAN . EQ . 1 . AND . NBADWT . LT . 20) 

2CALL  ADQUAN (NMAX , WAN , NBADWT , IRNERR , INITRN, AWSGDB , AWSGDG) 

C— -NORMALIZE  FULLY  ADAPTIVE  WEIGHTS 
IF(ISLC.EQ.O)  THEN 
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SAVE-O. ODO 
DO  33345  I-l.NEL 
SAVE-SAVE-CD ABS ( VAN (1,1)) --2 
33345  CONTINUE 

SAVESR-DSQRT (SAVE) 

DO  44456  I-l.NEL 

WAN (I , 1 ) -WAN ( I , 1 ) /S AVESR 

UANDB-20 . -DL0G10(CDABS (VAN (I , 1) ) ) 

IF(IR.EQ.l)  VRITE(6, 3599)1, WANDB 
IF(IR.EQ.l)  VRITE(8, 3599)1, WANDB 
3599  FORMAT ( IX , > I- ’,14 , 21 , ’ WANDB- > , II ,F12 . 3) 

44456  CONTINUE 
END  IF 

C— -TO  PRINT  NORMALIZED  WEIGHTS 
BIGWDB— 299.0 
DO  9102  I-l.NMAI 

IF(IR.EQ.l)  WRITE(6,6191)I .WAN (I , 1) 

IF(IR.EQ.l)  WRITE (8 ,6191) I ,WAN(I , 1) 

6191  FORMAT ( IX , ’ I ,WAN(I ,1)- ’ , I4,2X,2E12.5) 

WANDBM-20 . -DLOGIO (CDABS (WAN (1,1))) 

IF(WANDBM.GT. BIGWDB)  BIGWDB-WANDBM 
9102  CONTINUE 

DO  9112  I-l.NMAX 

WANDB-20 . -DLOGIO (CDABS (WAN (I , 1) ) ) -BIGWDB 
WANDG-DATAN2 (DIMAG (WAN (1,1)) ,DREAL(WAN (1,1))) /DCR 
IF(IR.EQ.l)  WRITE(6, 35991) I, WANDB, WANDG 
IF(IR.EQ.l)  WRITE(8, 35991) I. WANDB, WANDG 
35991  FORMAT (IX , 1 1- 1 ,14,21, 'WANDB , WANDG- ' , IX ,2F14 . 5) 

9112  CONTINUE 

C-- -COMPUTE  ADAPTED  INR 

CALL  INRT I 0 ( WAN , CO VNF , NMAX , WTCTR , CMPROD , A INRDB ) 

C---COMPUTE  CANCELLATION 

CANCDB-AINRDB-QINRDB 
WRITE! 6 , 3007 ) QINRDB , AINRDB , CANCDB 
WRITE ( 8 , 3007 ) 0 I NRDB , A I NRDB , CANCDB 
3007  FORMAT (IX, 'INR-  QUI.ADAP  .CANCEL- ’ ,2X .3F10 . 3 , 2X , ’  DB’) 
SUMC-SUMC+CANCDB 
C  — -COMPUTE  ADAPTED  INR  CH .  A 
DO  2255  I- 1, NMAX 
WANDMA (1,1) -WAN (1,1) 

IF ( I .GT.NMAX)  WANDMA (I , 1 )-(0. ODO .0 . ODO) 

2255  CONTINUE 

CALL  INRTIO (WANDMA , COVNF , NMAX , WTCTR . CMPROD , AINRAA) 

C-- -CANCELLATION  CH .  A 

CNCLNA-AINRAA-QINRAA 

WRITE (6 ,3738) QINRAA .AINRAA .CNCLNA 

WRITE(8 ,3733) OINRAA .AINRAA , CNCLNA 

3738  FORMATdX,  'SIDELOBE  CANCELLER  CH  .  A  INR-  OUI ,  ADAP  .CAN- ’  ,2X  ,  3F10 . 3) 
SUMAA-SUMAA+CNCLNA 
1829  CONTINUE 

AVECAN-SUMC/NRAN 
WRITE (6 , 4456 ) CANCDB . NRAN , A VEC AN 
4456  FORMAT ( IX , 'CANCDB ,NRAN , AVECAN- ’ ,F12 . 5 ,2X , 15 , 2X ,F12 . 5) 

WRITE ( 8 , 4456 ) CANCDB , NRAN , AVECAN 
AVECAA-SUMAA/NRAN 
WRITE(6 ,2220) A VEC A A 
WRITE (8 , 2220) A VEC A A 

222C  FORMATdX, 'AVE.  CANEL,  CH.  A,  -’.2X.F12.5) 

C  — -SECTION  TO  COMPUTE  ADAPTIVE  ARRAY  RADIATION  PATTERNS 
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CC  IFdNEAR.  EQ . 0 . OR . IANGLP . EQ . 0)  CO  TO  9990 
IF(INEAR. EQ.O)  GO  TO  9990 
BIGADB— 299 . 

IC-0 

DO  8919  IZ-1, NCOLZN 
DO  8919  IY-1, NROWYN 
DO  8919  II-l, NCOLXN 
IC-IC+1 
CSUMA-(0. ,0.) 

IFCISLC.EQ. 1)  CSUMA-CSUMA+DC0NJG(VAN(1 , 1) )-VXMANA(IC) 

DO  7921  IAX-l.NAUX 
IAXP1-IAX+1 

IF ( ISLC . EQ . 1 )  CSUMA-CSUMA+DCON JG ( WAN (I AXP1 , 1 ) ) • VAXCWA (I AX , IC) 

IF (ISLC . EQ . 0)  CSUMA-CSUMA+DCON JG (VAN ( I AX . 1 ) ) -VAXCVA ( IAX . IC) 

7921  CONTINUE 

IF (CABS ( CSUMA ) . EQ . 0 . )  CSUMA«(1 .E-10,0. ) 

PCHADB(  10-20.  «AL0G10(CABS(CSUMA)  ) 

IF(PCHADB(IC) .CT. BIGADB)  BIGADB-PCHADB(IC) 

CC  VRITE(6,4999)IC,PCHADB(IC) 

8919  CONTINUE 

C-- -NORMALIZE  ADAPTIVE  PATTERNS 

IF(IXZ.EQ.l)  VRITE(31. 18889) 

IF ( IXZ . EQ . 1 )  VRITE (31, 18888) XN , NCOLXN , DXN , ZN , NCOLZN , DZN 
IF(IXZ.EQ.l)  WRITE (3 1,4546) 

IC-0 

DO  3459  IZ-1, NCOLZN 
DO  3459  IY-1, NROWYN 
DO  3459  IX-1, NCOLXN 
IC-IC+1 

PCHADB ( IC) -PCHADB ( IC ) -BIGADB 
IF(IPCONN.EQ.O)  WRITE (31, -)IZ, IX .PCHADB (IC) 

IF ( IPCONN .EQ.O)  GO  TO  3459 
IF(IXY.EQ.l)  FH(IY , IX) -PCHADB (IC) 

/F(IXZ.EQ.l)  FMdZ, IX) -PCHADB (IC) 

IFdXZ.EQ.l. AND. IPCONN. EQ.l)  WRITE(31,-)  IZ ,  IX,FK(IZ,  XX) 

3459  CONTINUE 

WRITE  (25,-)  NCOLXN 
DO  1879  IDDD-1. NCOLXN 
WRITE (25,-) PCHADB (IDDD) 

1879  CONTINUE 

IF ( (NCOLXN . GT . 1 . OR . NROVYN . GT . 1 ) . AND . IPLOTN . EQ . 1 . AND . I ANGLP . EQ . 0) 
2CALL  PLOTTR( NCOLXN . NROWYN , NCOLZN , XC , YC , ZC , NPTSN , PCHADB , PCHADB , 
3PCHADB. PCHADB) 

IFdPCONN.  EQ.l.  AND.  IXY.  EQ.l) 

2CALL  PLCONT (NCOLXN , NROWYN , NCLXN2 , NRVYN2 , RX , RY , 

3FM.DMF , -50 . ,10. ,5,1) 

IF (IPCONN  EQ.l. AND . IXZ . EQ . 1 ) 

2CALL  PLCONT (NCOLXN , NCOLZN , NCLXN2 , NCLZN2 , RX , RZ , 

3FM.DMF.-50. ,10. ,5,1) 

9990  CONTINUE 
9999  CONTINUE 

CALL  DONEPL 
CALL  GETCP2 (CPUL) 

CPUT0T-CPUL-CPU1 
WRITE (6, 2006 )CPUTOT 

2006  FORMAT (IX, ’TOTAL  CPU  TIME- ’ , F15 . 2) 

STOP 

END 

C-— SUBROUTINE  TO  COMPUTE  IMPEDANCE  MATRIX 
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SUBROUTINE  ZMATRX (CTDP , STDP , IBLTSL , ICC1 , ICC ,CZ , Z) 

COMPLEX  Z(ICCl.ICC) ,CZ(1) 

COMMON/ A/DX , DY , NCOLX , NROWY , NEL , HZ , HL , ARAD , ZLO AD , ZCHAR 

COMMON  /GROUND/  IGRNDP 

COMMON  /CIRCLE/  ICIRC .RADIUS .HLS 

KB-NCOLX 

IDMB-NROWY 

C—COMPUTE  MUTUAL  IMPEDANCES  Z(1 , 1) ,Z(1 ,2) ,Z(1 .3) . Z(1,NEL). 

CC  CALL  RGDZMN( IGRNDP. ICC1, ICC, Z) 

CC  CALL  RCDZABCCTDP, STDP, IGRNDP, ICC1. ICC. Z) 

IF(ICIRC.EQ.O)  CALL  RGDZA2CCTDP , STDP , IGRNDP , ICC1 , ICC, Z) 

IF (ICIRC . EQ . 1 )  CALL  CADZA2 (CTDP , STDP , IGRNDP , ICC1 , ICC, Z) 
VRITE(8 , 10)  NEL, NCOLX. NROWY 
10  FORMAT (IX , ’NEL .NCOLX .NROWY* ’ ,315) 

ICOUNT-O 

DO  20  1-1, NCOLX 

DO  20  J-l, NROWY 

ICOUNT-ICOUNT+1 

WRITEC6 ,80)  I , J ,Z( 1 . ICOUNT) 

80  FORMATdX,  ’I, J,Z(1, ICOUNT)-'  ,2I4,iX,2E12.S) 

20  CONTINUE 

IDM-NEL 
ICC-NEL 

C---FILL  THE  IMPEDANCE  MATRIX 
IF(NCOLX.LE.l)  GO  TO  70 
IBLT-IBLTSL 
IDM1-IDMB 

IF(IBLT.EQ.O)  IDM1-NEL 
IFtNROWY.GT. 1)  CO  TO  40 
C-.-FILL  TOEPLITZ  MATRIX 
DO  30  1-2, NEL 
DO  30  J-I.NEL 
K-l+J-I 
Z(I , J)*Z( 1 ,K) 

30  CONTINUE 

GO  TO  SO 

40  CALL  BTOEPL ( IBLT , NB , IDMB , IDM1 , IDK , Z ) 

SO  CONTINUE 

IF ( IBLT .  EQ  .  0)  GO  TO  70 
DO  60  I-l.IDMB 
DO  60  J-l.IDM 
IC*( J- 1 ) -IDMB+I 
CZ(IC)-Z(I,J) 

CC  WRITE(8, 7878)1, J.Z(I.J) 

7878  FORMAT ( IX , ’ I , J ,Z(I , J)- ’ ,214 , 2X , 2E12 . 5) 

60  CONTINUE 

70  CONTINUE 

RETURN 
END 

C«* -SUBROUTINE  TO  COMPUTE  RECEIVE  BEAMFORMER  WEIGHTS 

SUBROUTINE  VRCVVT (NACOLX  ,  NAROWY  ,  NXDUM ,  NYDUM ,  EDGTDB  ,  WA  ,  WB) 
DIMENSION  WA( 1 ) ,WB(1) 

DIMENSION  WT(180) 

COMMON/ A  /DX  ,  DY ,  NCOLX  ,  NROWY  ,  NEL  ,  HZ ,  HL ,  ARAD  ,  ZLO  AD ,  ZCHAR 
COMMON  /CHEBY/  ICHEB , SLLDB 
PI-3  141S92654 

tti-z.-pi 

DCR-PI/180. 

CC  SLLDB-40. 
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TAP«10 . (EDCTDB/20 . ) 

DIP09260 

WRITE (6, •) TAP 

DIP09270 

AHP«(l.-TAP)/2. 

DIP09280 

NPDX-NCOLX -NACOLX -2»NXDUM 

DIP09290 

AXL-DI»(NAC0LX-1) 

DIP09300 

AYL*DY • (NAROWY - 1 ) 

DIP09310 

IFCTAP.NE. 1 . 0) FX-AXL/2 . *PI/ACOS  <  TAP) 

DIP09320 

IF(TAP.NE. 1 .0)FY-AYL/2 .*PI/ACOS(TAP) 

DIP09330 

VRITE(6,20) 

DIP09340 

20 

FORMAT (IX, ’BEFORE  CALL  CHEBWT’) 

DIP09350 

WRITEC6 ,30)SLLDB 

DIP09360 

30 

FORMAT ( IX , ’ SLLDB* ’ , 2X , F 10 . 2 ) 

DIP09370 

IF (ICHEB . Eq . 1 ) CALL  CHEBWT ( NACOLX , SLLDB , WT , RLOCS ) 

DIP09380 

VRITE(6 ,40) 

DIP09390 

40 

FORMAT! IX, ’AFTER  CHEBWT’) 

DIP09400 

DO  1000  IC*1,NEL 

DIP09410 

WA(  10*0.0 

DIP09420 

WB(IC)*0 . 0 

DIP09430 

1000 

CONTINUE 

DIP09440 

0*« COMPUTE  EFFECTIVE  DIPOLE  CENTER  COORDS.  FOR  BOTH  PHASE  CENTERS 

DIP09450 

XO— AXL/2. 

DIP09460 

YO—AYL/2. 

DIP09470 

IC*0 

DIP09480 

DO  80  1*1. NACOLX 

DIP09490 

DO  80  J*l, NAROWY 

DIP09S00 

IC-IC+1 

DIP09510 

X*X0+DX*(I-1) 

DIP09520 

Y»Y0+DY*(J-1) 

DIPO9S30 

TAPERX*1 . 0 

DIP09540 

TAPERY*1.0 

DIP09550 

IF(ICHEB.EQ.l)  GO  TO  70 

DIP09560 

IF(TAP.NE.l.O) TAPERX=COS ( P I *  X/FX ) 

DIP09570 

IF (TAP . NE . 1 . 0 . AND . FY . NE . 0 . 0)TAPERY-COS (PI*Y/FY) 

DIP09580 

6C 

TAPER-TAPERI*TAPERY 

DIP09590 

WT(IC)*TAPER 

DIP09600 

70 

CONTINUE 

DIP09610 

80 

CONTINUE 

DIP09620 

<>•.  TRANSFORM  FROM  SUB-APERTURES  TO  FULL  ARRAY 

DIP09630 

ibgna-nxdum.nrowy-nydum 

DIP09640 

IBGNB*IBGNA+NPDX*NROWY 

DIP09650 

IC-0 

DIFC9660 

DO  1010  IX*1, NACOLX 

DIP09670 

DO  1010  IY*1, NAROWY 

DIP09680 

IC-IC+1 

DIP09690 

IA-IBCNA* (IX-1 ) *NROWY+IY 

DIP09700 

IB*IBGNB+(IX-1)*NR0WY+IY 

DIP09710 

WA(IA)-WT(IC) 

DIP09720 

WB(IB)-WT(IC) 

DIP09730 

1010 

CONTINUE 

DIP09740 

RETURN 

DIP09750 

END 

DIP09760 

C*  “SUBROUTINE  TO  COVARIANCE  MATRIX  BASED  ON  NUMERICAL  INTEGRATION 

DIP09770 

C*“IN  THE  FREQUENCY  DOMAIN  ACCORDING  TO  SIMPSON’S  RULE 

DIP09780 

SUBROUTINE  COVSWC(VA ,VB .NCHAN .NFREQ .BWFHZ .COVAB) 

DIP09790 

COMPLEX  *16  VA(NCHAN, NFREQ) ,VB(NCHAN, NFREQ) 

DIP09800 

COMPLEX  *16  COVAB(NCHAN, NCHAN), CSUM.DCABF 

DIP09810 

REAL  *8  SWC(lOl) 

DIP09820 

DELT AF*BWFHZ / ( NFREQ - 1 ) 

DIP09830 

CALL  SIMWC (NFREQ , SWC) 

DIP09840 
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DO  10  ICH-l.NCHAN 
DO  10  JCH-l.NCHAN 
CSUM-CO .CDO ,0 .ODO) 

DO  20  IFR-1 .NFREQ 

DCABF-VA (ICH, IFR) »DCON JG ( VB ( JCH , IFR) ) 

CSUM-CSUM+DCABF*SVC ( IFR) 

CC  WR1TE(6,6767)IFR,CSUM 

6767  FORMATClX, ’IFR.CSUM-’ ,I4,2X,2E12.5) 

20  CONTINUE 

COVABCICH , JCH) -DELTA! /2 . *CSUM 
C«»»NEV  LINE  TO  NORMALIZE  COVAB 

COVABCICH , JCH) -COVAB (ICH , JCH) /BWFHZ 
CC  WRITE (6 , 4567) ICH , JCH , COVAB ( ICH , JCH) 

4567  FORMATClX, ’ICH. JCH, COVABCICH, JCH)-’ ,214, 2X.2E12. 5) 

10  CONTINUE 

RETURN 
END 

C- ••SUBROUTINE  TO  GENERATE  SIMPSON'S  1/3  RULE  WEIGHTING  COEF . 
C»*»INTGERAL  F(X)DX- (DELTAX/3 . ) • (F(1)+4»F(2)+2»F(3)+4»F(4>*  . ,+FfOOD)) 

C**«  THE  SERIES  142424  . 1  ARE  SIMPSON’S  COEF. 

SUBROUTINE  SIKVCCNCOEF ,SWC) 

REAL  »8  SWCCNCOEF) 

DO  10  N-l.NCOEF 
XNN-FLOATCN) 

NN-N/2 
TT-XNN/2 . 

DIF-TT-FLOAT (NN) 

NC-2 

IFCDIF.Eq.O.)  NC=4 
IFCN.EQ. l.OR.N.EQ.NCOEF)  NC-1 
SWCCN)-NC 
10  CONTINUE 
RETURN 
END 

C** •SUBROUTINE  TO  COMPUTE  RECEIVED  VOLTAGES  DUE  TO  JAMMER  SOURCES 
SUBROUTINE  V J AMMR (N JK . NEL . PWRJDB , ICC1 , ICC , PS , CZ , Z , CWTA 
1 ,NB , IDMB , NAUX , IAUXA , ZLOAD , X J AMIN , Y  J AMIN , Z J AMIN , 

2IFR , NFR , CTPR . STPR , CTDP , STDP . VMAINA , VAUXA) 

COMPLEX  «16  CWTACNEL , 1) 

COMPLEX  PS(l) ,CZ(1) ,Z(ICC1,ICC) 

COMPLEX  VJMC180) .VREFJM 

COMPLEX  CJ.CSUMA 

COMPLEX  *16  VMAINA (NJM, NFR) 

COMPLEX  *16  VAUXA (NAUX ,  NJM ,  NFR) 

DIMENSION  PWRJDB(l) .PWRJ(IO) 

DIMENSION  XJAMIN(l) .YJAMIN(l) .ZJAMINC1) 

INTEGER  IAUXA ( 1 ) 

COMMON/PCFWTR/PCDXIN 

COMMON  / B/  NGEN.IGEN, THETAS, PHIS, IMUT.IBLTSL.IPATRN 

COMMON  /D/  FGHZ , RL AMD A , IWL , IS , NSCANS 

COMMON  /WRITE/  IWR 

COMMON  /CIRCLE/  ICIRC .RADIUS ,HLS 

PI-3.141592654 

DCR-PI/180 . 

CJ-(0. ,1.) 

CINMTR-0 . 0254 
ISYM-0 

C« ••CONVERT  DB  TO  POWER  (RELATIVE  TO  NOISE) 

DO  10  II-l.NJK 


DIP098S0 
DIP09860 
DIP09870 
DIP09880 
DIP09890 
DIP09900 
DIP09910 
DIP09920 
DIP09930 
DIf 09940 
DIP09960 
DIP09960 
DIP09970 
DIP09980 
DIP09990 
DIP10000 
DIP10010 
DIP10020 
DTP1003C 
DIP10040 
DIP10050 
DIP10060 
DIP10070 
DIP10080 
DIP10090 
DIP10100 
DIP10110 
DIP10120 
DIP10130 
DIP10140 
DIP10150 
DIP10160 
DIP10170 
DIP101S0 
DIP10190 
DIP10200 
DIP10210 
DIP10220 
DIP10230 
DIP10240 
DIP10250 
DIP10260 
DIP10270 
DIP10280 
DIP10290 
DIP10300 
DIP10310 
DIP10320 
DIP10330 
DIP10340 
DIP10350 

DIP10360 

DIP10370 

DIP10380 

DIP10390 

DIP10400 

DIP10410 
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PWRJ(II)-10 . --(PWRJDBOD/IO . ) 

WRITE  (6 , 66)  II ,  PWRJ  ( II ) 

WRITE(8, 66)11, PWJU(II) 

66  F0RMAT(1X , ’ II ,PWRJ  IN  POWER- ’ , 14 ,21 ,E12 . 5) 

10  CONTINUE 

IF(IFR.CT.l)  GO  TO  65 
VRIT£(6 ,20) IWL 
DO  6544  IJAM-l.NJM 
WRITEC6  40) 

WRITE(6,50)XJAMIN(IJAM) ,YJAMIN(I JAM) ,ZJAMIN(I JAM) 

5644  CONTINUE 
20  F0RMAT(1X.’IVL«’,I4) 

WRITEC6,30) 

30  FORMAT (IX, ’CHANGE  NEAR  FIELD  JAMMER  POSITIONS  (INCHES)?,  ICh-1’ 
READ(6,-)ICH 
IF(ICH.EQ.O)  CO  TO  55 
DO  8887  IJAM-l.NJM 
WRITE(6 ,40) 

40  FORMAT (IX, ’XJAMTN. YJ4MIV  ZIL'IIV  =  ’) 

READ(5  ,*  ,  END-55)  XJAMIN(  I  JAM)  .YJAMIN(IJAM)  ,ZJAMIN(I  JAM) 

WRITE (6, 50) XJ AMIN (I JAM) , YJAMIN(IJAM) ,ZJAMIN(IJAM) 

60  FORMAT ( II, 6F10. 3, 21,215) 

8887  CONTINUE 
55  CONTINUE 

PCDX-PCDXIN-CINMTR 
WRITE (6 ,4757) IFR 

4757  FORMAT (IX, ’FREQ-  INDEX,  IFR-’, 14) 

60  CONTINUE 

C--ALL  DIMENSIONS  IN  METERS 
DO  180  IPHACN-1,1 
WRITE (6 ,2223) IPHACN 
2223  FORMAT (IX, ’IPHACN-’ ,14) 

IREFDP— PCDX/2- *(IPHACN-1)»PCDX 
C---PERFORM  JAMMER  SOURCE  SCAN 
DO  180  IJAM-l.NJM 

X-X JAMIN( I JAM) •CINMTR+PCDX- ( IPHACN- 1 ) 

Y-YJAMIN(IJAM)*CINMTR 
ZPOS-ZJAMIN(I JAM) -CINMTR 
WRITE ( 6 , 6688 ) X , Y , ZPOS , IREFDP 
6688  FORMAT ( IX, 'X.Y. ZPOS. IREFDP  (METERS)-’ ,4F12 .4) 

IF(ICIRC.EQ.O) 

2CALL  NFDPX2(CTPR,STPR.CTDP.STDP, X.Y, ZPOS, IREFDP, 1 , VJM , VREFJM) 
IF(ICIRC.EQ.l) 

2CALL  NFDPC2 (X.Y, ZPOS, IREFDP, 1, VJM, VREFJM) 

C---NORMALIZE  INCIDENT  JAMMER  POWER 
DO  70  INORM-1 ,NEL 

VJM ( INORM) -VJM ( INORM) / VREFJM-SQRT ( PWRJ ( I JAM) ) 

WRITE (6 , 77 ) INORM , V JM (INORM) 

77  FORMAT (IX, ’INORM, VJM (I NORM)  VOLTAGE- ’, 14 ,2X , 2E12 . 5) 

70  CONTINUE 

IF(IMUT.EQ.O)  GO  TO  90 

C—SOLVE  EACH  SYSTEM  OF  EQUATIONS  FOR  THE  UNKNOWN  CURRENTS 
IF(IBLTSL.EQ.l)  GO  TO  80 
112-1 

IFUJAM.GT.l.  OR.  IPHACN.  GT.l)  112-2 
VRITE(6,6110) 

6110  FORMAT (IX, ’CALL  CROUT  IN  VJAMKER 1 ) 

CALL  CROUT (Z, VJM, ICC1, ICC, ISYM.IWR, 112, NEL) 

GO  TO  90 


DIP10420 

DIP10430 

DIP10440 

DIP10450 

DIP10460 

DIP10470 

DIP1048C 

DIP10490 

DIP10500 

DIP10610 

DIP10520 

DIP10530 

DIP10550 

DIP10560 

DIP10570 

DIP10580 

DIP10590 

DIP10600 

DIP10610 

DIP10620 

DIP10630 

DIP10640 

DIP10660 

DIP10670 

DIP10680 

DIP10700 

DIP10710 

DIP10720 

DIP10730 

DIP10740 

DIP10750 


DIP10790 


DIP10820 

DIP10830 

DIP10840 

DIP10850 

DIP10860 

DIP10870 

DIP10880 

DIP10890 

DIP1O900 

DIP10910 

DIP10920 

DIP10930 

DIP10940 

DIP10950 

DIP10960 
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80  I ENTRY-4 

IF (X JAN . GT  1 . OR. IPHACN. GT . 1)  I ENTRY-3 
CALL  BLTSOL(CZ , VJM.PS .NB , IDMB .IENTRY) 

IF(NEL.LT.40)  CALL  CNORM(VJM.NEL) 

90  CONTINUE 

C---COMPUTE  RECEIVED  VOLTAGES  FOR  PRESENT  SOURCE  POSITION 
DO  100  IEL-l.NEL 
VJM(IEL)-VJMCIEL) -ZLOAD 
CC  WRITE (6 ,4456) IEL, VJM(IEL) 

4456  FORMATdX, ’IEL,VJH(IEL)  RECEIVED  VOLT.-’ ,I4,2X,2E12.S) 
100  CONTINUE 

C-.-PERFORM  BEAN  FORMATION  FOR  MAIN  A 
C— PHASE  CENTER  A 

CSUNA-(O.ODO.O.ODO) 

DO  110  KC-l.NEL 
CSUMA-CSUMA+VJMCKC) -CWTA(KC ,1) 

110  CONTINUE 

VNAINA ( I JAM , IFR) -CSUMA 

WRITE (6 , 2222) I JAM , IFR , VMAINA ( I J AM , IFR) 

2222  FORMAT (IX, ’I JAM, IFR, VMAINA-  (AFTER  8 . F . ) ’ ,214 ,2X ,2E12 . 5) 
C... COMPUTE  AUXILIARY  CHANNEL  VOLTAGES 
DO  7000  IAUX-l.NAUX 
VAUXA ( IAUX , I J AM , IFR) -VJM(IAUXA( I AUX) ) 

7000  CONTINUE 
180  CONTINUE 

DO  9000  IJAM-1 .NJAMS 

CCC  WRITE(6 ,4433) I JAM , IFR, VMAINA (I JAM , IFR) 

4433  FORMAT (IX , ’ VJAMMR : I JAM, IFR, VMAINA- ’ , 214 , 2X , 2E12 . 5) 

DO  9001  IAUX-1 ,NAUX 

WRITEC6 , 3333) I AUX , I JAM , IFR , VAUXA ( IAUX , I JAM , IFR) 

3333  FORMAT ( IX , ’VJAMMR : IAUX , I JAM, IFR, VAUXA- 1 ,314 ,2X ,2E12 . 5) 
9001  CONTINUE 
9000  CONTINUE 
RETURN- 
END 

C—SUBROUTINE  TO  COMPUTE  INTERFERENCE  TO  NOISE  RATIO 
SUBROUTINE  INRTIO(WT ,COV ,NEL , WTCTR . CMPROD .DBINR) 

IMPLICIT  REAL  *8  (A-H.O-Z) 

COMPLEX  -16  WT(NEL.l) ,COV(NEL,NEL) ,WTCTR(1 ,NEL) 

COMPLEX  -16  CMPROD ( 1 ,NEL) ,CPR0D1 (1,1) ,CPR0D2( 1 ,1) ,CINR 
COMMON  /NORMAL/  INRNOR 
WRITE (12, 2222) 

2222  FORMATdX,  ’INSIDE  INRTIO  SUBROUTINE’) 

DO  110  I-l.NEL 
CC  WRITE(12, 445)1, WT(I,1) 

445  FORMATdX,  ’I,WT(  1,1)-’ , 2X ,  14 ,2X ,2E12 . 5) 

110  CONTINUE 

CALL  CONJTR(WT,NEL.l, WTCTR) 

DO  111  I-l.NEL 

CC  WRITE ( 12 ,666) I , WT(I , 1) , WTCTR (l , I) 

666  FORMATdX,  'I.  WT,  WTCTR-'  ,2X,  14, 2X.4E12. 5) 

111  CONTINUE 

CALL  CMMULTt WTCTR, COV, 1 .NEL.NEL, CMPROD) 

-ALL  C"MUtTf CMPROD, WT,  1  ,NEL ,  1  .CPR0D1) 

CALL  CMMULT(WTCTR ,  GT , 1 ,  NEL , 1 , CPR0D2) 

WRITE  ( 6 , 333 )  CPROD 1  (1 , 1 )  ,  CPR0D2  (1.1) 

WRITE(8 ,333)CPR0D1 (1 , 1 ) ,CPR0D2(1 , 1) 

333  FORMAT (1 X , ’ CPROD 1 , CPR0D2- ’ , 4E 1 2 . 5 ) 

C---INR  NORMALIZED 


DIP10970 
DIP10980 
DIP 10990 
DIPllOOO 
DIP11010 
DIP11020 
DIP11030 
DIP11040 
DIP11050 
DIP1106 j 
DIP11070 
DIP11080 
DIP11090 
DIPlllOO 
DIP11110 
DIP11120 
DIPU130 
DIP11140 
DIP11150 
DIP11160 
DIPU170 
DIP11180 
DIP11190 
DIP11200 
DIP11210 
DIP11220 
DIP11230 
DIP11240 
DIP11250 
DIP11260 
DIP11270 
DIP11280 
DIP11290 
DIP11300 
DIP11310 
DIP11320 
DIP11330 
DIP11340 
DIP11350 
DIP11360 

DIP11370 

DIP11380 

DIP11390 

DIP11400 

DIP11410 

DIP11420 

DIP11430 

DIP11440 

DIP11450 

DIP11460 

DIP11470 

DIP11480 

DIP11490 

DIP11500 


DIP11520 

DIP11530 
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cc 

INRNOR-O 

WRITE (6, 7739) INRNOR 

DIP11550 

VRITE(8.7739)INRN0R 

DIP11660 

7739 

FORMAT (IX, ’INR  NORMALIZATION  PARAMETER,  INRNOR-’ ,14) 

DIP11670 

IF ( INRNOR . EQ . 1 )  CINR-CPROD 1(1,1) /CPR0D2 (1,1) 

DIP11680 

C***INR  SOT  NORMALIZED  FOR  INRNOR-O 

IF ( INRNOR. EQ.O)  CINR-CPR0D1 (1 , 1) 

DIP11600 

DBINR-10. *DIOG10(CDABS(CINR) ) 

DIP11610 

RETURN 

DIP11620 

END 

DIP11630 

SUBROUTINE  CMMULT(A,B,L,M,N,C) 

DIP11640 

COMPLEX  *16  A(L,M) ,B(M,N) ,C(L,N) 

DIP11650 

DO  20  I-l.L 

DIP11660 

DO  20  J-l.N 

DIP11670 

C(I,J)-DCMPLX(0.0D0,0.0D0) 

DIP11680 

DO  20  K-l.M 

DIP11690 

C(I,J)*C(I,J)+(A(I,K)»B(K,J)) 

DIP11700 

20 

CONTINUE 

DIP11710 

RETURN 

DIP11720 

END 

DIP11730 

C***SUBROUTINE  TO  COMPUTE  CONJUGATE  TRANSPOSE  OF  A  MATRIX 

DIP11740 

SUBROUTINE  CON JTR ( A , L , M , ACTR) 

DIP11750 

COMPLEX  *16  A(L,M) ,ACTR(M,L) 

DIP11760 

DO  10  1*1, M 

DIP11770 

DO  10  J«1,L 

DIP11780 

ACTR(1 , J)“DCONJG(A( J , I) ) 

DIP11790 

CC 

VRITE(12 ,7766) I , J , A( J , I) , ACTR (I , J) 

DIPU800 

7766 

FORMAT ( IX , ' I , J , A ( J  I' , ACTR(I , J)*’ , 2X , 214 , 2X ,4E12. 5) 

DIPU810 

10 

CONTINUE 

DIP11820 

RETURN 

DIP11830 

END 

DIP11840 

SUBROUTINE  DCMINV ( A , L , M , IDM , NEQ) 

DIP11850 

COMPLEX  *16  A(IDM. IDM) .BIGA, HOLD. DET 

DIP11860 

INTEGER  L(IDM) ,M(IDM) 

DIP11870 

N*NEQ 

DIP11880 

DET-DCMPLX ( 1 . ODO , 0 . ODO ) 

DIP11890 

DO  80  K-l.N 

DIP 11900 

L(K)«K 

DIP11910 

M(K)«K 

DIP11920 

BIGA*A(K,K) 

DIP11930 

DO  20  J-K.N 

DIP11940 

DO  20  I-K,N 

DIP11950 

10 

IF (CDABS(BIGA) -CDABS (A(I ,J)))15,19,19 

DIP11960 

15 

BIGA*A(I,J) 

DIP11970 

LOO-I 

DIP11980 

M(K)*J 

DIP11990 

19 

CONTINUE 

DIP12000 

20 

CONTINUE 

DIP12010 

J-L(K) 

DIP12020 

IF(J-K)35,35,25 

DIP12030 

25 

CONTINUE 

DIP12040 

DO  30  I-l.N 

DIP12050 

HOLD— A(K.I) 

DIP12060 

A(K.I) 'A(J.I) 

DIP12070 

30 

A(J,I)-HOLD 

DIP12080 

35 

I-M(K) 

DIP12090 

IF(I-K)  45,45,38 

DIP12100 

38 

CONTINUE 

DIP12110 

Do  J*1,N 

DIP12120 

90 


T 


H0LD«-A(J,K) 

DIP12130 

AC J ,K)-AC J , I) 

DIP12140 

40 

AC J , I) -HOLD 

DIP12150 

45 

CONTINUE 

DIP12160 

DO  55  I-l.N 

DIP12170 

IF(I-K)50, 55.50 

DIP12180 

60 

ACI,K)-A(I,K)/C-BIGA) 

DIP12190 

65 

CONTINUE 

DIP12200 

DO  65  I-l.N 

DIP12210 

DO  65  J-l.N 

DIP12220 

IF  ( I -K)  60, 64 ,60 

DIP12230 

60 

IF(J-K)62,64,62 

DIP12240 

62 

ACI,J)*A(I,K)*A(K,J)+A(I 

,J) 

DIP12250 

64 

CONTINUE 

DIP12260 

65 

CONTINUE 

DIP1227C 

DO  75  J»1,N 

DIP12280 

IF(J-K)70,75,70 

DIP12290 

70 

A(K,J)=A(K,J)/BIGA 

DIP12300 

75 

CONTINUE 

DIP12310 

DET-DET-BIGA 

DIP12320 

A(K,K)=1 .ODO/BIGA 

DIP12330 

80 

CONTINUE 

DIP12340 

K-N 

DIP12350 

100 

K-K-l 

DIP12360 

IF  00  150 , ISO , 105 

DIP12370 

105 

I-LCK) 

DIP12380 

IF(I-K)120,120,108 

DIP12390 

108 

CONTINUE 

DIP12400 

DO  110  J-l.N 

DIP12410 

HOLD-AU.K) 

DIP12420 

AC  J  ,K.'“-AC  J ,  I) 

DIP12430 

110 

AC J , I)*HOLD 

DIP12440 

120 

J-MCK) 

DIP12450 

IF ( J-K) 100,100,125 

DIP12460 

125 

CONTINUE 

DIP12470 

DO  130  I-l.N 

DIP12480 

HOLDBACK, I) 

DIP12490 

ACK.I)—  AC  J  ,1) 

DIP12500 

130 

A CJ, I) -HOLD 

DIP12510 

GO  TO  100 

DIP12520 

150 

RETURN 

DIP12530 

END 

DIP12540 

C-*- SUBROUTINE  CNORM 

DIP12550 

0»*[ 

COMPUTES  A  NORMALIZED  COMPLEX 

COLUMN  VECTOR  SUCH  THAT  THE 

DIP12560 

C---MAXIMUM  ELEMENT  HAS  UNITY  MAGNITUDE. 

DIP12570 

C 

DIP12580 

C---PRINTS  MAGNITUDE  AND  PHASE 

OF 

NORMALIZED  VECTOR 

DIP12590 

c 

DIP12600 

SUBROUTINE  CNORMCV.N) 

DIP12610 

c*** 

V  IS  THE  INPUT  COMPLEX 

COLUMN  VECTOR 

DIP12620 

c*** 

N  IS  THE  LENGTH  OF  V 

DIP12630 

COMPLEX  VCD  ,SS 

DIP12640 

CNOR-O . 0 

DIP12650 

DO  10  K-l.N 

DIP12660 

SA-CABSCVCK) ) 

DIP12670 

IFCSA.GT  CNOR)  CNOR-SA 

DIP12680 

10 

CONTINUE 

DIP12690 

IFCCNOR.LE.O. )  CNOR- 1.0 

DIP12700 

DO  30  K-l ,N 

DIP12710 
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SS-V(K) 

SA-CABS(SS) 

SNOR-SA/CNOR 

PHR-O. 

IF(SA.GT.O.)  PHR-ATAN2 (AIMAG (SS) ,REAL(SS) ) 

PH-57. 29578-PHR 
WRITE(8,20)  K,SNOR,SA,PH 
20  FORMAT (IX, 15 ,F10 . 6,3X,E15 . 3.F10 .0) 

30  CONTINUE 
RETURN 
ENG 

SUBROUTINE  GETCP2(RCPU) 

REAL  ETIKE,TARRAY(2) ,RCPU 
TIME-ETIME (TARRAY ) 

RCPU-T ARRAY ( 1 ) 

RETURN 

ENG 

C-* -SUBROUTINE  TO  COMPUTE  INGUCED  VOLTAGE  BETWEEN  PROBE  ANG 
C-. -DIPOLE  ARRAY  ELEMENTS 
C— *  RING  ARRAY 

SUBROUTINE  NFDPC2  (XP  ,  YP  ,ZP  ,XREF ,  IREF .  V  ,  VREF) 

C— DISTANCES  ARE  IN  METERS 

COMPLEX  V(l) ,VD,VR1,VR£F 
COMPLEX  -16  D2ABG 

COMMON/A/DX , DY , NCOLX , NROWY .NEL.HZ.HL, ARAD , ZLO AD , ZCHAR 
COMMON  /CIRCLE/  ICIRC .RADIUS ,KLS 
C— NOTE:  AA.BB,...  FOR  PROBE  11,22....  FOR  DIPOLE 
PI-3.1415926535 
DCR-PI/180. 

C-— PROBE  DIMENSIONS 
XAA-XP 
IBB-XP 
XCC-XP 
YAA-YF-HLS 
YBB-YP 
YCC-YP+HLS 
785  ZAA-ZP 
ZBB-ZP 
ZCC-ZP 
IC-0 

C— — CIRCULAR  DIPOLE  ARRAY  ELEMENTS 
DELPHI-360 . /NEL 
DO  10  IX-l.NEL 
PHI -DELPHI - (IX- 1) 

ID-RADIUS-COS(PHI-DCR) 

ZD-RADIUS-SIN (PHI-DCR) 

DIST«SQRT((XP-XD)  —  2+(ZP-ZD)—  2) 

IF(DIST.LT.ARAD)  XD-XD+ARAD 
Xll-XD 
I22-XD 
X33-XD 
Yll-YP-HL 
Y22-YP 
Y33-YP+HL 
Zll-ZD 
Z22-ZD 
Z33-ZC 
9446  IC-IC+1 

CC  WRITE(6 ,2233)  IX ,  IY  ,XD ,  YD 


DIP12720 

DIP12730 

DIP12740 

DIP12750 

DIP12760 

DIP12770 

DIP12780 

DIP12790 

DIP12800 

DIP12810 

DIP12820 


YDI00010 


YDI00060 

YDI00070 

YDI00120 


YDI00260 

YDI00270 

YDI00280 

YDI00300 


YDI00330 

YDI00430 


YDI00490 


YD 100650 
YDI00660 


YD 100670 


2233  FORMAT (IX ,  ’ IX , IY ,XD , YD* ’,21,214, 2X ,2F12 .6) 

CC  WRITE (6 , 7854) XAA , IBB , ICC , YAA , YBB , YCC , ZAA , ZBB , ZCC 

7854  F0RMAT(1X, ’XYZABC-’ .9E12.5) 

CC  WRITE (6, 7855) XI 1 ,X22 .X33 , Yll , Y22 , Y33 ,Z11 ,Z22 ,Z33 

7855  FORMAT ( IX , ’XYZ123- ’ ,9E12 . 5) 

CALL  DSZABC ( XAA , XBB , ICC , Y AA , YBB , YCC , ZAA , ZBB , ZCC . XI 1 , X22 , X33 , 
2Y 1 1 , Y22 , Y33 , Z 1 1 , Z22 , Z33 , DZABG ) 

VD-DZABG 

CC  WRITE(6,33)IC,VD 

33  FORMAT ( II , ’  IC , VD- ’ ,2X , X5,2E12 .5) 

V(IC)«VD 

VAMPDB-  20 . *AL0G10(CABS(V (IC) ) ) 

VPKASE«ATAN2(AIMAG( V(IC) ) ,REAL(V(IC) ) )*180 . /3 . 141592654 
CC  WRITE  (6 , 4455)  IC ,  VAMPDB ,  VPHASE 

4455  FORMAT (IX, ’ IC ,VAMDB .VPHASE* ’ ,2X , 14 ,2X ,2F12 .2) 

1C  CONTINUE 

C---COMPUTE  REFERENCE  VOLTAGE  (FICTITIOUS  ELEMENT  AT  XREF,  Y-0) 
C---SET  XREF-RADIUS,  THIS  MAKES  ELEMENT  1  THE  REFERENCE 
IREF-RADIUS 

C»**SET  IREF-ARAD,  (REF.  CLOSE  TO  THE  ORIGIN) 

XREF-ARAD 

Xll-XREF 

X22-XREF 

X33-XREF 

Yll— HL 

Y22-0.0 

Y33=HL 

Zll-0.0 

Z22-0 . 0 

Z33-0.0 

9267  CALL  DSZABG (XAA , XBB , ICC , YAA , YBB , YCC , ZAA , ZBB , ZCC , XI 1 ,X22,X33 , 
2Y11 , Y22 , Y33 , Zll , Z22 , Z33 , DZABG) 

VR1-DZABG 
VREF-VR1 
99  RETURN 
ENL 

C-* -SUBROUTINE  TO  COMPUTE  MUT.  IMPED.  BETWEEN  STRAIGHT  DIPOLES 
C---ARRANGED  IN  A  RING  ARRAY  (CIRCLE) 

SUBROUTINE  CADZA2 (CTDP , STDP , IGRNDP , ICC1 , ICC , Z) 

COMPLEX  ZMA , ZABG , Z(ICC1 , ICC) 

COMMON  /A/  DX , DY , NCOLX . NROWY , NEL , HZ , HL , ARAD , ZLOAD , ZCHAR 
COMMON  /CIRCLE/  ICIRC, RADIUS, HLS 
C— -ALL  DIMENSIONS  IN  METERS 
PI-3 . 1415926535 
DCR-PI/180 . 

C---FIXED  POSITION  FOR  ELEMENT  1 
IAA-RADIUS 
IBB-RADIUS 
XCC-RADIUS 
YAA— HL 
YBB-0.0 
YCC-HL 
ZAA-0.0 
ZBB-0.0 
ZCC-O. 0 
YU-YAA 
Y22-YBB 
Y33-YCC 


YDI00680 

YDI00690 

YDI00700 


YDI00780 

YDI00790 

YDI00800 

YDI00810 

YDI00820 

YDI00830 


YDI00870 


YDI00980 


YDI01000 

YDI01010 

YDI01020 

YDI01080 

YDI01090 
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C— VARIABLE  POSITION  FOR  RING  ARRAY  ELEMENTS 
DELPHI -360 . /NEL 
IC-0 

DO  20  I-l.NEL 
IC-IC+1 

PHI -DELPH I*(I-1)-DCR 
ID-RADIUS-COS  (PHI ) 

ZD- RADI’JS-SIN  (PHI ) 

111-ID 

I22-XD 

X33-XD 

IF(I.EQ.l)  I11-X11+ARAD 

IF(I.EQ.l)  I22-X22-ARAD 

IF(I.EQ.l)  I33-I33+ARAD 

Zll-ZD 

Z22-ZD 

Z33-ZD 

CC  VRITE(6,87)IC 

87  F0RMAT(1X , ' IC- ’ , 15) 

CC  WRITE ( 6 , 88 ) IAA , XBB , XCC , YAA , YBB , YCC , Z AA , ZBB , ZCC 

88  FORMAT (IX, ’XYZABC- ’ ,9E12 .4) 

CC  WRITE (6 , 89)111, X22,X33,Y11,Y22,Y33,Z11,Z22,Z33 

89  F0RMAT(1X, ’XYZ123-’ .9E12.4) 

ZMA-ZABG ( XA ‘ , XBB , XCC , Y AA , YBB , YCC , ZA A , ZBB , ZCC , X 1 1 , 122 , X33 , 

2Y1 1 , Y22 , Y33 ,Z1 1 , Z22 , Z33) 

Z(1,IC)-ZMA 

IF(IC.EQ.l)  Zd,l)-Z(l,l)+ZLOAD 
20  CONTINUE 

DO  40  I-l.NEL 

IF ( I . LT . 9) WRITE (6,30)I,Z(1,I) 

WRITE (8 ,30) I , Z(1 , I) 

30  FORMAT (IX , 1 Z(1 , ’ , 14 , ■ )- ■ ,2E12 . 5) 

40  CONTINUE 

RETURN 
END 

•  «**»••••{ il*  zabgenlosa  .f  —  —  — 

C— PROGRAM  TO  CALCULATE  MUTUAL  IMPEDANCE  BETWEEN  TWO  DIPOLES 
C---WITH  ARBITRARY  LENGTH  AND  ORIENTATION.  A  PIECEWISE- 
C— SINUSOIDAL  CURRENT  DISTRIBUTION  IS  ASSUMED. 

COMPLEX  FUNCTION  ZABGUl  ,X2 ,X3,  Y1 ,  Y2 , Y3 ,Z1  ,Z2  ,Z3  ,XA,XB ,XC , YA , YB 
2,ZA,ZB,ZC) 

COMPLEX  Pll,P12,P21,P22,011,Q12,021,522,Rll,R12,R2l,R22 
COMPLEX  SI 1 , S12 , S21 , S22 , JCOM , GAM , CGDS , SGDS , SGDT , ETA ,EP3 
COMPLEX  EGDS.EGDT 
COMMON  /?/  FHZ ,ER3 ,SIG3 ,TD3 
C— ALL  DIMENSIONS  IN  METERS 
PI-3 . 141592654 
TPI-2.-PI 
B-TPI 

JCOM- (0. ,1.) 

E0-8.8S4E-12 
U0-1.2566E  3 
OMEGA-TPI-FHZ 

IF(SIG3 . LT . 0 . )EP3-ER3*E0»CMPLX (1 . , -TD3) 

IF(TD3 .LT .0 . )EP3«CMPLX(ER3*E0 , -SIG3/0MEGA) 

ETA-CSQRT (U0/EP3 ) 

GAM-0MEGA«CSQRT(-U0-EP3) 


ZAB00010 

ZAB00020 

ZAB00030 

YCZAB00040 

ZAB00050 

ZAB00060 

ZAB00070 


ZAB00090 


WIR01670 

WIR01720 

WIR01730 

WIR01740 

WIR01760 
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AM-0.0001 

IF  (CABS  (GAM*  AH)  .CT.0.06)  WRITE(6 ,7923)  AM 
7923  FQRMAT(1X,  ’CABS  (GAM*  AM)  IS  GREATER  THAN  0.06,  AM-',E14.5) 
INT-0 
XBA-XB-XA 
YBA-YB-YA 
ZBA-ZB-ZA 
X21-I2-X1 
Y21«Y2-Y1 
Z21-Z2-Z1 

DS-SQRT(XBA*XBA+YBA*YBA+ZBA*ZBA) 

DT-SQRT (I2,.*I21+Y21*Y21*Z21*Z2l) 

CC  DSK«B*DS 

CC  DTK-B*DT 

CC  CGDS«CMPLX(COS(DSK) ,0.0) 

CC  SGDS-CMPLX (0 . O.SINCDSK) ) 

CC  SGDT-CMT  '.(0.0,  SIN  (DTK) ) 

C***F0R  LOSS)  MEDIUM  THE  NEXT  LINES  ARE  APPROPRIATE 
EGDS«CEXP(GAM*DS) 

EGDT-CEXP(GAM*DT) 

CGDS- (EGDS+1 . /EGDS) /2 . 

SGDS-CEGDS- 1 . /EGDT) /2 . 

SGDT-(EGDT-l./EGDT)/2. 

CC  WRITE(6,1345)XA,XB,YA,YB,ZA,ZB 

1345  FORMAT ( IX , ’XYZAB- 1 ,6E14 .5) 

CC  WRITE (6 ,1346)X1,X2,Y1,Y2,Z1,Z2 

1346  FORMAT ( IX , ’XYZ12- ’ ,6E14 . 5) 

CALL  GGS(XA,YA,ZA,XB,YB,ZB,X1,Y1 ,Z1 ,X2,Y2,Z2,AM, 

2DS , CGDS , SGDS , DT , SGDT , INT , ETA , GAM , PI 1 ,P12 , P21 , P22) 

CALL  GGS (XA , YA ,ZA ,XE , YB ,ZB ,X2 , Y2 , Z2 ,X3 , Y3 , Z3 , AM , 

2DS, CGDS,  SGDS. DT,  SGDT,  IKT,  ETA,  GAM.  QU,Q12,Q21,Q22) 

CALL  GGSCXB , YB , ZB ,XC , YC ,ZC ,X1 , Y1 ,Z1 ,X2 , Y2 ,Z2 , AM , 

2DS, CGDS, SGDS, DT. SGDT, TNT, ETA, GAM, R1 1 ,Rl2 ,R21 ,R22) 

CALL  GGS (XB , YB .ZB ,XC , YC ,ZC ,X2 , Y2 , Z2 ,X3, Y3 ,Z3 , AM, 

2DS, CGDS, SGDS. DT, SGDT, INT, ETA, GAM, Sll .S12.S21 ,S22) 
ZABG-P22+Q21-R12-S11 
CC  WRITE (6 ,7898) ZABG 

7898  FORMAT (IX, ’EXITING  ZABG  WITH  ZABG- ’ , 2E14 . 5) 

nnjR« 

END 

•  He  dzabgnloss . f ******* 

C****DOUBLE  PRECISION  VERSION 

C... PROGRAM  TO  CALCULATE  MUTUAL  IMPEDANCE  BETWEEN  TWO  DIPOLES 
C**»WITH  ARBITRARY  LENGTH  AND  ORIENTATION.  A  PIECEWISE- 
C***SINUSOIDAL  CURRENT  DISTRIBUTION  IS  ASSUMED. 

SUBROUTINE  DSZABG (SX1 , SX2 , SX3 , SY 1 , 3Y2 , SY3 , SZ 1 , SZ2 , SZ3 , 

2SXA ,SXB ,SXC ,SYA ,SYB , Sv~ ,SZA ,SZB ,SZC ,DZABG) 

IMPLICIT  REAL*8  (A-N.O-Z) 

COMPLEX* 16  P11.P12.P21 ,P22 ,011 , Q12 , Q21 ,022 ,R11,R12,R21 ,R22 

COMPLEX *16  Sll, S12.S21.S22, 3COM, GAM, CGDS , SGDS , SGDT , ETA 

COMPLEX* 16  DZABG.EP3, EGDS, EGDT 

REAL*4  SX1 , SX2 . SX3.SY1 , SY2 , SY3.SZ1 , SZ2 ,SZ3 

REAL*4  SX A , SXB , SXC , S YA , S YB , SYC , SZA , SZB , SZC 

REAL*4  FHZ,ER3,SIG3,TD3 

COMMON  /F/  FHZ,ER3 ,SIG3 ,TD3 

JC0M*(0 . DO , 1 . DO) 

PI*3 . 1415926535898D0 
TPI*2.0D0*PI 


ZAB00130 

ZAB00140 

ZABOOISO 

ZAB00160 

ZAB00170 

ZAB00180 

ZAB00190 

ZAB0D200 

ZAB00210 

ZAB00220 

ZAB00230 

ZAB00240 

ZAB00250 

ZAB00260 


ZAB00270 
ZAB00280 
ZAB00290 
ZAB00300 
ZAB00310 
ZAB00320 
ZAB 00330 
ZAB00340 
ZAB00350 


ZAB00360 

ZAB00370 


DZA00020 

DZA00030 

DZA00040 

DZA00050 

DZA00060 

DZA00070 

DZA00080 

DZA00090 

DZA00110 

DZA00120 


DZA00140 


B-TPI 


EO-8.864D-12 

DO-1.2666D-6 

OMEGA*TPI »FHZ 

WIR01670 

cc 

WRITE (6 ,2843) OMEGA 

WIR01680 

2843 

FORMAT (11, ’OMEGA-’ ,E12. 6) 

W1R01690 

IFCSIG3  LT .0 .ODO)EP3*ER3*EO»DCMPLX( 1 .ODO , -TD3) 

WIR01720 

IF  ( TD3 .  L  T .  0 .  ODO )  EP3*DCMPLX  <  ER3*  EO ,  -S I G3 /OMEGA ) 

WIR01730 

CC 

WRITE ( 6 , 7 755 ) ER3 , EO , EP3 

7755 

FORMAT (11, 'ER3 ,E0,EP3* ’ ,4E12 . 5) 

ETA»CDSQRT(U0/EP3) 

WIR01740 

CAM-OMEGA*CDSQRT ( -U0*EP3) 

WIR01750 

C*«COMPUTE  GAMMA  BY  EQUATION  IN  XAYT  PAGE  333 

CC 

GAM-JCOM*OMEGA»CDSQRT(UO»EP3)*CDSQRT(1.DO-JCOM*SIG3/'(JMEGA*EP3)) 

CC 

WRITE (6, 8888) GAM 

8888 

FORMAT (IX, ’KYAT  GAM* * ,2E12 . 5) 

AM-0.0001DC 

INT*0 

DZA00180 

X1-SX1 

DZA00190 

X2-SX2 

DZA00200 

X3-SX3 

DZA00210 

Y1*SY1 

DZA00220 

Y2-SY2 

D7.A00230 

Y3*SY3 

DZA00240 

Z1-SZ1 

DZA00250 

Z2-SZ2 

DZA 00260 

Z3*SZ3 

DZA00270 

IA-SXA 

T7A00280 

XB-SXB 

DZA 00290 

IC-SXC 

DZA00300 

YA-SYA 

DZA00310 

YB-SYB 

DZA00320 

YC*SYC 

DZA00330 

ZA-SZA 

DZA00340 

ZB-SZB 

DZA00350 

ZC*S7C 

DZA00360 

IBA-XB-XA 

DZA00370 

YBA-YB-YA 

DZA00380 

ZBA*Z8-ZA 

DZA00390 

X21-X2-X1 

DZA00400 

Y21-Y2-Y1 

DZA00410 

Z21-Z2-Z1 

DZA00420 

DS*DSQRT (XBA*XBA*YBA*YBA*ZBA*ZBA) 

DZA00430 

DT-DS0RT(X21*X21+Y21*Y21+Z21*Z21) 

DZA 00440 

CC 

DSK*B*OS 

DZA00450 

CC 

DTK*B*DT 

DZA00460 

cc 

CGDS*DCMPLX(DCOS(DSK) , 0 . ODO) 

DZA00470 

cc 

SGDS*DCMPLX (0 . ODO.DSIN (DSK) ) 

DZA00480 

cc 

SGDT-DCMPLX ( 0 . ODO , DS IN ( DTK ) ) 

DZA00490 

EGDS-CDEIP(GAM*DS) 

EGDT  *CDEXP ( G AM  *DT ) 

CGDS* (EGDS-H  . ODO/EGDS ) /2 . ODO 

SGDS*(EGDS- 1 . ODO/EGDS) /2 . ODO 

SCDT- ( EGDT- 1 . ODO/EGDT ) /2 . ODO 

CALX  DGGS ( XA , YA , ZA , XB , YB , ZB , X 1 , Y 1 , Z 1 , X2 , Y2 , Z2 , AM , 

DZA 00500 

2DS, CGDS, SGDS.DT.SGDT.INT, ETA, GAM, Pll ,P12 ,P21 ,P22) 

DZA00510 

CALL  DGGStXA , YA ,ZA , XB , YB , ZB , X2 , Y2 , Z2 , X3 , Y3 , Z3 , AM , 

DZA00S20 

2DS, CGDS, SGDS.DT.SGDT.INT, ETA, GAM, Qll ,Q12,Q21 ,Q22) 

DZA00530 
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CALL  DGGS(XB,YB,ZB,XC,YC,ZC,X1,Y1,Z1,X2,Y2,Z2,AM, 
2DS.CGDS .SGDS.DT.SGDT, INT,ETA,GaM,R11 .R12.R21 ,R22; 

CALL  DGGS(XB.YB,ZB,XC.YC  ZC ,X2 , Y2 ,Z2 ,X3 , Y3 ,Z3 , AM , 
2GS , CGLS , SGDS , DT,SGDT, INT ,£TA ,GAM,$11 , S12,S21 ,S22) 
DZABG«P22+Q21+R12+si'. 

CC  VRITE(6,8899)DZABG 

8899  FORMAT (IX, ’EXITING  DZABG,  DZATG- ■ .2E14 . 5) 

RETURN 

END 


DZA00S40 

DZAOOSSO 

DZA0GS60 

DZA00570 

DZA00580 


DZA00590 

DZA00600 
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